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CHAPTER  1 

WHAT  IS  A  SITE  SEISMICITY  STUDY? 


INTRODUCTION 

The  objective  of  a  seismicity  study  is  to  quantify  the  level  and  characteristics  of  ground 
motion  shaking  that  pose  a  risk  to  a  site  of  interest.  The  approach  taken  in  this  work  is  to  use 
the  historical  epicenter  data  base  in  conjunction  with  available  geologic  data  to  form  a  best 
estimate  of  the  probability  distribution  of  site  ground  motion. 

This  report  presents  techniques  that  have  been  automated  into  a  procedure  to  compute: 

•  Regional  earthquake  recurrence  parameters 

•  Regional  probability  distribution 

•  Fault  earthquake  magnitude  recurrence  parameters 

•  Probability  distribution  of  site  acceleration  from  each  fault 

•  Total  probability  distribution  of  site  acceleration 

•  Causative  magnitudes  and  separation  distances  associated  with  acceleration 

•  Response  spectra  based  on  site  soil  conditions  and  causative  events 


OUTLINE  OF  PROCEDURE 

The  procedures  were  developed  as  computer  programs  designed  to  run  on  standard 
desktop  DOS-based  computers.  System  requirements  include  640K  of  memory,  a  math 
coprocessor  chip,  and  a  hard  disk.  A  CDROM  is  required  to  use  the  recommended  epicenter 
data  base. 

The  procedure  consists  of: 

•  Evaluating  tectonics  and  geologic  settings 

•  Specifying  faulting  sources 

•  Determining  site  soil  conditions 

•  Determining  the  geologic  seismic  slip  rate  data 
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•  Specifying  the  epicenter  search  area  and  search  of  data  base 

•  Specifying  and  formulating  the  site  seismicity  model 

•  Developing  the  recurrence  model 

•  Determining  the  maximum  source  events 

•  Selecting  the  motion  attenuation  relationship 

•  Computing  individual  fault/source  seismic  contributions 

•  Summing  the  sources 

•  Determining  the  site  matched  spectra  for  causative  events 

This  report  will  present  a  summary  and  discussion  of  the  technology  for  each  of  the  elements  of 
the  analysis.  The  user’s  manuals  are  given  in  Appendixes  A  and  B  and  are  based  on  Ferritto 
(1993). 

APPLICATION 

The  procedures  were  developed  subject  to  the  following  limitations: 

•  The  exposure  period  or  life  of  the  structure  is  50  years. 

•  Return  times  of  events  of  interest  are  not  appreciably  longer  than  about  one  in  a 
thousand  years.  This  procedure  is  not  intended  to  predict  events  such  as  the 
10,000-year  event  with  high  accuracy. 


REFERENCE 

Ferritto,  J.M.  (1993).  Seismic  hazard  analysis,  Naval  Civil  Engineering  Laboratory,  User’s 
Guide  UG-0027.  Port  Hueneme,  CA,  1993. 
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CHAPTER  2 

EARTHQUAKE  ENGINEERING  FUNDAMENTALS 


PLATE  TECTONICS 

The  United  States  is  located  on  the  North  American  plate,  the  western  portion  of  which 
meets  the  Pacific  plate.  The  interaction  of  these  two  plates  is  responsible  for  the  high  seismic 
activity  that  has  in  the  past  and  continues  now  to  take  place  in  the  Western  United  States.  Plate 
tectonic  theory  has  explained  much  of  the  geologic  activity.  Also  explained  in  terms  of  plate 
tectonic  theory  is  the  seismic  activity  experienced  in  the  Central  and  Eastern  United  States.  This 
midplate  activity  can  be  very  destructive. 

Figure  2-1  shows  a  cross  section  of  the  earth.  The  lithosphere,  composed  largely  of 
basalt,  extends  to  an  average  depth  of  about  100  km.  Below  the  lithosphere  is  the  asthenosphere 
that  extends  to  a  depth  of  400  km.  Because  its  upper  portion  is  partially  molten,  seismic  velocity 
in  that  region  is  decreased.  The  lithospheric  plates  are  able  to  "float"  on  this  plastic  layer.  Not 
all  of  the  asthenosphere  is  molten;  however,  there  is  a  rigid  portion. 

Most  seismic  activity  is  located  at  plate  boundaries  and,  therefore,  boundaries  are  of 
considerable  interest.  There  are  three  major  types  of  interaction  between  adjacent  plates:  (1) 
spreading  boundaries,  (2)  converging  boundaries,  and  (3)  transforming  boundaries.  Figure  2-2 
illustrates  the  three  kinds  of  boundaries. 


GEOLOGIC  FAULTS  AND  EARTHQUAKES 

Since  the  San  Francisco  earthquake  of  1906  and  the  subsequent  work  on  the  elastic 
rebound  theory  of  earthquakes,  general  agreement  has  been  reached  on  the  close  relationship 
between  earthquakes  and  geologic  faults.  Most  tectonic  earthquakes  that  cause  major  structural 
damage  are  associated  with  fracture  on  a  fault.  Plate  motion  causes  stress  in  the  earth’s  rock 
crust.  Earthquakes  occur  when  the  strength  of  the  fault  can  no  longer  withstand  the  stress  that 
has  built  up.  Fault  plane  solutions  and  earthquake  mechanism  studies  have  contributed  to  a 
consistent  picture  of  the  earthquake  generation  process  that  satisfactorily  explains  most  of  the 
observed  facts. 

A  fault  is  a  rupture  in  the  earth  along  which  opposite  faces  have  been  displaced.  The 
basic  kinds  of  faults  are  illustrated  in  Figure  2-3  and  are  defined  as  follows: 

1 .  Strike  Slip.  Strike  is  the  direction  along  a  fault,  and  strike  slip  refers  to  displacement 
along  this  line.  Right  lateral  or  left  lateral  refers  to  the  direction  of  movement  of  the  opposite 
side  when  one  faces  the  fault. 

2.  Normal.  A  normal  fault  refers  to  movement  of  one  side  of  the  fault  away  from  the 
other,  producing  tension. 
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3.  Thrust.  A  thrust  fault  refers  to  movement  of  one  side  toward  or  over  the  other  side 
producing  compression. 

Most  faults  are  combinations  of  strike  slip  and  normal  or  thrust  movement.  The  fault 
plane  itself  can  be  curved  and  blocks  can  be  rotated  relative  to  each  other.  The  fault  trace  is  the 
line  of  the  fault  along  the  ground  surface.  The  strike  of  a  fault  is  measured  from  north  in 
degrees.  The  dip  of  a  fault  is  used  to  measure  the  slope  of  the  fault  plane  with  the  surface 
(Figure  2-4). 

Some  misunderstandings  have  occurred,  and  perhaps  some  significant  differences  of 
opinion,  about  the  direct  relationship  between  geologic  faults  and  the  earthquake  hazard.  There 
is  evidence  dating  from  1906  to  suggest  that  destructive  ground  shaking  is  not  necessarily  at  a 
maximum  in  the  immediate  vicinity  of  the  causative  fault  (Hudson,  1972).  More  often  than  not, 
the  maximum  destructive  ground  shaking  is  miles  from  the  fault,  as  explained  by  a  number  of 
the  features  of  the  generation  and  propagation  of  seismic  waves.  Classical  photographs  of  the 
1906  major  movements  along  the  San  Andreas  fault,  for  example,  show  horizontal  surface 
displacements  of  as  much  as  S  meters  passing  close  to  a  small  wood-frame  house  that  received 
no  significant  damage.  Similarly,  during  the  San  Fernando  earthquake  of  9  February  1971,  a 
2-meter  vertical  fault  scarp  passed  directly  through  a  wooden  barn  just  a  short  distance  from  a 
single-story  residence.  The  bam  was  severely  damaged,  but  no  significant  structural  damage  to 
the  house  was  noted  (Hudson,  1972).  The  San  Fernando  earthquake  also  furnished  numerous 
examples  of  surface  faulting  passing  through  heavily  populated  areas.  Although  severe  structural 
deformation,  with  a  resulting  economic  loss,  occurred  in  numerous  cases,  catastrophic  collapses 
leading  to  loss  of  life  and  serious  injury  were  not  directly  associated  with  these  surface  breaks. 
Hazardous  collapses  were  in  all  cases  the  consequence  of  severe  ground  shaking,  which  is 
pervasive  over  a  large  area  and  is  not  limited  to  the  vicinity  of  faults. 

The  focus  or  hypocenter  is  the  point  within  the  earth’s  crust  where  the  initial  rupture 
occurs  and  from  which  the  first  waves  are  released.  The  projection  of  this  point  to  the  ground 
surface  is  the  epicenter.  The  epicenter  and  hypocenter  do  not  necessarily  indicate  the  center  of 
the  total  energy  release  of  the  earthquake,  but  rather  the  point  where  the  seismic  energy  waves 
were  first  created.  For  small  earthquakes,  the  center  of  total  energy  release  and  the  epicenter 
are  not  far  apart  because  the  fault  break  length  is  short;  however,  this  is  not  the  case  for  large 
earthquakes.  The  majority  of  earthquakes  in  the  United  States  have  had  relatively  shallow  focal 
depths  (0  to  40  km).  In  California,  earthquakes  have  occurred  in  regions  where  surface  fault 
patterns  were  clearly  visible.  In  the  Puget  Sound  area,  earthquakes  are  focused  at  deeper 
locations  within  the  earth’s  crust  so  that  a  surface  rupture  is  not  observable.  In  the  eastern 
United  States,  the  relationship  to  surface  rupture  in  general  has  not  been  closely  identified  (Bolt, 
1970;  Newmark  and  Rosenblueth,  1971). 

A  fault  undergoing  tectonic  creep,  or  one  with  abrupt  displacement,  causes  changes  in  the 
terrain  it  crosses.  Very  distinctive  patterns  are  produced  where  active  faults  cross  streams,  such 
as  landslides.  The  ongoing  geologic  process  causes  scarps,  trenches,  sag  ponds,  and  stream 
offsets.  Figure  2-5  shows  a  landform  with  an  active  fault  (Wesson,  et  al.,  1975). 

Estimates  of  the  maximum  size  and  frequency  of  earthquakes  on  a  fault  are  based  on  the 
geologically  determined  slip  rate  and  the  historic  record  of  ground  deformation  (where  available), 
the  seismic  history  of  the  fault  and  surrounding  tectonic  region,  a  geological  evaluation  of  the 
tectonic  setting,  and  empirically  derived  relationships  between  earthquake  magnitude  and  fault 
length. 
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SURFACE  EFFECTS  OF  FAULT  MOVEMENTS 


Wher.  faults  are  considered,  the  assumption  is  commonly  made  that  the  creation  of 
entirely  »cw  faults  by  an  earthquake  is  unlikely  (Krinitzsky,  1974).  Significant  surface  faults 
and  their  activity  can  be  found  by  proper  geologic  investigations.  Cluff,  Slemmons,  and 
Waggoner  (1970)  have  studied  the  character  of  typical  surface  effects  by  faults.  These  are 
illustrated  in  Figures  2-6  through  2-8,  as  reported  by  Krinitzsky  (1974). 


EARTHQUAKE  MAGNITUDE 

Engineers  can  define  a  design  earthquake  for  a  site  in  terms  of  the  earthquake  magnitude, 
M,  and  the  strength  of  ground  motion.  Factors  that  influence  the  selection  of  a  design 
earthquake  are  the  length  of  geologic  fault  structures,  the  relationship  of  the  fault  to  the  regional 
tectonic  structure,  the  geologic  history  of  displacement  along  the  structure,  and  the  seismic 
history  of  the  region. 

The  design  earthquake  in  engineering  terms  is  a  specification  of  levels  of  ground  motion 
that  the  project  is  required  to  survive  successfully  with  no  loss  of  life  and  acceptable  damage  and 
loss  of  service.  A  design  earthquake  on  a  statistical  basis  considers  the  probability  of  the 
recurrence  of  a  historical  event. 

Earthquake  magnitudes  can  be  specified  in  terms  of  a  design  level  earthquake  that  can 
reasonably  be  expected  to  occur  during  the  life  of  the  structure.  As  such,  this  represents  a 
service  load  that  the  structure  must  withstand  without  significant  structural  damage  or  interruption 
of  a  required  operation.  A  second  level  of  earthquake  magnitude  is  a  maximum  credible  event 
for  which  the  structure  must  not  collapse;  however,  significant  structural  damage  can  occur.  The 
inelastic  behavior  of  the  structure  must  be  limited  to  ensure  the  prevention  of  collapse  and 
catastrophic  loss  of  life. 

The  selection  of  a  magnitude  level  may  be  based  on: 

1 .  Known  design-level  and  maximum-credible  earthquake  magnitudes  associated  with 
a  fault  whose  seismicity  has  been  estimated. 

2.  Specification  of  probability  of  occurrence  for  a  given  life  of  the  structure  (such 
as  having  a  10  percent  chance  of  being  exceeded  in  25  years). 

3.  Specification  of  required  level  of  ground  motion  as  in  a  code  provision. 

4.  Fault  length  empirical  relationships. 

Earthquake  magnitude  can  be  related  to  length  of  fault  for  shallow  depth  earthquakes.  Data  have 
been  plotted  by  Seed,  et  al.  (1969),  Krinitzsky  (1974),  Housner  (1965),  and  Tocher  (1958)  to 
provide  the  curves  indicated  (Figure  2-9).  It  is  important  to  note  that  in  some  regions, 
correlations  of  these  types  are  of  little  value  since  many  of  the  important  geologic  features  can 
be  deeply  buried  by  weathered  materials. 

Magnitude  as  measured  on  the  Richter  scale  is  calculated  from  a  standard  earthquake,  one 
which  provides  a  maximum  trace  amplitude  of  1  Mm  on  a  standard  Wood-Anderson  torsion 
seismograph  at  a  distance  of  100  km.  Magnitude  is  the  logt0  of  the  ratio  of  the  amplitude  of  any 
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earthquake  at  the  standard  distance  to  that  of  the  standard  earthquake.  Each  full  numeral  step 
in  the  scale  (two  to  three,  for  example)  represents  an  energy  increase  of  about  32  times. 
Experience  with  past  earthquakes  is  presently  the  only  useful  basis  for  relating  fault  length  and 
motion  to  magnitudes  of  associated  earthquakes. 


FAULT  LENGTH  AND  EARTHQUAKE  MAGNITUDE 

A  useful  insight  into  the  relationship  between  earthquake  magnitude  and  length  of 
observed  fault  slippage  is  presented  by  Iida  (1965).  He  groups  faults  of  all  types  and  shows  a 
wide  spread  of  points  for  which  upper  and  lower  boundaries  are  drawn,  lida’s  wide  spread  of 
values  should  be  kept  in  mind  when  one  considers  the  linear  relationships  that  have  been 
suggested  by  numerous  authors. 

Fault  movements  below  a  magnitude  of  5  are  usually  contained  in  the  subsurface.  If 
active  fault  movement  is  found  at  a  site,  even  short  movement,  it  should  be  viewed  as  evidence 
for  an  earthquake  capability  of  greater  than  5  (Krinitzsky,  1974). 

It  is  important  to  note  the  spread  in  the  data.  Krinitzsky  (1974)  concludes,  "Any  fault 
break  in  competent  rock,  no  matter  how  small,  should  be  taken  as  indicative  of  the  capability  for 
at  least  a  5.4  earthquake. "  It  is  important  that  the  local  seismic  history  and  the  behavior  of  other 
analogous  faults  be  considered. 

Consideration  should  be  given  to  the  possibility  of  not  identifying  all  the  faults  in  a  region 
that  may  be  active.  This  is  especially  true  in  the  Central  and  Eastern  United  States.  To  account 
for  this  a  "floating  earthquake"  (one  that  may  be  assumed  capable  of  occurring  anywhere  in  the 
region)  should  be  considered  (Krinitzsky,  1974). 

Mark  and  Bonilla  (1977)  evaluated  data  to  develop  relationships  between  surface  fault 
displacement  and  earthquake  magnitude.  More  recent  data  by  Coppersmith  (1991)  are  shown 
in  Figures  2-10  and  2-11. 


GEOLOGICALLY  DETERMINED  SLIP  RATES 

The  offset  of  distinctive  rock  units  establishes  the  rate  of  fault  movement  within  fairly 
wide  bounds.  Commonly  these  offsets  average  the  rate  of  movement  over  millions  of  years,  and 
sudden  slip  cannot  be  distinguished  from  creep.  Data  for  the  San  Andreas  fault  suggest  an 
average  slip  rate  of  1  to  2  cm/yr  over  the  last  20  million  years.  But  to  predict  movements  in  the 
immediate  future,  the  most  recent  hundreds  to  thousands  of  years  are  the  most  important 
(Wesson,  et  al.,  1975). 

The  history  and  rate  of  fault  movement  have  been  obtained  within  this  brief  time  period 
in  a  few  special  circumstances  in  southern  California,  using  absolute  age  dating  techniques. 

Figure  2-12  is  a  simplified  sketch  of  a  trench  wall  showing  vertical  deformation  of 
initially  flat-laying  sediments  and  sedimentary  contacts  associated  with  predominantly  horizontal 
movement  on  the  Coyote  Creek  fault  of  southern  California  (from  Clark,  et  al.,  1972).  The 
trench,  dug  shortly  after  the  1968  Borrego  Mountain  earthquake,  crosses  a  branching  break  of 
the  fault  zone  along  which  about  50  mm  of  vertical  displacement  and  about  the  same  amount  of 
horizontal  displacement  took  place  during  the  earthquake.  Deposits  at  points  A,  B,  and  C  were 
dated  radiometrically.  The  vertical  displacement  of  the  sedimentary  contacts  plotted  against  the 
age  of  the  corresponding  deposits  yields  an  average  rate  of  vertical  deformation  of  about  0.5 
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mm/yr  for  the  past  3,000  years.  This  suggests  a  recurrence  interval  for  earthquakes  the  size  of 
the  1968  event  of  about  200  years,  Figure  2-13  (Clark,  et  al.,  1972). 


RECURRENCE  INTERVALS  FROM  GEOLOGIC  SLIP  RATES 

Wallace  (1970)  presented  an  approach  that  has  been  used  by  Lamar,  et  al.  (1973).  The 
recurrence  interval  at  a  point  can  be  estimated  by: 


R 


X 


(2-1) 


where: 


D 

S 

C 

S  -  C 


=  displacement  per  seismic  event 
=  long-term  slip  rate 
=  tectonic  creep  rate 
=  average  seismic  slip  rate 


Lamar,  et  al.  (1973)  present  the  following  quoted  discussion: 

"The  following  assumptions  are  made:  (1)  Slip  on  faults  occurs  incrementally 
during  earthquakes  and  will  continue  at  the  same  rate  as  that  determined  from 
geodetic  data  and  offset  of  geologic  units.  (2)  Elastic  strain  accumulates  between 
earthquakes;  the  displacement  during  an  earthquake  represents  the  release  of  this 
accumulated  elastic  strain.  (3)  Tectonic  creep  is  aseismic  slip  which  reduces  the 
accumulation  of  elastic  strain  available  for  release  during  earthquakes. 
...Recurrence  intervals  determined  by  [Equation  2-1]  represent  a  long-term 
average;  there  is  however,  evidence  of  significant  local  (Ambraseys,  1970)  and 
worldwide  (Davies  and  Brune,  1971)  time  variations  in  the  level  of  seismic 
activity." 


For  most  faults,  creep  can  be  evaluated.  Therefore,  as  an  expedient,  Equation  2-1  is 
simplified  as: 


R 


X 


D 

S 


(2-2) 


Equation  2-2  is  appropriate  when  the  rupture  length  is  large  compared  to  the  distance  of  the  site 
to  the  fault.  When  this  is  not  the  case,  Equation  2-2  is  multiplied  by  the  ratio  of  length  of 
rupture  to  total  fault  length  to  account  for  the  spatial  distribution. 
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APPLICATION  OF  SLIP  RATE  TO  COMPUTE  RECURRENCE  DATA 


Lamar,  et  al.  (1973)  have  investigated  the  occurrence  data  for  six  faults  in  the  southern 
California  area.  Their  data  are  summarized  in  Table  2-1  and  illustrate  the  concept.  The 
displacement,  D,  used  in  Equation  2-2  to  calculate  recurrence  intervals  at  a  given  point  on  a  fault 
were  derived  from  Bonilla  and  Buchanan  (1970),  except  in  the  case  of  the  White  Wolf  and  Sierra 
Madre  faults.  For  these  two  faults  historic  displacements  were  available,  and  they  deviated 
significantly  from  the  least-squares-fit  curve  for  reverse  faults. 

Quoting  from  Lamar,  et  al.  (1973): 

"The  determination  of  long-term  slip  rates  and  recurrence  intervals  provides  a  new 
approach  to  earthquake  hazard  assessment.  The  results  can  be  strikingly  different 
from  those  based  on  the  historic  earthquake  record.  For  example,  the  ...  Garlock 
fault,  which  has  not  caused  damaging  historic  movement,  may  have  accumulated 
sufficient  elastic  strain  for  a  destructive  earthquake.  On  the  other  hand,  historic 
ruptures  on  faults  such  as  the ...  San  Fernando,  have  released  [some]  accumulated 
elastic  strain,  so  that  a  destructive  earthquake  [may  be  less  probable]  for  the  next 
few  hundred  years  [depending  on  the  amount  of  strain  release  and  buildup].  The 
recurrence  intervals  in  [Table  2-1]  must  be  considered  tentative  and  are  subject  to 
revision  as  new  information  becomes  available.  For  the  most  part,  slip  rates  are 
poorly  known,  and  the  curves  relating  magnitude  to  surface  displacement  and 
rupture  length  are  based  on  meager  data  with  considerable  scatter.  More  accurate 
age-dates  of  offset  strata  on  faults  are  needed,  and  additional  studies  following 
earthquakes  throughout  the  world  are  required  to  refine  the  empirical  relations 
between  magnitude,  surface  displacement  and  rupture  length.  This  research  offers 
the  prospect  of  more  quantitative  assessments  of  earthquake  risk." 

Since  the  early  work  of  Wallace  (1970),  more  emphasis  has  been  placed  on  use  of 
geologic  data.  The  historic  seismicity  record  in  the  United  States  and  other  areas  is  generally 
too  short  to  fully  define  the  recurrence  of  particular  individual  faults  for  low  probability  events. 
Fault  slip  rates  derived  from  geologically  defined  intervals  afford  the  opportunity  of  spanning 
several  cycles  of  large  earthquakes  on  a  fault.  Coppersmith  and  Youngs  (1990)  note  that  the  best 
geologic  units  for  assessing  slip  rate  for  recurrence  purposes  are  late-Quartemary  or  Holocene 
units.  Assessing  slip  rates  over  relatively  young  units  will  avoid  averaging  out  long-term  changes 
in  the  slip  rate  from  regional  changes  in  tectonic  stress. 


SEISMIC  MOMENT 

Seismic  moment  has  been  used  in  conjunction  with  slip  rate.  Seismic  moment,  M0,  is 
a  means  of  describing  the  size  of  an  earthquake  in  terms  of  physical  parameters: 


M0  =  |iAfD 


(2-3) 
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where:  n  =  rigidity  or  shear  modulus  (taken  as  3  x  1011  dyne/cm2) 

Ar  =  rupture  area  on  the  fault  plane 
D  =  average  displacement  over  slip  surface 

The  total  seismic  moment  rate  can  be  estimated  using  the  above  formulation  substituting  the  total 
fault  plane  area  and  the  average  slip  rate  along  the  fault  instead  of  the  displacement.  Thus,  the 
seismic  moment  rate  provides  a  link  between  geologic  data  and  seismicity.  Seismic  moment  rates 
determined  from  slip  data  can  be  compared  with  seismic  moment  rates  based  on  seismicity  data 
(Youngs  and  Coppersmith,  1985). 

Seismic  moment,  M0,  can  be  related  to  magnitude,  m,  as  follows: 


log  Me  =  C  m  +  d 


(2-4) 


Hanks  and  Kanamori  (1979)  report  that  c  =  1.5  and  d  =  16.1  in  California.  The  moment 
magnitude,  m,  is  considered  equivalent  to  local  magnitude  when  in  the  range  of  3  to  7  and  to 
surface  wave  magnitude  when  in  the  range  of  5  to  7.5. 
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Table  2-1 

Recurrence  Intervals  on  Selected  Faults  in  Southern  California 
(from  Lamar,  et  al.,  1973,  edited  by  Southern  California  Section 
of  the  Association  of  Engineering  Geologists) 


Fault 

Slip 

Rate 

(cm/yr) 

Length 

(km) 

Recurrence 

Intervale  (yr)  al 
a  Point  on  Fault 
<*,> 

Recurrence 

Intervale  (yr)  over 

Length  of  Fault 
<*0 

M  6 

M  7 

M  8 

M  6 

M  7 

M  8 

Noithweat  Trend  Right-Slip 

■ 

San  And  real  (aouthem  segment) 

500 

10 

40 

200* 

0.3-1 

3-10 

40-100* 

San  Jacinto  Fault  System 

B 

440 

too* 

400* 

2000 

4-10* 

40-100* 

400-1000 

Whittier-  Elsinore  -  Ague  Calienle- 
Laguna  Salads 

0.08 

260 

300 

2000 

6000 

20-90 

200-900 

3000-6000 

Northeast  Trend  Left-Slip 

Big  Pine 

0.2 

95 

100 

600* 

3000 

20-100 

300-600* 

3000 

Garlock 

0.8 

250 

30 

200 

600 

2-10 

30-90 

300600 

■ 

White  Wolf 

0.04 

53 

1000 

2000* 

4000 

200-900 

1000-2000* 

4000 

Sierra  Mad  re  Fault  System 

m 

90 

100* 

300* 

800 

30-100* 

100-300* 

800 

•Most  likely,  bated  on  historic  record. 
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CONTINENTAL  CRUST 


BASALTIC 
OCEAN  CRUST « 


Figure  2-1.  Cross  section  of  earth  structure  (from  Walper,  1976). 


TRANSFORM 
FAULT 


a.  Divergent  or  spreading 
boundary 


b.  Convergent  or  collision 
boundary 


c.  Transform  or  transcursion 
boundary 


Figure  2-2.  The  three  kinds  of  boundaries  between  plates 
(from  Walper,  1976). 
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b  NORMAL  FAULT 


a  STRIKE-SLIP  FAULT 
A*  =  STftlKC  SLi* 


*•  a  $L‘* 
ac  = thro* 

•C  «  mCAvC 


f 


C  THRUST 'REVERSE'  FAULT 

a«  «  »Lt# 

ac  »  throw 

•C  *  HIAVC 


d.  ROTATIONAL  NORMAL  FAULT  WITH  STRIKE 
SLIP  COMPONENT 

A*  a  NtT  SL*»  AO  «  TMItOW 

AK  OR  Cl  «  STRtKl  su»  OC  *  MCAvC 
AC  OH  (•  9  d\+  1L«A 


Figure  2-3.  Fault  terminology  (from  Krinitzsky,  1974). 


Figure  2-4.  Fault  measurement  (from  Krinitzsky,  1974). 
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Figure  2-5  Landforms  developed  along  recently  active  strike-slip 
faults  (after  U.S.  Geological  Survey,  1975). 


15 


Figure  1-6.  Damage  associated  with  movement  along  strike-slip 
fault  (from  Krinitzsky,  1974). 


b.  DURING 


e.  AFTER 


Figure  2-7.  Damage  from  normal  fault 
(from  Krinitzsky,  1974). 


b.  DURING 


c.  AFTER 


Figure  2*8.  Damage  from  thrust  (reverse)  fault 
(from  Krinitzsky,  1974). 


LENGTH,  KlUJkiTERS 


0.S  1  S  10  SO  100  soo 

LENGTH  (MILES)  OF  SURFACE  RUPTURE.  MAIN  FAULT 


Figure  2-9.  Earthquake  magnitude  versus  length  of  surface  rupture 

(from  Krinitzsky,  1974). 
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Average  Displacement  (m) 


Figure  2-10.  Earthquake  magnitude  versus  average  surface  displacement, 
(from  Coppersmith,  Proceedings  Fourth  International 
Conference  on  Seismic  Zonation,  1991) 
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Surface  Rupture  Length  (km) 


Figure  2-11.  Earthquake  magnitude  versus  fault  surface  rupture  length, 
(from  Coppersmith,  Proceedings  Fourth  International 
Conference  on  Seismic  Zonation,  1991) 
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DEPTH,  IN  METRES 


EAST 


HORIZONTAL  DISTANCE.  IN  METRES 

Figure  2-12.  Sketch  of  trench  wall  (from  Clark  et  al.,  1972). 


205y«ars 

Figure  2-13.  Estimated  recurrence  data  (from  Clark  et  al.,  1972). 
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OtPTM.  IN  FC€T 


CHAPTER  3 

EPICENTER  DATA  BASE 


INTRODUCTION 

The  U.S.  Geological  Survey  National  Earthquake  Information  Center,  Denver,  Colorado 
has  produced  "The  Global  Hypocenter  Data  Base"  CDROM  which  contains  parameters  for  more 
than  438,000  earthquake  events.  Seven  world- wide  and  12  regional  earthquake  catalogs  were 
assembled  to  produce  this  data  base,  spanning  a  time  period  from  2100  B.C.  through  1990. 
Useful  data  for  the  United  States  is  generally  constrained  to  the  period  when  instruments  were 
available  to  compute  event  magnitude.  Each  earthquake  is  detailed  where  data  are  available  with 
date,  origin  time,  location,  4  magnitude  estimates,  intensity,  and  cultural  effects. 


DATA  BASE  SEARCHING 

A  computer  program,  EPIC,  is  available  for  searching  the  CDROM.  EPIC  makes  data 
available  to  information  users  via  a  user-defined  search  request.  The  request  determines  which 
steps  are  necessary  to  produce  the  desired  output.  It  includes  a  search  method,  combination  of 
search  elements,  and  a  destination  output  format,  including: 

•  Global 

•  Rectangular  grid 

•  Circle 

•  Irregular  grid 
Search  elements  include: 

•  Date  interval  (year/month) 

•  Magnitude  range 

•  Catalog  selection 

•  Depth  interval 

•  Intensity  range 

•  Magnitude  range  or  intensity  interval 


22 


•  Duplicate  earthquake  eliminator 

•  Cultural  effects  (casual  i  ties/ damage) 

•  Information  (fault  plane/moment  tensor  solutions) 

•  Associated  seismic  phenomena  (e.g.,  faulting,  ground  response) 

An  automated  plotting  package  that  produces  seismicity  maps  in  multicolor  or 
monochrome  is  incorporated  into  the  EPIC  software.  The  data  to  be  mapped  are  extracted  from 
the  selected  data  and  plotted  in  a  global  or  regional  format.  The  map  may  be  delayed  or 
cancelled  at  user’s  discretion. 

The  availability  of  the  CDROM  data  base  of  epicenters  and  EPIC  software  greatly 
facilitates  creation  of  the  historical  epicenter  subset  required  for  use  of  the  automated  seismicity 
analysis  tools  developed  and  presented  in  this  report.  Any  data  base  can  be  used  provided  the 
data  are  in  the  same  format  as  the  CDROM  and  EPIC  data.  Details  are  presented  in  the  EPIC 
user’s  manual,  which  will  not  be  repeated  here. 


DATA  BASE  DISCUSSION 

A  number  of  data  fields  for  some  events  are  unfilled  because  the  information  is  not 
available.  Information  on  cultural  effects,  intensity,  and  other  phenomena  associated  with  the 
event  has  been  included  for  earthquakes  in  the  United  States.  This  information  has  sometimes 
been  entered  for  non-United  States  earthquakes,  particularly  since  May  1968,  although  significant 
gaps  still  exist. 

The  quality  of  epicenter  determinations  varies  significantly  with  the  time  period  studied. 
Before  1900,  locations  are  usually  noninstrumentally  determined  and  are  given  as  the  center  of 
the  macroseismic  effects.  Most  instrumental  epicenters  prior  to  1961,  excluding  local 
earthquakes  in  California,  were  located  to  the  nearest  1/4  or  1/2  degree  of  latitude  and  longitude. 
Reliable  information  on  the  quality  of  many  epicenter  determinations  is  lacking.  Beginning  in 
1960,  epicenters  have  been  determined  by  computer,  and  the  accuracy  is  generally  better. 
However,  although  stated  to  tenths  or  hundredths  of  a  degree,  the  location  accuracy  is  usually 
a  few  tenths  of  a  degree.  Since  May  1968,  the  latitude  and  longitude  values  for  most  events 
have  been  listed  to  three  decimal  places.  This  precision  is  not  intended  to  reflect  the  accuracy 
of  the  location  of  events  except  for  local  California  earthquakes  and  special  epicenter 
determinations.  Where  several  sources  have  determined  an  epicenter  for  the  same  earthquake, 
one  solution  has  been  designated  as  the  most  reliable.  Usually  it  is  the  source  believed  to  contain 
the  best  data  set  for  the  earthquake.  In  some  cases,  data  from  two  sources  were  combined  to 
provide  a  more  complete  record. 

Magnitudes  from  a  number  of  different  sources  are  included  in  the  earthquake  data  file. 
Gutenberg  and  Richter  (1954)  and  Richter  (1958)  discuss  the  development  of  the  magnitude 
scale.  Many  magnitudes  published  by  Gutenberg  and  Richter  (1954)  were  later  revised  by 
Richter  (1958).  The  revised  magnitudes  are  used  in  the  file  even  though  the  source  is  identified 
as  Gutenberg  and  Richter  (1954).  The  concept  of  earthquake  magnitude  is  not  restricted  to  one 
value.  Several  definitions  are  possible,  depending  on  which  seismic  waves  are  measured.  Three 
different  magnitude  scales,  BODY  WAVE  (MB),  SURFACE  WAVE  (MS),  and  LOCAL  (ML), 
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are  distinguished  in  this  file.  In  addition,  another  data  field,  OTHER  MAGNITUDE,  was 
included  when  it  was  unclear  which  scale  was  used.  Richter  (1958)  and  other  modem 
seismology  references  provide  detailed  discussions  of  this  topic.  The  different  scales  do  not  give 
exactly  comparable  results,  and  different  values  frequently  are  given  for  the  same  earthquake. 
It  is  common  practice  to  average  the  individual  magnitudes  from  different  stations  to  get  a  more 
uniform  value  within  each  scale  (MB,  MS,  and  to  a  lesser  extent  ML). 

In  general,  the  file  contains  earthquakes  of  magnitude  4.0  or  less  only  for  the  United 
States  region  and  for  areas  within  dense  seismic  station  networks.  However,  no  claim  is  made 
for  the  statistical  homogeneity  of  these  events.  Inclusion  of  earthquakes  of  magnitude  4.0  to  5.0 
also  is  influenced  by  the  proximity  of  seismic  stations  to  the  source  or  epicenter. 

A  maximum  intensity  is  listed  for  many  of  the  earthquakes.  Each  is  assigned  according 
to  the  Modified  Mercalli  Intensity  Scale  of  1931.  Some  of  these  values  have  been  converted 
from  reported  intensities  on  other  scales. 


CORRECTION  OF  EPICENTERS 

Early  records  of  earthquakes  may  show  inaccurate  epicentral  locations.  Also,  magnitudes 
or  intensities  may  differ  from  values  that  would  now  be  assigned.  As  an  example,  Krinitzsky 
(1974)  reports  the  following: 

"The  location  of  a  listed  earthquake  of  15  May  1909  in  southern 
Saskatchewan  is  shown  in  [Figure  3-1).  It  has  a  maximum  intensity  of  VIII  and 
is  the  largest  recorded  event  in  this  general  area.  Its  location  was  made  when 
seismograph  stations  were  few,  and  those  that  operated  were  far  less  accurate  than 
they  are  today.  The  intensity  of  VIII  was  assigned  on  the  basis  of  the  large  felt 
area.  The  event  was  checked  by  referring  to  the  newspaper  accounts  of  this  time. 

This  was  not  too  formidable  a  task.  State  archives,  state  university  collections, 
and  national  archives  often  have  collections  of  local  newspapers  gathered  in 
central  depositories.  One  may  request  microfilm  copies  of  these  papers  for  the 
dates  of  interest,  review  them,  and  assign  intensity  values  based  on  their 
descriptions.  This  exercise  provided  intensities  for  more  than  50  communities 
though  this  was,  and  still  is,  a  thinly  populated  area.  The  resulting  iso-intensity 
map  shows  there  was  no  intensity  VIII.  The  greatest  was  VI.  Also,  the  region 
of  VI  had  its  center  to  the  east  by  about  1  degree  from  the  NOAA  location.  In 
[Figure  3-1],  the  revised  location  is  shown  in  relation  to  three  other  earthquakes 
of  1956,  1968,  and  1972.  Their  locations  were  accurate  to  begin  with  because 
of  better  instrumental  capability.  One  is  associated  with  a  seismically  interpreted 
fault  that  also  agrees  with  a  geologically  mapped  fault.  Its  trend  is  toward  the 
three  other  events.  Thus,  the  revised  location  for  the  1909  event  helps  to  interpret 
a  fault  trend." 


SEISMIC  ARRAYS 

Seismometers  have  been  installed  near  known  active  faults  to  record  microearthquakes. 
Figures  3-2  and  3-3  show  the  location  of  recorded  microearthquakes  for  a  year,  and  major  fault 
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zones  in  California.  The  events  recorded  range  in  magnitude  from  0.5  to  1 .5.  The  figures  show 
areas  where  the  microearthquakes  closely  trace  the  faults.  There  are  also  regions  where  few 
microearthquakes  occurred.  The  San  Andreas  fault  north  of  the  Sargent  fault  exhibits  less 
activity  than  southern  portions,  perhaps  indicating  it  is  locked.  However,  sufficient 
microearthquakes  have  occurred  to  show  the  continuation  of  the  fault.  It  may  be  reasoned  that 
small  failures  might  occur  before  a  major  rupture  occurs;  alternatively,  the  large  number  of 
microearthquakes  demonstrate  active  creep  which  may  be  sufficient  to  prevent  sizable  strain 
accumulation  and  preclude  a  large  event. 

Since  earthquakes  are  associated  with  faults,  it  might  be  thought  that  they  should  precisely 
overlay  the  fault  location.  This  is  not  the  case  because  the  distribution  of  seismometers  is  uneven 
and  has  changed  with  time.  There  are  limitations  in  the  accuracy  of  the  techniques  used  to  locate 
epicenters,  principally  from  variations  in  the  assumed  propagation  velocities.  Further  explanation 
for  the  location  of  epicenters  being  off  their  associated  fault  comes  from  the  simplified  model 
used  to  locate  them.  This  is  explained  in  Figure  3-4.  Note  that  the  center  of  earthquake  energy 
is  located  at  the  focus.  For  an  inclined  fault,  the  surface  location  (epicenter)  is  a  distance 
removed  from  the  surface  fault  location.  It  is  only  in  vertical  faults  that  one  might  expect  the 
epicenter  to  lay  on  the  fault. 

Krinitzsky  (1974)  concludes  that  earthquakes  can  be  related  to  existing  faults  and  that  the 
possibility  of  formation  of  new  faults  should  not  be  considered  in  design.  Large  earthqu  s 
require  fault  breaks  of  considerable  distance.  The  uncertainties  that  occur  in  the  association  of 
earthquakes  with  faults  can  occur  only  for  small  events.  Generally,  in  the  Western  United  States 
the  extent  of  geologic  investigation  precludes  the  omission  of  a  large  fault  remaining  unknown. 
However,  there  are  uncertainties  associated  with  eastern  earthquakes.  For  example,  causative 
faults  responsible  for  the  New  Madrid  earthquake  of  181 1  and  1812  have  not  yet  been  identified. 
This  may  be  the  result  of  insufficient  geologic  investigations.  The  importance  of  considering  the 
extent  and  quality  of  geologic  investigations  is  evident. 


LIMITATIONS  TO  HISTORIC  DATA 

A  period  of  demonstrated  quiescence  over  a  geological  time  period  indicates  inactivity  of 
the  fault  and  probable  continued  inactivity.  However,  inactivity  over  a  period  of  historic 
recording  (50  to  100  years)  does  not  imply  future  inactiviy.  Rather,  it  may  point  to  a  region 
which  $s  locked  and  through  which  a  major  fault  rupture  may  propagate.  A  number  of 
earthquakes  producing  damage  in  southern  California  occurred  on  faults  lacking  historic  activity. 
Caution  must  be  exercised  to  recognize  that  the  accuracy  of  an  incomplete  data  base  is  very 
limited  when  extrapolated  for  return  periods  greatly  exceeding  the  length  of  the  period  of 
recorded  data.  Furthermore,  aftershocks  must  be  distinguished  from  main  shocks.  An  area 
having  recently  undergone  a  large  event  releasing  strain  built  up  for  hundreds  or  thousands  of 
years  is  probably  safe  against  a  large  release  in  the  near  future.  Thus,  a  recent  large  event  on 
a  fault  might  actually  indicate  safety  in  the  immediate  future,  rather  than  an  indication  of 
activity.  A  single  event  by  itself  cannot  give  an  accurate  measure  of  return  time. 

Krinitzsky  (1974)  is  quoted  below: 

"In  the  United  States  the  history  of  earthquake  activity  is  greatly  truncated. 

At  best  the  record  covers  less  than  350  years.  This  it  does  in  very  few  places. 

For  most  of  the  country  it  is  about  150  years.  This,  however,  is  general 
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earthquake  history  and  not  the  record  of  movement  on  specific  faults  or  in  fault 
zones.  In  other  parts  of  the  world  the  record  for  both  earthquake  history  and 
interpreted  fault  activity  is  better. 

N.N.  Ambraseys  (1971)....  has  studied  the  record  of  damaging 
earthquakes  for  the  past  two  millennia  in  some  portions  of  the  Near  East.  [Figure 
3-5]....  shows  the  cumulative  distribution  during  this  time  of  damaging 
earthquakes  along  (a)  the  Border  Zone,  which  is  a  northeast-southwest  trend  of 
faults  that  extends  toward  the  Dead  Sea  along  the  border  of  Israel;  and  (b)  the 
Anatolian  Zone  that  is  roughly  an  east-west  zone  of  faulting.  It  intersects  with 
the  Border  Zone.  Ambraseys  points  out  that  during  the  first  five  centuries  the 
Border  Zone  was  quiescent  while  the  Anatolian  Zone  was  active.  For  the 
following  six  centuries  the  pattern  was  reversed  only  to  be  reversed  again  during 
the  eleventh  century.  In  this  case  there  is  a  factor  of  dependence  of  activity  along 
one  zone  on  that  of  the  other.  The  two  join  to  encompass  a  miniplate,  and 
movement  tends  to  shift  from  one  side  to  another  on  a  cyclical  basis.  This  tends 
to  discredit  the  validity  of  probabilistic  projections  based  only  on  a  short  history. " 

Also,  as  noted  in  Krinitzsky  (1974),  there  are  likely  to  be: 

....  "variations  in  rates  of  slippage  along  various  segments  of  any  one  long  active 
fault  or  fault  zone.  This  is  known  to  be  the  case  along  the  San  Andreas.  During 
its  relatively  short  historic  period,  it  was  noted  that  major  earthquakes  moved 
from  place  to  place  along  the  fault.  Portions  that  once  moved,  notably  the 
segment  that  slipped  during  the  San  Francisco  earthquake  of  1906,  have  become 
locked  while  slippages  occurred  elsewhere.  [Figure  3-6]  is  a  schematic  statement 
of  variations  in  rates  of  slippage  in  inches  per  year....  The  slippage  is  that  which 
has  occurred  during  an  interval  of  60  years.  Though  this  is  instructive  of 
irregularities  in  the  rates  of  movement,  it  is  not  intended  as  a  guide  to  the  future. 
If  anything,  the  future  is  very  likely  to  be  different.  The  segment  between 
Cholame  and  Camp  Dix  appears  to  be  in  a  locked  position.  Stresses  are  building 
up.  One  day  this  segment  will  rupture  suddenly." 
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Figure  3-1.  Revision  of  epicentral  location  based  on  association  with  fault 

(from  Krinitzsky,  1974). 
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Figure  3-2.  Micro  earthquakes  indicating  fault  trend,  see  Figure  3-3 

(from  Krinitzsky,  1974). 


Figure  3-3.  Major  fault  zones  in  central  California 
(rom  Krinitzsky,  1974). 
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Figure  3-4.  Epicenters  (from  US  Gelogical  Survey ,  1975) 
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Figure  3-5.  Time  distribution  of  damaging  earthquakes  in  Turkey 

(from  Ambraseys,  1971) 
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Figure  3-6.  Variation  in  rate  of  slippage  along  san  Andreas  Fault 

(from  Krinitzsky,  1974). 


CHAPTER  4 

ESTIMATING  EARTHQUAKE  RECURRENCE 


EXPONENTIAL  MAGNITUDE  DISTRIBUTION 

A  fundamental  step  in  the  estimation  of  seismic  hazard  is  the  definition  of  recurrence 
interval  of  possible  earthquakes  as  quantification  of  site  exposure.  In  1954  Guttenberg  and 
Richter  developed  an  exponential  frequency  magnitude  relationship: 


log  N(m)  =  a  -  bm 


(4-la) 


where:  N  (m)  =  cumulative  number  of  earthquakes  greater  than  m 
a  =  constant 
b  =  constant 

Equation  4- la  can  be  written  in  the  form  of  an  exponential  distribution: 


N  (m)  »  exp  (a  -pm) 


(4-lb) 


where:  a  =  a  ln(10) 

B  =  b  ln(10) 

A  lower  bound,  m,,  can  be  selected  as  an  arbitrary  reference  point.  The  following  can  be 
developed: 

N(m)  =  N(m,)  exp(-P(m-m,))  (4-2) 


where:  rrij  =  arbitrary  reference  magnitude 

Equations  4-1  and  4-2  are  constrained  by  an  upper  limit  magnitude  associated  with  the  capability 
of  a  specific  fault  to  generate  such  an  event  based  on  the  fault’s  length  and  maximum  rupture 
possible.  The  physical  limitations  of  an  upper  limit  truncate  the  magnitude  distribution. 


31 


To  estimate  frequency  with  less  error,  a  magnitude  frequency  relation  in  the  form  of: 


log10  N(m)  =  a0  *  bj  m  ♦  b2m2 


(4-3) 


has  been  used.  This  quadratic  form  has  an  advantage  of  not  overestimating  the  occurrence  of 
large  events  and  avoids  the  discontinuity  of  the  function  from  a  truncated  linear  frequency 
distribution. 

Another  generalized  formulation  of  the  Guttenberg-Richter  relation  is  of  the  form: 

In  N(m)  -  A  -  B  exp  (a  m) 


CHARACTERISTIC  MAGNITUDE 

Coppersmith  and  Youngs  (1990)  report  that  recent  geologic  studies  of  late  quaternary 
faults  strongly  suggest  that  the  exponential  recurrence  model  is  not  appropriate  for  expressing 
earthquake  recurrence  on  individual  faults.  Their  studies  suggest  many  individual  faults  tend  to 
generate  essentially  the  same  size  ot  characteristic  earthquakes  having  a  relatively  narrow  range 
of  magnitudes  at  or  near  the  maximum.  This  conclusion  is  based  upon  evaluating  the  amount 
of  displacement  per  event  for  studies  of  the  Wasatch  and  south-central  San  Andreas  faults.  This 
implies  that  earthquake  recurrence  does  not  conform  to  an  exponential  recurrence  model  but 
rather  one  that  has  a  variable  b  value.  The  type  of  geologic  recurrence  interval  data  developed 
for  the  Wasatch  and  San  Andreas  faults  are  not  generally  available  for  most  faults  in  the  Western 
United  States.  Fault  slip  rate  data  are  required. 

Youngs  and  Coppersmith  (1985)  note  that  when  geologically  derived  recurrence  intervals 
for  characteristic  earthquakes  are  compared  with  relationships  derived  from  seismicity  data,  a 
marked  mismatch  occurs.  The  characteristic  earthquake  was  found  to  include  a  band  of  events 
of  about  one-half  magnitude  width.  For  events  less  than  the  characteristic  event  magnitude,  the 
exponential  recurrence  behavior  was  found  to  be  a  satisfactory  representation.  The  increment 
between  the  minimum  characteristic  magnitude  and  the  portion  of  the  recurrence  curve  showing 
exponential  behavior  at  recurrence  rates  greater  than  the  rate  for  characteristic  events  is  about 
one  magnitude  in  width.  The  magnitude  range  showing  nonexponential  behavior  is  about  1.5 
magnitudes  in  width.  Figure  4-1  shows  the  generalized  function.  To  simplify  application,  Amc 
is  0.5  magnitude  units,  m’  is  set  at  mu  -  Amc,  and  the  value  of  n  (mc)  =  n  (m’-l). 


EVENT  RETURN  TIME 

The  Poisson  process  has  been  used  to  describe  earthquake  occurrences  and  represents  a 
basic  model  with  only  a  single  parameter  to  define.  A  Poisson  process  is  a  continuous  time, 
integer-valued  counting  process  with  stationary  independent  increments.  This  means  the  number 
of  events  occurring  in  an  interval  of  time  depends  only  on  the  length  of  the  interval.  The 
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probability  of  an  event  occurring  in  the  interval  is  independent  of  the  history.  The  probability 
distribution  of  the  number  of  earthquakes  is  given  by: 


(n,  At) 


n! 


(4-4) 


where:  PN  =  probability  of  occurrence  of  n  events  of  a  given  magnitude  range 
k  =  mean  rate  of  occurrences  per  unit  of  time  t 
The  expected  value  of  the  numbers  of  earthquakes  in  time  t  is: 

E(t)  =  kt  (4-5) 


A  characterization  of  the  Poisson  process  is  that  the  time  between  events  is  independent, 
identically  and  exponentially  distributed  with  a  constant  rate  of  occurrence,  k.  The  density 
function  of  the  time  between  events  is: 

f(t)  =  XeA‘  (46) 


A  consequence  of  Equation  4-6  is  the  hazard  rate  for  the  Poisson  process  is  constant.  In  general, 
the  hazard  rate  is  defined  as: 


h(t) 


f(t) 

1-F(t) 


(4-7) 


where  F(t)  is  the  cumulative  distribution  function  of  the  time  between  events.  The  quantity 
h(t)dt  is  the  conditional  probability  of  an  event  occurring  in  (t,  t+dt)  given  there  are  no  events 
in  the  interval  (0,t).  For  the  Poisson  process, 


h(t) 


keXt 

1  -  (1  -  e**) 


k 


(4-8) 


The  constant  hazard  rate  implies  that  the  occurrence  of  an  event  in  (t,t+dt)  is  independent  of  the 
time  since  the  last  event.  This  means  that  whether  an  earthquake  of  size  m  just  occurred  or 
whether  there  is  a  significant  gap  in  events,  the  probability  of  an  event  occurring  is  independent 
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of  past  history.  Physically,  the  energy  released  during  event  m  does  not  affect  the  reservoir  of 
stored  energy  available  for  subsequent  events. 

A  significant  advantage  of  the  Poissan  process  is  that  it  is  based  only  on  the  epicenter  data 
base.  Thus,  when  additional  information  is  lacking,  it  may  be  applied.  Other  models,  as  will 
be  discussed  in  subsequent  sections,  will  require  additional  parameters  which  may  not  be 
available,  and  will  require  substantial  effort  to  acquire. 

Weinbull  Distribution 

An  improvement  over  the  Poisson  distribution  is  the  Weinbull  distribution  suggested  by 
Chou  and  Fisher  (1975): 


P(t)  =  1  -  e'*T 


(4-9) 


where  u  and  y  are  scale  and  shape  parameters,  respectively.  Several  methods  are  available  to 
estimate  the  parameters.  The  maximum  likelihood  method  is  recommended  because  it  utilizes 
the  available  information  in  the  most  appropriate  manner.  The  shape  factor  y  is  estimated  by 
solving  the  general  equation: 
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and  the  scale  parameter  n  is  obtained  by  equating: 
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(4-11) 


where  n  is  the  size  of  the  sample  and  t,  is  the  time  interval  involved. 

A  graphical  method  of  plotting  historic  earthquake  data  is  very  useful  and  widely  used 
in  practice.  A  new  random  variable  is  introduced  as  Z  =  In  (/utv),  and. 


F(Z)  =  1  -  e  *2 


(4-12) 
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The  list  of  earthquake  occurrence  data  is  arranged  in  groups  of  intensity  or  magnitude  ranges. 
Within  each  group  of  earthquake  events  of  the  same  magnitude  range,  the  data  are  ordered  by 
time  intervals  between  occurrences,  with  the  most  frequent  first.  The  plotting  position  of  any 
data  point  within  a  group  is  [t^,  F(i)]  where  t,  is  the  i“  longest  time  interval  and, 


F(i) 


i 

n  +  1 


(4-13) 


where:  i  =  position  of  the  list;  i.e.,  first,  second,  etc. 
n  =  total  number  of  events  in  the  group 

The  parameters  p  and  y  are  determined  by  the  intercept  and  slope,  respectively,  of  the  plotted 
data  following  the  relationship: 


Z.  =  In  p  +  y  to  t4 


(4-14) 


where  t,  is  the  ith  longest  time  interval  and  Zx  is  given  by: 


lnln 


1 

1  -  F(Z) 


(4-15) 


or 


(4-16) 


When  y  =  1,  the  Weibull  distribution  is  equivalent  to  the  Poisson  distribution. 

Semi-Markov  Process 

A  Markov  or  semi-Markov  process  can  be  defined  as  a  process  in  which  the  occurrences 
of  earthquakes  make  transitions  from  one  range  of  earthquake  magnitude  to  each  of  several  other 
ranges.  The  transitions  are  probabilistic  and  have  a  one-step  memory  and  the  probability  of 
moving  to  a  given  magnitude  depends  on  the  preceding  magnitude.  A  semi-Markov  process  is 
also  characterized  by  a  probabilistic  holding  time  between  successive  transitions.  The  probability 
that  the  holding  time  between  two  successive  earthquakes  is  equal  to  a  given  value  depends  on 
the  magnitude  of  the  two  events. 

The  semi-Markov  process  is  consistent  with  the  physical  understanding  of  the  earthquake 
process.  That  process  consists  of  a  gradual  uniform  accumulation  and  periodic  release  of 
significant  strain  energy  within  short  periods  of  time  following  an  earthquake  of  large  magnitude. 
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The  occurrence  of  another  such  sized  event  at  the  same  location  is  less  likely.  As  the  time 
without  occurrence  of  a  large  magnitude  event  increases,  so  does  the  probability  of  the 
occurrence  of  such  an  event.  The  size  of  the  next  large  event  and  the  holding  time  to  that  event 
are  influenced  by  the  amount  of  strain  energy  released  in  the  previous  event  and  the  time  during 
which  strain  energy  has  been  accumulating. 

Definition  of  the  semi-Markov  model  requires  specification  of  the  most  recent  earthquake 
for  a  fault  and  the  elapsed  time  since  that  event  for  various  magnitude  levels.  Transition 
probabilities  for  each  magnitude  to  other  magnitudes  must  be  specified.  A  probability 
distribution  of  holding  times  between  the  occurrence  of  two  successive  events  must  be 
determined.  The  exposure  time  or  period  of  interest  must  be  specified. 

Woodward  Clyde  (1982)  describes  development  of  such  a  model  that  uses  a  Poisson 
process  for  events  below  a  specified  magnitude  and  a  semi-Markov  process  for  events  above  that 
level.  This  model  requires  a  Bayesian  approach  for  parameter  estimation.  The  ability  to  use  this 
class  of  model  depends  on  the  ability  to  quantify  input  data.  The  Bayesian  approach  will  be 
discussed  in  the  next  section.  The  Markov  process  utilizes  subsets  of  the  epicenter  data  base  for 
earthquake  occurrences.  This  analysis  of  windows  into  the  data  requires  a  large  catalog  of  events 
which  are  often  not  available.  For  this  reason,  semi-Markov  models  incorporate  other 
techniques. 

Bayesian  Process 

One  method  that  has  been  used  to  supplement  data  is  the  use  of  Bayes’  Theorem. 
Subjective  information  can  be  combined  with  historical  data  to  develop  parameters  for  a  seismic 
model.  Since  the  time  interval  or  holding  times  between  large  magnitude  earthquakes  may  be 
several  hundred  years,  the  historical  seismicity  data  above  are  not  sufficient  to  provide  reliable 
estimates  of  parameters  for  a  semi-Markov  model.  A  Bayesian  procedure  with  both  historical 
seismicity  data  and  subjective  expert  input  has  been  used.  From  the  seismicity  data,  estimates 
of  transition  time  can  be  made  from  statistical  examination  of  the  time  interval  of  the  magnitudes 
of  events  which  follow  an  event  of  magnitude  m;  for  all  increments  of  i.  The  holding  times 
between  earthquakes  of  magnitude  m,  and  mj  can  be  evaluated.  Subjective  expert  input  takes  the 
form  of  specification  of  fractiles  of  similar  transition  probabilities  and  holding  times. 

The  Woodward  Clyde  (1982)  study  was  done  for  a  specific  region.  Obtaining  data 
required  is  a  major  impediment  from  general  adoption  of  Bayesian  techniques  into  an  easily 
useable  general  model  for  engineering  applications. 

Kiremidjian,  et  al.  (1990)  also  developed  a  Bayesian  procedure  using  a  semi-Markov 
model.  They  note  that  at  plate  boundaries  there  are  two  types  of  forces:  relatively  uniform 
continuous  tectonic  forces  and  time-dependent  forces  from  the  asthenosphere.  The  coupling  of 
the  two  forces  causes  nonlinear  stress  accumulation.  They  developed  a  random  slip  stochastic 
model  to  represent  the  nonlinear  stress  accumulation.  To  estimate  the  interarrival  time 
distribution,  a  Weibull  distribution  is  used.  The  parameters  for  the  distribution  are  difficult  to 
determine  and  subjective  geophysical  information  is  used  to  supplement  the  data.  This  model 
has  been  applied  to  the  Middle  America  Trench  and  again  is  not  thought  adaptable  for  general 
engineering  application. 


RENEWAL  MODELS 


A  renewal  process  consists  of  a  set  of  independent,  positive  random  variables  with 
identical  distribution.  In  a  Poisson  process,  the  sequence  of  arrival  times  is  a  renewal  process 
with  exponentially  distributed  random  variables.  In  a  Markov  process,  the  time  between  entry 
to  a  fixed  state  form  a  renewal  process.  Whenever  an  event  occurs,  the  process  starts  over.  This 
process  attempts  to  characterize  the  underlying  physical  process  of  strain  buildup  and  release. 

One  approach  using  this  process  is  to  assume  a  Poisson  process  -  -ith  a  constant  occurrence 
rate  over  some  time  period.  If  an  earthquake  of  the  same  set  magniti  e  has  not  occurred  during 
this  period,  the  rate  of  occurrence  increases  to  a  larger  value.  After  an  event,  the  rate  decreases 
to  the  original  value. 

Cornell  and  Winterstein  (1988)  describe  a  renewal  model.  A  general  analytical  form  for 
the  frequency-magnitude  relation  in  the  form  of  an  approximate  Weibull  distribution  is  used: 


ln[N(m)]«oc  -  [  p  (m  -  m^m(  (4-17) 


where:  a  =  In  [N(mt)]  reflects  the  recurrence  rate  of  significant  events,  mt 

Vm  =  (var  [M]i/2/(E  [M]  -  mt)  is  the  coefficient  of  variation  of  M  -  mt 
6  =51 

The  log-linear  relation  of  the  Guttenberg-Richter  law  follows  from  Equation  4-17  when  Vm  = 
1 .  The  duration  between  events  is  not  constrained  to  be  exponentially  distributed  as  implied  by 
the  Poisson  process.  This  permits  inclusion  of  characteristic  time  or  event  models.  There  is  a 
dependence  between  time  between  events  and  the  size  of  the  last  event. 

The  characteristic  interevent  time  is  represented  by  a  Weibull  model. 

Anagnos  (1993)  developed  a  model  using  a  Markov  renewal  process  to  describe 
earthquake  states.  Since  specific  stress  release  data  are  not  available,  specific  magnitude  ranges 
are  selected  and  correlated  with  slip.  It  assumes  stress  builds  up  linearly  at  constant  rate,  and 
change  depends  only  on  the  present  state.  An  associated  semi-Markov  process  is  used  to 
establish  the  duration  of  the  visits  to  the  states.  Given  that  an  earthquake  of  given  magnitude 
just  occurred,  the  associated  semi-Markov  process  enters  a  specific  state  and  remains  in  that  state 
for  some  time  until  accumulated  stress  reaches  a  threshold  level.  The  process  then  moves  onto 
the  next  state.  Event  probabilities  are  defined  for  each  transition  state.  The  stress  release, 
occurrence  of  an  earthquake,  and  drop  to  another  state  are  random  events.  The  shape  of  the 
release  distribution  is  based  on  the  seismicity  of  the  fault  or  region  defined  by  the  Guttenberg- 
Richter  frequency-magnitude  relation.  A  characteristic  repeat  event  may  also  be  used  to  define 
the  magnitude  distribution.  A  Weibull  distribution  is  used  to  compute  the  holding  time 
distribution,  which  is  the  probability  of  the  time  to  the  next  event.  The  mean  and  variance  of 
the  Weibull  distribution  can  be  estimated  from  interarrival  time  data  and  slip  rate  data. 
Probability  of  occurrence  of  an  earthquake  greater  than  a  specific  magnitude  and  the  expected 
number  of  events  can  be  calculated  using  recursive  relations  for  the  associated  discrete  time  semi- 
Markov  process. 
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Application  of  this  model  requires  specification  of  the  lower-bound  fault  magnitude,  mean 
magnitude  and  associated  standard  deviation,  and  upper-bound  magnitude.  The  parameters  of 
the  Weibull  distribution  for  time  to  the  next  event  must  be  established.  This  may  be  based  on 
data  from  a  least  squares  regression  analysis  of  magnitude,  or  must  be  developed  for  the  specific 
site.  The  time-predictable  model  reflects  the  dependence  on  the  time  since,  and  the  size  of,  the 
last  event.  When  the  gap  in  event  occurrence  becomes  much  larger  than  the  mean  interevent 
time,  the  behavior  of  the  process  approaches  a  Poisson  process.  For  small  gaps,  the  Poisson 
process  gives  higher  probability  of  occurrence  estimates.  However,  if  the  gap  is  sufficiently 
large,  the  Poisson  model  will  underestimate  the  hazard  from  future  events.  The  time  predictable 
model  exhibits  a  correlation  between  interevent  time  and  the  preceding  magnitude  level,  although 
both  remain  independent  of  subsequent  magnitude.  The  time  predictable  renewal  process  gives 
higher  hazard  predictions  than  a  simple  renewal  process. 


DISCUSSION 

The  Poisson  process  has  a  distinct  advantage  in  its  ability  to  be  used  with  available  data 
to  characterize  its  single  parameter.  The  independence  assumptions  associated  with  the  Poisson 
process  do  not  permit  it  to  characterize  the  underlying  physical  process  of  strain  buildup  and 
release.  The  time  dependence  feature  is  lacking.  It  was  found  to  be  unconservative  when  a  long 
gap  occurred  since  the  last  event  occurred  on  a  fault.  Cornell  and  Winterstein  (1988)  give  an 
excellent  evaluation  of  the  limitations  and  applicability  of  the  Poisson  model.  They  studied  a 
broad  set  of  models  with  temporal  and  magnitude  dependence,  including  time  and  slip  predictable 
models.  They  considered  agreement  acceptable  for  engineering  hazard  studies  when  results 
agreed  within  a  factor  of  three  for  a  50-year  time  window  and  magnitude  levels  with  annual 
exceedance  probabilities  of  0.001  or  less.  According  to  Cornell  and  Winterstein  (1988): 

"The  Poisson  model  has  been  commonly  used  for  several  reasons.  These 
include:  (1)  some  successful  comparisons  of  its  predictions  with  observations..., 

(2)  rather  broad  acceptance  that,  lacking  evidence  to  the  contrary,  the  model  is  not 
unreasonable  physically  (especially  for  the  less-than-the-largest  events  that  may 
govern  hazard);  and,  more  formally,  (3)  the  fact  that  the  sum  of  non-Poissonian 
processes  may  be  approximately  Poisson.  But  perhaps  most  importantly,  it  is  the 
simplest  model  that  captures  the  basic  elements  of  the  problem.  ...Significantly, 
these  parameters  of  the  standard  hazard  model  are  those  that  the  engineering 
seismologist  commonly  estimates  and  is  therefore  best  prepared  to  specify. 

Should  an  alternative  model  be  considered,  questions  arise,  first,  as  to  which 
alternative  model  should  be  considered  and,  second,  as  to  how  in  practice  to 
estimate  both  the  additional  model  parameters  and  the  initial  conditions  (e.g.,  size 
and  time  of  the  last  significant  event)  upon  which  non-Poissonian  predictions  may 
depend.  Therefore,  the  practical  application  on  non-Poissonian  models  requires 
much  more  detailed  knowledge  of  specific  tectonic  features.  If  the  engineering 
conclusions  are  not  substantively  different,  the  implied  effort  may  not  be  justified. 

Cases  in  which  the  Poisson  estimate  is  insufficient  are  limited  practically 
to  those  in  which  the  hazard  is  controlled  by  a  single  feature  for  which  the  elapsed 
time  since  the  last  significant  event  exceeds  the  average  time  between  such  events. 
Moreover,  this  situation  creates  a  problem  only  if  there  is  reason  to  believe  that 


38 


the  fault  displays  strongly  regular,  "characteristic  time"  behavior.  In  particular, 
the  Poisson  estimate  will  generally  be  adequate  if  the  mean  interevent  time 
between  significant  events  exceeds  either  the  seismic  "gap"  (elapsed  time  since  the 
last  such  event)  or  the  length  of  the  historical  record,  whichever  is  less.  For 
strongly  regular  earthquakes,  the  mean  gap  length  under  random  entry  is  roughly 
half  of  the  mean  interevent  time;  therefore,  this  gap  with  higher  than  Poisson 
hazard  may  be  rather  unusual  in  engineering  design  practice. 

Finally,  in  many  practical  situations,  two  or  more  features  will  be 
important  hazard  contributions  at  a  particular  site.  In  these  cases,  the  combined 
hazard  is  better  estimated  by  the  Poisson  model  than  is  the  hazard  from  any  single 
feature." 

Lomnitz  (1989),  in  a  discussion  of  the  Cornell  and  Winterstein  (1988)  paper,  makes  the 
following  comment: 

"Cornell  and  Winterstein  found  that  the  fit  to  the  Poisson  model  improved 
as  the  number  of  discrete  sources  increased.  If  there  are  two  or  more  faults,  "the 
combined  hazard  is  better  estimated  by  the  Poisson  model  than  is  the  hazard  from 
any  single  feature."  This  result  is  not  altogether  unexpected.  An  elegant 
theorem...  proves  that  the  sum  of  any  number  of  random  point  processes  tends  to 
a  Poisson  process.  The  larger  the  number  of  arbitrary  component  processes,  the 
better  the  Poisson  fit. 

Hence,  it  is  not  quite  fair  to  state  that  the  effects  of  temporal  and 
magnitude  dependence  "are  ignored  in  the  conventional  Poisson  earthquake 
model"  (Cornell  and  Winterstein,  1988).  The  model  is  not  that  unsophisticated. 

The  Poisson  process  is  a  limiting  process  for  the  sum  of  many  point  processes  - 
all  of  which  exhibit  time  and  magnitude  dependence! 

For  example,  the  earthquake  hazard  in  North  China  is  governed  by  a  few 
large  faults,  both  on  land  and  in  the  Gulf  of  Bo,  plus  an  unknown  number  of 
small  faults.  The  capital  city  of  Beijing,  which  appears  to  be  somewhat  removed 
from  the  major  faults,  may  be  adequately  planned  on  the  basis  of  a  Poissonian 
earthquake  hazard  (which  in  effect  means  that  the  ground  conditions  dominate  the 
hazard).  But  I  would  be  concerned  about  applying  the  same  criterion  to  a  site  on 
the  Tangshan  Fault  -  event  though  a  significant  earthquake  occurred  on  that  fault 
as  recently  as  27  July  1976." 

For  engineering  hazard  studies  using  the  historical  data  base  and  available  slip  data,  a 
Poisson  model  may  be  used  as  a  starting  point.  Geologic  data  should  then  be  used  to  adjust 
recurrence  data  computed  from  the  historical  data  base.  Characteristic  event  data  can  easily  be 
incorporated,  and  this  should  be  done  where  required.  This  represents  an  "engineering"  solution 
for  a  range  of  studies  where  the  exposure  time  is  50  years  or  less  and  the  risk  levels  of  interest 
are  in  the  range  of  0.001 .  Studies  where  the  risk  range  is  0.0001  require  more  detailed  analysis. 
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CHAPTER  5 

GROUND  MOTION  ESTIMATION 


INTRODUCTION 

This  chapter  will  review  ground  motion  attenuation  equations  which  are  used  to  determine 
levels  of  acceleration  as  a  function  of  distance  from  the  source  and  magnitude  of  the  earthquake. 
The  increased  installation  of  strong  motion  data  recording  equipment  has  provided  an 
accumulation  of  earthquake  records.  Correlations  have  been  made  of  peak  acceleration  with 
distance  for  various  events.  These  equations  allow  us  to  estimate  the  ground  motion  at  a  site 
from  a  specified  event  and  the  uncertainty  associated  with  the  estimate.  This  estimation  is  a  key 
step  in  any  seismic  hazard  analysis.  There  are  a  number  of  attenuation  equations  that  have  been 
developed  by  various  researchers.  As  the  data  set  expands  and  more  data  points  are  available, 
agreement  among  researchers  improves. 


JOYNER  AND  BOORE 

Joyner  and  Boore  (1988)  developed  an  attenuation  equation  based  on  a  regression  analysis 
of  a  carefully  selected  set  of  events.  The  events  are  restricted  to  moment  magnitude,  M,  greater 
than  5  and  less  than  7.7  and  shallow  fault  rupture.  Their  equation  is: 


log  y  =  a  ♦  b  (M  -  6)  +  c  (M  -  6)2  ♦  d  (log  (r))  +  k  r  +  s 
r  =  (r2  +  h2)w 


(5-1) 


where  a,  b,  c,  d,  k,  s,  and  h  are  given  in  Table  5-1  for  estimating  quantities  corresponding  to 
the  randomly  oriented  horizontal  component,  and  in  Table  5-2  for  estimating  quantities 
corresponding  to  the  larger  of  the  two  horizontal  components. 

Tables  5-1  and  5-2  also  give  estimates  of  the  standard  deviation  (a)  of  an  individual 
prediction  of  log  y  using  the  equations. 
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CROUSE,  ET  AL. 


Crouse,  et  al.  (1984)  developed  the  following  equation  for  peak  horizontal  acceleration 
and  horizontal  pseudovelocity  response  at  S  percent  damping.  It  was  developed  from  data 
recorded  at  deep  soil  sites  (generally  greater  than  60  m  in  thickness)  during  shallow  crustal 
earthquakes  in  southern  California  (Crouse  1984,  1987;  Vyas,  et  al.,  1988): 


Iny  =  a  ♦  b  Ms  +  c  +  d  ln(r  +  1)  +  kr  (5-2) 


where:  y  =  peak  horizontal  acceleration  (gal)  or  horizontal  pseudovelocity  response  (cm/s) 
Ms  =  surface-wave  magnitude 

r  =  closest  distance  (km)  from  rupture  surface  to  recording  site 

Coefficients  a,  b,  c,  d,  and  k  are  given  in  Table  5-3.  Both  horizontal  components  were 
used  so  that  the  values  of  y  predicted  by  Equation  5-2  correspond  to  the  randomly  oriented 
horizontal  component. 


SADIGH,  ET  AL. 

Sadigh,  et  al.  (1986)  developed  an  equation  for  peak  horizontal  acceleration  and 
horizontal  pseudoacceleration  response  at  5  percent  damping  from  data  from  the  Western  United 
States  supplemented  by  significant  recordings  of  earthquakes  at  depth  less  than  20  km  from  other 
parts  of  the  world.  Both  horizontal  components  were  used: 


Iny  =  a  +  b  M  +  c,  (8.5  -  M)®2  «-  d  In  [r  +  h,  exp  (hj  M)]  (5-3) 


where:  y  =  peak  horizontal  acceleration  (g)  or  horizontal  pseudoacceleration  (g) 

M  =  moment  magnitude 

r  =  the  closest  distance  (km)  to  the  rupture  surface 
a,  b,  Cj,  d,  and  h4  are  given  in  Table  5-4. 

The  values  in  Table  5-4  were  derived  for  strike-slip  earthquakes.  To  obtain  estimates  for 
reverse-slip  events,  the  strike-slip  estimates  should  be  increased  by  20  percent.  Sadigh,  et  al. 
(1989)  used  additional  earthquake  data  through  1988  to  develop  a  set  of  coefficients  for  short- 
period  horizontal  ground  motion  at  rock  sites  in  reverse-slip  earthquakes  (Table  5-5).  Data  from 
both  horizontal  components  were  used  in  developing  the  equations,  and  the  results  apply  to 
reverse-slip  earthquakes.  The  results  should  be  reduced  by  17  percent  to  give  estimates  for 
strike-slip  events  and  by  9  percent  to  give  estimates  for  reverse-oblique  slip  events.  The  shape 
of  response  spectra  computed  from  Equation  5-3  does  not  change  with  distance  for  die 
coefficients  in  either  Table  5-4  or  5-5. 
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DONOVAN  AND  BORNSTEIN 


Donovan  and  Bomstein  (1978)  developed  the  following  equation  for  peak  horizontal 
acceleration  from  the  Western  United  States  data.  Both  horizontal  components  were  used: 

y  =  a  exp(b  M)  (r  «■  25)d 
a  =  2,154,000(r)'2l° 

(5-4) 

b  =  0.046  +  0.445  log(r) 
d  =  2.515  +  0.486  log(r) 


where:  y 
M 
r 


peak  horizontal  acceleration  (gal) 
any  magnitude 

distance  (km)  to  the  energy  center,  default  at  a  depth  of  5  km 


Table  5-6  gives  the  standard  deviation  of  the  natural  logarithm  of  an  individual  prediction  of  y. 


CAMPBELL 

Campbell  (1987,  1989)  developed  equations  for  estimating  peak  acceleration,  peak 
velocity,  and  pseudovelocity  at  5  percent  damping.  He  used  a  worldwide  data  set  including 
earthquakes  as  recent  as  1987  based  on  the  following  criteria: 

"(1)  the  largest  horizontal  component  of  peak  acceleration  was  at  least  0.02  g;  (2) 
the  accelerograph  triggered  early  enough  to  record  the  strongest  phase  of  shaking; 

(3)  the  magnitude  of  the  earthquake  was  5.0  or  larger;  (4)  the  closest  distance  to 
seismogenic  rupture  was  less  than  30  or  50  km,  depending  on  whether  the 
magnitude  of  the  earthquake  was  less  than  or  greater  than  6.25;  (5)  the  shallowest 
extent  of  seismogenic  rupture  was  no  deeper  than  25  km." 

Records  from  instruments  on  the  abutments  or  toes  of  dams  were  excluded,  as  were  records  from 
"hard-rock"  sites  and  shallow-soil  sites,  which  were  defined  as  sites  with  1  to  10  m  of  soil 
overlying  rock. 

Campbell’s  (1989)  equation  is: 
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In  y  =  a  +  bM+dln[r+h,  expChj  M) 

3 

-*■  q  F  fj  tanh  [£j  (M  +  f,)]  +  g,  tanMgj  D)  ♦  1,  Kj 

i-J 


(5-5) 


where:  y  =  ground  motion  parameter  of  interest,  the  vertical  component  or  the  mean  of  two 
horizontal  components. 

M  =  surface  wave  magnitude  Ms  if  both  local  magnitude  ML  and  Ms  are  greater  than 
or  equal  to  6.0  or  ML  if  both  Ms  and  ML  are  less  than  6.0. 

r  =  the  shorest  distance  (km)  to  the  zone  of  seismogenic  rupture,  identified  from 
spatial  distribution  of  aftershocks,  from  earthquake  modeling  studies,  from 
regional  crustal  velocity  studies,  and  from  geodetic  and  geologic  data. 

D  =  depth  to  basement  rock  (km). 

F  and  K  are  defined  in  Table  5-7. 

Values  for  the  coefficients  of  Equation  5-5  are  given  in  Tables  5-8  and  5-9.  Also  given 
in  Table  5-8  are  the  values  of  the  standard  deviation  of  an  individual  prediction  of  In  y.  The 
standard  deviations  are  substantially  less  than  those  in  Tables  5-1  and  5-2  after  conversion  from 
natural  to  common  logarithms.  The  shape  of  response  spectra  computed  from  Equation  5-5  does 
not  change  with  distance. 


IDRISS 

Idriss  (1985,  1987)  developed  the  following  for  the  randomly  oriented  horizontal 
component  of  peak  horizontal  acceleration: 

In  y  =  In  a  +  d  ln(r  +  20)  (5-0) 

where:  y  =  peak  horizontal  acceleration  (g) 

M  =  surface-wave  magnitude  for  M  greater  than  or  equal  to  6  and  local  magnitude 
otherwise. 

r  -  closest  distance  (km)  to  the  source  for  M  greater  than  6  and  hypocentral 
distance  otherwise. 

a  and  d  are  given  in  Table  5-10. 


45 


Idriss  proposed  that  peak  acceleration  be  used  to  scale  the  response  spectral  shapes  for  different 
site  conditions  with  magnitude-  and  period-dependent  correction  factors.  The  shape  of  response 
spectra  computed  by  Idriss’  method  does  not  change  with  distance. 


COMPARISON  OF  EQUATIONS 

According  to  Joyner  and  Boore  (1988): 

"...  .to  properly  compare  the  different  relationships,  adjustments  must  be  made  for 
the  different  definitions  of  distance.  [Figure  5-1]  compares  peak  horizontal 
acceleration  for  the  randomly  oriented  horizontal  component  at  magnitude  6.5  as 
estimated  by  Donovan  and  Bomstein  (1978),  Joyner  and  Boore  (1988),  Idriss 
(1987),  and  Campbell  (1989).  The  definition  of  distance  used  in  [Figure  5-1]  is 
the  closest  distance  to  the  vertical  projection  of  the  rupture  on  the  surface  of  the 
earth.  The  curves  of  Donovan  and  Bomstein  and  Campbell  were  adjusted 
assuming  a  source  depth  of  5  km.  The  curve  shown  for  Idriss  is  that  for  deep  soil 
sites.  The  curve  shown  for  Campbell  is  that  for  strike-slip  earthquakes  recorded 
at  free-field  sites. 

At  short  distances,  where  it  matters  the  most,  the  different  relationships 
agree  to  within  a  fraction  of  the  uncertainty  of  an  individual  estimate  as  given  by 
any  of  the  authors.  This  suggests  that  the  short-distance  estimates  at  magnitude 
6.5  are  controlled  by  the  data.  The  differences  at  large  distances  are  not  of  much 
practical  importance." 

Figure  5-1  gives  the  same  comparison  for  magnitude  7.5.  The  agreement  at  short  distance  is  not 
as  good  as  at  magnitude  6.5,  reflecting  the  scarcity  of  data  points,  but  it  is  within  the  uncertainty 
of  an  individual  estimate. 

From  Figure  5-1,  the  Donovan  and  Bomstein  equation  yields  an  almost  upper  bound  and 
is  a  conservative  estimate.  This  equation  was  selected  for  further  study.  The  study  included  the 
following  attenuation  equations: 

1.  McGuire  (1978) 

2.  Trifunac  and  Brady  (1975) 

3.  Campbell  (west)  (1982) 

4.  Campbell  (east)  (1982) 

5.  Donovan  and  Bomstein  (1978) 

6.  Joyner  and  Boore  (1988) 

Figure  5-2  shows  the  attenuation  relationships  along  with  events  selected  from  the  catalog 
of  strong  motion  records  cited  by  Chang  (1978).  Records  on  buildings  and  vertical  records  were 
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excluded.  As  can  be  seen,  there  is  significant  scatter.  The  standard  deviation  in  equation 
prediction  to  date  is  given  as: 


Equation 

Mean 

Deviation*1 

Standard 

Deviation 

McGuire 

-0.005 

0.089 

Trifunac  and  Brady 

+0.008 

0.105 

Campbell  (west) 

-0.059 

0.111 

Campbell  (east) 

-0.041 

0.101 

Donovan  and  Bomstein 

-0.042 

0.104 

Joyner  and  Boore 

-0.056 

0.111 

■"Minus  =  equation  underes:: mates. 


Another  data  set  used  by  Joyner  and  Boore  (1988)  is  shown  in  Figure  5-3.  The  mean  and 
standard  deviations  of  the  equations  are: 


Equation 

Mean 

Deviation* 

Standard 

Deviation 

McGuire 

+0.001 

0.0901  | 

Trifunac  and  Brady 

+0.076 

0.1885  I 

Campbell  (west) 

-0.027 

Campbell  (east) 

-0.013 

0.0913 

Donovan  and  Bomstein 

-0.014 

0.0844 

Joyner  and  Boore 

-0.021 

0.0830 

■"Minus  =  equation  underestimates. 


The  McGuire  equation  has  the  least  overall  bias.  The  standard  deviations  except  for  the 
Trifunac  and  Brady  equations  are  statistically  about  the  same.  From  the  above,  there  is  little 
difference  between  McGuire’s  equation  and  Donovan’s  equation. 

The  significance  of  selection  of  the  earthquake  attenuation  equation  is  primarily  in  the 
determination  of  separation  distances.  McGuire  and  Donovan,  for  example,  define  the 
relationship  in  hypocentral  distance;  Trifunac  and  Brady  in  terms  of  epicentral  distance;  and 
Campbell  and  Joyner  and  Boore  in  terms  of  the  closest  distance  to  the  fault.  For  risk  analysis 
of  faults  where  the  site  tends  to  be  located  near  one  end  of  the  fault,  the  selection  of  a 
relationship  based  on  the  shortest  distance  will  tend  to  give  higher  values.  To  illustrate  the 
significance  of  the  attenuation  relationship,  a  study  was  made  of  a  site  0.5  miles  from  the 
Hayward  fault.  Results  for  the  225-year  return  time  are: 


Joyner  and  Boore 
Campbell 
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0.535  g 
0.540  g 


McGuire  0.396  g 

Donovan  and  Bomstein  0.382  g 

The  equations  shown  in  Figure  5-1  represent  the  current  consensus  of  top  researchers  in 
the  field  and  show  good  agreement. 


DATA  FROM  THE  LOMA  PRIETA  EARTHQUAKE 

The  Loma  Prieta,  California  earthquake  of  17  October  1989  brought  a  significant  addition 
to  the  strong  motion  data  base.  Figure  5-4,  with  equations  from  Joyner  and  Boore  (1988)  and 
Loma  Prieta  data,  shows  the  larger  of  the  two  horizontal  components  of  peak  acceleration  plotted 
against  the  closest  distance  to  the  vertical  projection  of  the  rupture  on  the  surface  of  the  earth. 
The  solid  curves  in  Figure  5-4  are  Joyner  and  Boore  (1988)  for  the  larger  of  two  horizontal 
components  for  a  moment  magnitude  of  6.9  (Kanamori  and  Helmberger,  1990),  and  the  dashed 
curves  are  the  curves  corresponding  to  plus  and  minus  one  standard  deviation  of  an  individual 
estimate.  The  recorded  data  are  plotted  with  different  symbols  for  rock,  alluvium,  and  bay-mud 
sites.  The  values  for  the  three  site  categories  are  factors  of  1.6,  1.8,  and  4.5  higher  on  the 
average  than  the  estimates  from  the  equations,  which  were  derived  from  a  data  set  that  included 
rock  and  alluvium  sites,  but  not  bay-mud  sites.  The  larger  values  for  rock  and  alluvium  may 
represent  ordinary  earthquake-to-earthquake  variability,  but  the  values  for  bay  mud  clearly 
represent  a  local  site  effect. 


DISCUSSION 

Table  5-11  summarizes  the  earthquake  attenuation  equations  reviewed  above.  The 
earthquake  epicenter  data  base  records  a  surface  magnitude,  a  local  magnitude,  a  body 
magnitude,  and  another  magnitude.  There  is  uncertainty  in  each  magnitude  since  events  often 
do  not  have  all  magnitudes  established.  Averages  of  the  individual  magnitudes  from  each  station 
are  made  to  get  a  more  uniform  value.  Specification  of  a  moment  magnitude  would  complicate 
use  of  the  existing  epicenter  data  base.  All  the  equations  give  similar  values  (Figure  5-1). 
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TABLE  5-1.  COEFFICIENTS  IN  THE  EQUATIONS  OF  JOYNER  AND  BOORE 
FOR  THE  RANDOMLY  ORIENTED  HORIZONTAL  COMPONENT  OF 
_  GROUND  MOTION 


Period  (») 

a 

6 

c 

h 

d 

k 

» 

Vso 

e 

*\otv 

Pseudovelocity  response  (cm 

/«) 

0.1 

2.16 

0.25 

-0.06 

11.3 

-1.0 

-0.0073 

-0.02 

0.28 

0.15 

2.40 

.30 

-  .08 

10.8 

-1.0 

-  .0067 

-  .02 

.28 

0.2 

2.46 

.35 

-  .09 

9.6 

-1.0 

-  .0063 

-  .01 

.28 

0.3 

2.47 

.42 

-  .11 

6.9 

-1.0 

-  .0058 

.04 

590 

-0.28 

.28 

0.4 

2.44 

.47 

-  .13 

5.7 

-1.0 

-  .0054 

.10 

830 

-  .33 

.31 

0.5 

2.41 

.52 

-  .14 

5.1 

-1.0 

-  .0051 

.14 

1020 

-  .38 

.33 

0.75 

2.34 

.60 

-  .16 

4.8 

-1.0 

-  .0045 

.23 

1410 

-  .46 

.33 

1.0 

2.28 

.67 

-  .17 

4.7 

-1.0 

-  .0039 

.27 

1580 

-  .51 

.33 

1.5 

2.19 

.74 

-  .19 

4.7 

-1.0 

-  .0026 

.31 

1620 

-  .59 

.33 

2.0 

2.12 

.79 

-  .20 

4.7 

-1.0 

-  .0015 

.32 

1620 

-  .64 

.33 

3.0 

2.02 

.85 

-  .22 

4.7 

-0.98 

.0 

.32 

1550 

-  .72 

.33 

4.0 

1.96 

0.88 

-0.24 

4.7 

-0.95 

0.0 

0.29 

1450 

-0.78 

0.33 

Peak  acceleration  (g) 

0.43 

0.23 

0.0 

8.0 

-1.0 

-0.0027 

0.0 

0.28 

Peak  velocity  ( 

cm/s) 

2.09 

0.49 

Til 

mam 

-1.0 

-0.0026 

0.17 

1190 

-0.45 

0.33 

©  Reprinted  from  Joyner  and  Boore  (1988)  with  permission  from  American  Society  of  Civil  Engineers 


TABLE  5-2.  COEFFICIENTS  IN  THE  EQUATIONS  OF  JOYNER  AND  BOORE 
FOR  THE  LARGER  OF  TWO  HORIZONTAL  COMPONENTS 
OF  GROUND  MOTION _ 


Period  (s) 

a 

b 

c 

h 

d 

* 

s 

Vso 

e 

°io»V 

Pseudovelocity  response  (cm/s) 

0.1 

2.24 

10.6 

-1.0 

-0.0067 

-0.06 

0.27 

0.15 

2.46 

.34 

-  .10 

10.3 

-1.0 

-  .0063 

-  .05 

.27 

2.54 

.37 

-  .11 

9.3 

-1.0 

-  .0061 

-  .03 

.27 

2.56 

.43 

-  .12 

7.0 

-1.0 

-  .0057 

.04 

650 

-0.20 

.27 

0.4 

2.54 

.49 

-  .13 

5.7 

-1.0 

-  .0055 

.09 

870 

-  .26 

.30 

2.53 

.53 

-  .14 

5.2 

-1.0 

-  .0053 

.12 

1050 

-  .30 

.32 

2.46 

.61 

-  .15 

4.7 

-1.0 

-  .0049 

.19 

1410 

-  .39 

.35 

2.41 

.66 

-  .16 

4.6 

-1.0 

-  .0044 

.24 

1580 

-  .45 

.35 

1.5 

2.32 

.71 

-  .17 

4.6 

-1.0 

-  .0034 

.30 

1780 

-  .53 

.35 

.75 

-  .18 

4.6 

-1.0 

-  .0025 

.32 

1820 

-  .59 

.35 

2.17 

.78 

-  .19 

4.6 

-1.0 

.0 

.29 

1620 

-  .67 

.35 

.  Euk 

4.6 

•0.98 

.  0.0 

0.24 

1320 

-0.73 

0.35 

Peak  acceleration  (y) 

0.49 

0.23 

0.0 

8.0 

-1.0 

-0.0027 

0.0 

0.28 

Peak  velocity 

(cm/s) 

2.17 

0.49 

HZ3 

-1.0 

-0.0026 

0.17 

1190 

-0.45 

0.33 

Reprinted  from  Joyner  and  Boore  (1988)  with  permission  from  American  Society  of  Civil  Engineers 
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TABLE  5-3.  COEFFICIENTS  IN  THE  EQUATIONS  OF  CROUSE 
_ HORIZONTAL  COMPONENT  OF  GROUND  MOTION 


Period  (s) 

a 

b 

e 

d 

k 

Pseudovelocity  response  (an/*) 

0.05 

-  2.44178 

0.84826 

-0.02579 

0.59914 

0.10 

-  0.61623 

-  .00999 

.68673 

0.20 

-  4.47801 

-  .11673 

.64716 

0.40 

1.17453 

•  .04411 

•  .47398 

.62089 

0.60 

-  6.02161 

2.66493 

-  .15619 

-  .52586 

.62275 

1.00 

-  5.89916 

2.48235 

-  .13036 

-  .52261 

.62745 

2.00 

-11.48576 

4.01914 

-  .23152 

-  .56791 

.63277 

2.50 

•12.33454 

4.15828 

-.23359  , 

-  .56280 

.66459 

4.00 

-14.90528 

4.54962 

-  .24999 

•  .32351 

6.00 

-14.77796 

4.33959 

-0.23491 

•0.20849 

Bull 

0.79595 

Peak  acceleration  (?) 

2.48456 

0.73377 

•0.01509 

•0.50558 

•0.00935 

0.58082 

(P  Reprinted  from  Joyner  and  Boore  (1988)  with  permission  from  American  Society  of  Civil  Engineers 


TABLE  5-4.  COEFFICIENTS  IN  THE  EQUATIONS  OF  SADIGH 


M  <  6.5 

M  >6.5 

Period  (s) 

a 

b  Cl 

ca  d 

hi 

hi 

<hnif 

hi 

•Tin  » 

Pseudoacceleration  response  (g)  at  soil  sites 

0.1 

-2.024 

1.1  0.007 

2.5  -1.75 

0.8217 

0.4814 

1.332  -  0.148M 

0.3157 

0.6286 

0.37 

0.2 

-1.696 

1.1  .0 

2.5  -1.75 

.8217 

.4814 

1.453  -  0.162M 

.3157 

.6286 

.40 

0.3 

-1.638 

1.1  -  .008 

2.5  -1.75 

.8217 

.4814 

1.486  -  0.164M 

.3157 

.6286 

.42 

0.5 

-1.659 

1.1  -  .025 

2.5  -1.75 

.8217 

.4814 

1.584  -  0.176M 

.3157 

.6286 

.44 

1.0 

-1.975 

1.1  -  .060 

2.5  -1.75 

.8217 

.4814 

1.62  -  0.18M 

.3157 

.6286 

.45 

2.0 

-2.414 

1.1  -  .105 

2.5  -1.75 

.8217 

.4814 

1.62  -  0.18M 

.3157 

.6286 

.45 

4.0 

-3.068 

1.1  -0.160 

2.5  -1.75 

0.8217 

0.4814 

1.62  -  0.18M 

0.3157 

0.6286 

0.45 

Peak  acceleration  (g)  at  soil  sites 

-2.611 

1.1  0.0 

2.5  -1.75 

0.8217 

0.4814 

1.26  -  0.14M 

0.3157 

0.6286 

0.35 

Pseudoacceleration  response  (g)  at  rock  sites 

0.1 

-0.688 

1.1  0.007 

2.5  -2.05 

1.353 

0.406 

1.332  -  0.148M 

0.579 

0.537 

0.37 

0.2 

-0.479 

1.1  -  .008 

2.5  -2.05 

1.353 

.406 

1.453  -  0.162M 

.579 

.537 

.40 

0.3 

-0.543 

1.1  -  .018 

2.5  -2.05 

1.353 

.406 

1.486  -  0.164M 

.579 

.537 

.42 

0.5 

-0.793 

1.1  -  .036 

2.5  -2.05 

1.353 

.406 

1.584 -0.176M 

.579 

.537 

.44 

1.0 

-1.376 

1.1  -  .065 

2.5  -2.05 

1.353 

.406 

1.62  -  0.18M 

.579 

.537 

.45 

2.0 

-2.142 

l.l  -  .100 

2.5  -2.05 

1.353 

.406 

1.62  -  0.18M 

.579 

.537 

.45 

4.0 

-3.177 

1.1  -0.150 

2.5  -2.05 

1.353 

0.406 

1.62  —  0.18M 

0.579 

0.537 

0.45 

Peak  acceleration  (g)  at  rock  sites 

-1.406 

1.1  0.0 

2.5  -2.05 

1.353 

0.406 

1.26  —  0.14M 

0.579 

0.537 

0.35 

0  Reprinted  from  Joyner  and  Boore  (1988)  with  permission  from  American  Society  of  Civil  Engineers 
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TABLE  5-5.  COEFFICIENTS  IN  THE  ALTERNATIVE  EQUATIONS 
OF  SADIGH  ET  AL.(  1989}  FOR  THE  RANDOMLY  ORIENTED 
HORIZONTAL  COMPONENT  OF  SHORT-PERIOD 


GROUND  MOTION  AT  ROCK  SITES  IN  REVERSE-FAULT  EARTHQUAKES 


Period  (s) 

a 

b 

Cl 

ej 

d 

hi 

Pseudoacceleration  response  (g)  M  <  6.5 

0.1 

0.218 

mm 

0.006 

2.5 

-2.1 

3.656 

0.250 

1.29  —  0.14M 

0.2 

wm 

1.00 

-0.010 

2.5 

-2.1 

3.656 

0.250 

1.31  —  0.14M 

1.00 

-0.023 

2.5 

-2.1 

3.656 

0.250 

1.33  -  0.14M 

Wm 

1.00 

-0.040 

2.5 

-2.1 

3.656 

0.250 

1.35-0.14M 

1.0 

BH 

1.00 

-0.064 

2.5 

-2.1 

3.656 

0.250 

1.36  —  0.14M 

Peak  acceleration  ( g ) 

M  <  6.5 

-0.442 

1.00 

-2.1 

3.656 

0.250 

1.27  —  0.14M 

Pseudoacceleration  response  (g)  M  >  6.5 

0.1 

-0.432 

1.10 

0.006 

2.5 

-2.1 

0.616 

0.524 

0.38 

0.2 

-0.232 

1.10 

-0.010 

2.5 

-■2.1 

0.616 

0.524 

0.40 

-0.248 

1.10 

-0.023 

2.5 

-2.1 

0.616 

0.524 

0.42 

0.5 

-0.469 

1.10 

-0.040 

2.5 

-2.1 

0.616 

0.524 

0.44 

-1.059 

1.10 

-0.064 

2.5 

-2.1 

0.616 

0.524 

0.45 

Peak  acceleration  (9)  M  >  6.5 

-1.092 

1.10 

-2.1 

0.616 

0.524 

0.36 

Reprinted  from  Joyner  and  Boore  draft  U.S.  Geological  Survey  paper  in  publication. 


TABLE  5-6  .TANDARD  DEVIATION  GIVEN  BY  DONOVAN  AND  BORNSTEIN  (1978) 
FOR  THE  NATURAL  LOGARITHM  OF  AN  INDIVIDUAL  ESTIMATE 
OF  PEAK  HORIZONTAL  ACCELERATION 


Peak  acceleration 

0.01 

0.05 

0.10 

0.15 

Standard  deviation  of 

natural  logarithm  of 
peak  acceleration 

0.50 

0.48 

0.46 

0.41 

(€)  Reprinted  from  Joyner  and  Boore  (1988)  with  permission  from  American  Society  of  Civil  Engineers 


TABLE  5-7.  DEFINITION  OF  VARIABLES  IN  THE  EQUATIONS  OF  CAMPBELL  (1989) 


Fault  type 


1  reverse 
0  strike-slip 


Building  effects 


%  _  1  embedded  buildings  3-11  stories 
1  0  other 

ft  _  1  embedded  buildings  greater  than  11  stories 

0  other 

l  =  1  nonembedded  buildings  greater  than  2  stories 

3  0  other 


Reprinted  from  Joyner  and  Boore  draft  U.S.  Geological  Survey  paper  in  publication. 
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TABLE  5-8.  COEFFICIENTS  IN  THE  EQUATIONS  OF  CAMPBELL  (1989) 


Period  (s) 

a  6 

hj  d  9/1  h 

h 

01  02 

&ln  y 

Mean  horizontal  pseudovelocity  response  at  5  percent  damping  (cm/s) 

-0.648  1.08 

0.311 

0.597  -1.81  0.382 

0.42 

-0.379  1.08 

.311 

.597  -1.81  .382 

.44 

0.251  1.08 

.311 

.597  -1.81  .382 

.46 

0.754  1.08 

.311 

.597  -1.81  .382 

.48 

mSm 

1.424  1.08 

.311 

.597  -1.81  .382 

.50 

1.788  1.08 

.311 

.597  -1.81  .382 

.50 

/■S'-i 

2.170  1.08 

.311 

.597  -1.81  .382 

.50 

m 

2.009  1.08 

.311 

.597  -1.81  .382  0.425  0.570 

-4.7 

.50 

m  : 

1.930  1.08 

.311 

.597  -1.81  .382  0.685  .570 

-4.7 

.50 

0.75 

1.612  1.08 

.311 

.597  -1.81  .382  1.27  .570 

-4.7 

.50 

1.0 

1.268  1.08 

.311 

.597  -1.81  .382  1.74  .570 

-4.7 

.50 

1.5 

0.487  1.08 

.311 

.597  -1.81  .382  2.43  .570 

-4.7 

0.344  0.553 

.50 

0.040  1.08 

.311 

.597  -1.81  .382  2.83  .570 

-4.7 

.469  .553 

.50 

-0.576  1.08 

.311 

.597  -1.81  .382  3.17  .570 

-4.7 

.623  .553 

.50 

mm 

-0.766  1.08 

0.311 

0.597  -1.81  0.382  3.08  0.570 

-4.7 

0.857  0.553 

0.50 

Mean  horizontal  peak  acceleration  (9) 

-2.470  1.08 

0.311 

0.597  -1.81  0.382 

0.421 

Mean  horizontal  peak  velocity  (cm/s) 

-1.974  1.34 

0.00935 

1.01  -1.32  0.327 

1.16  0.0776 

0.395 

Vertical  pseudovelocity  response  at  5  percent  damping  (cm/s) 

0.04 

-2.082  0.978 

0.0536 

0.674  -1.45  0.239 

0.62 

-1.634  .978 

.0536 

.674  -1.45  .239 

.62 

-0.903  .978 

.0536 

.674  -1.45  .239 

.62 

0.1 

-0.488  .978 

.0536 

.674  -1.45  .239 

.62 

0.15 

-0.125  .978 

.0536 

.674  -1.45  .239 

.62 

0.2 

0.157  .978 

.0536 

.674  -1.45  .239 

.62 

0.3 

0.356  .978 

.0536 

.674  -1.45  .239 

.62 

0.4 

0.188  .978 

.0536 

.674  -1.45  .239  0.214  0.546 

-4.7 

.62 

0.5 

0.038  .978 

.0536 

.674  -1.45  .239  0.435  .546 

-4.7 

.62 

0.75 

-0.035  .978 

.0536 

.674  -1.45  .239  0.719  .546 

-4.7 

.62 

1.0 

-0.448  .978 

.0536 

.674  -1.45  .239  1.37  .546 

-4.7 

.62 

1.5 

-1.287  .978 

.0536 

.674  -1.45  .239  2.18  .546 

-4.7 

0.344  0.553 

.62 

2.0 

-1.580  .978 

.0536 

.674  -1.45  .239  2.36  .546 

-4.7 

.469  .553 

.62 

3.0 

-1.741  .978 

.0536 

.674  -1.45  .239  2.24  .546 

-4.7 

.623  .553 

.62 

4.0 

-1.975  0.978 

0.0536 

0.674  -1.45  0.239  2.46  0.546 

-4.7 

0.857  0.553 

0.62 

Vertical  peak  acceleration  (9) 

-4.003  0.978 

0.0536 

0.674  -1.45  0.239 

0.569 

Vertical  peak  velocity  (cm/s) 

-4.336  1.72 

0.00594 

1.14  -1.51  0.337 

0.520 

Reprinted  from  Joyner  and  Boore  draft  U.S.  Geological  Survey  paper  In  publication. 
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TABLE  5-9.  ADDITIONAL  COEFFICIENTS  FOR  BUILDING  EFFECTS 
_ IN  THE  EQUATIONS  OF  CAMPBELL  (1989) _ 

Mean  horizontal  component  Vertical  component 

Period  (s)  /i  /2  /a  /i  f2  I3 


Pseudovelocity  response  at  5  percent  damping  (cm/s) 


0.04 

-0.180 

-0.489 

BSE9H 

-0.180 

-0.489 

-0.083 

-0.712 

-0.180 

-0.489 

-0.206 

-0.582 

-0.371 

0.1 

-0.180 

-0.489 

-0.197 

-0.650 

-0.370 

0.15 

-0.180 

-0.489 

-0.392 

0.2 

-0.180 

-0.489 

-0.392 

0.3 

-0.180 

-0.489 

-0.392 

0.4 

-0.180 

-0.489 

-0.347 

0.5 

-0.180 

-0.489 

-0.153 

0.75 

-0.180 

-0.489 

-0.347 

-0.180 

-0.219 

-0.278 

1.5 

-0.180 

0.074 

0.284 

0.619 

-0.180 

0.072 

0.437 

0.992 

0.218 

0.391 

0.663 

0.291 

0.691 

1.15 

0.503 

0.759 

0.085 

0.722 

1.10 

Peak  acceleration  (g) 

-0.180 

-0.489 

-0.392 

Peak  velocity  (cm/s) 


_ _ 0.366 _ 0,388 

Reprinted  from  Joyner  and  Boore  draft  U.S.  Geological  Survey  paper  in  publication. 


TABLE  5-10.  PARAMETERS  IN  THE  EQUATIONS  OF  IDRISS  (1987)  FOR 
THE  RANDOMLY  ORIENTED  HORIZONTAL  COMPONENT  OF  PEAK  ACCELERATION  (g) 


M 

Rock  and  stiff  soil  sites 
a  d 

Deep  soil  sites 
a  d 

ffln  y 

4.5 

606 

-2.57 

189 

-2.22 

0.70 

5.0 

617 

-2.46 

195 

-2.13 

.58 

5.5 

452 

-2.28 

147 

-1.97 

.48 

6.0 

282 

-2.07 

98 

-1.79 

.42 

6.5 

164 

-1.85 

61.6 

-1.60 

.38 

7.0 

91.7 

-1.63 

37.2 

-1.41 

.35 

7.5 

49.8 

-1.41 

22 

-1.22 

.35 

8.0 

28.5 

-1.21 

13.7 

-1.05 

.35 

8.5 

15.9 

-1.01 

8.4 

-0.88 

0.35 

£)  Reprinted  from  Joyner  and  Boore  (1988)  with  permission  from  American  Society  of  Civil  Engineers 
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Table  5-11 

Comparison  of  Equations 


Equation 

Magnitude 

Distance 

a 

Problems 

Joyner  & 
Boore 

M0 

Shortest 

distance 

0.28 

Range  magnitude 
limited  5  to  7.5 

Crouse 

Ms 

Shortest 

distance 

0.58 

Uses  Ms 

Sadigh 

M0 

Shortest 

distance 

0.35  to 
1.26-0. 14M 

Uses  M0 

Donovan 

M 

Dist.  to 
hypocenter 

0.4  to  0.5 

Campbell 

ml/ms 

Shortest 

distance 

0.421 

Uses  Ms  and  ML 

Idriss 

Ms 

Shortest 

distance 

0.35  to 

0.70 

Range  limited 

M  =  1.5  to  8.5 

=  Moment  magnitude 
=  Surface  magnitude 
=  Local  magnitude 
=  Any  magnitude 
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o 


o 

M  0.01 


z 

o 

fsj  0.01 


DISTANCE,  KM 


Figure  5-1. 

Comparison  of  different  reUtiooshipa  for  peak  harisontal  acceleration  at  magnitude  6.5  (a)  and  7.3 
(b).  OB,  from  Donovan  and  Bernstein  (1978);  I,  from  Idris*  (1987)  for  deep  soil  eitee;  JB,  from  Joyner  and  Boot* 
(1982),  rodneed  by  13  percent  ao  as  to  approximate  the  vaht*  for  the  randomly  nri  rnttd  horizontal  component  ;C.  from 
Campbell  (1989)  for  a  strike-slip  — —w-p.v.  recced ed  at  a  free- field  site.  For  magnitude  6-3  the  cum  from  Cambell 
is  coincident  with  that  from  Joyner  and  Boost.  The  distance  plotted  is  the  doe  set  distance  to  the  vertical 

projection  of  the  rupture  on  the  surface  of  the  earth.  The  rams  of  Donovan  and  Bornetem  and  those  of  Campbell 
are  adjusted  assuming  a  source  depth  of  S  km. 


Reprinted  from  Joyner  and  Boore  (1988)  with  permission  from  American  Society  of  Civil  Engineers 
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8 


Figure  5-3.  Attenuation  equation  and  strong  motion  data,  magnitude  6.5. 

Data  set  from  Joyner  and  Boore 


2 


0  1 
Log  Distance  (km) 

Figure  5-4.Loma  Prieta  earthquake  and  Joyner  and  Boore 
equation  (from  draft  USGS  paper) 
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CHAPTER  6 

DEVELOPING  A  SEISMIC  MODEL 


INTRODUCTION 

This  chapter  discusses  the  process  of  constructing  a  seismic  model  and,  using  data 
presented  in  previous  chapters,  describes  a  new  procedure  for  computing  the  probability  of  site 
ground  motion.  As  mentioned  previously,  our  objective  is  to  develop  an  engineering  procedure 
that  will  give  results  of  sufficient  accuracy  to  estimate  site  motion  for  events  with  return  time  on 
the  order  of  1 ,000  years.  As  an  engineering  procedure,  it  is  intended  to  use  historical  data  and 
geologic  data  where  available.  Procedures  that  require  panels  of  experts  or  extensive  geologic 
investigation  are  not  considered  practical  or  feasible  for  limited  scope  engineering  investigations. 
Nonetheless,  given  that  such  data  exist,  the  conclusions  regarding  seismic  activity  should  be 
useable. 


BUILDING  A  SEISMIC  SOURCE  MODEL 

As  noted  by  Coppersmith  (1991),  many  elements  of  seismic  source  characterization 
depend  on  the  tectonic  environment.  In  the  Western  United  States,  the  tectonic  environment  is 
such  that  earthquakes  are  associated  on  known  faults.  However,  in  the  Eastern  United  States, 
the  causative  geologic  structures  are  generally  not  known.  A  seismic  model  must  be  based  on 
the  knowledge  of  the  local  area.  It  can  consist  of  an  area  source  zone  for  eastern  sites  or  a 
detailed  fault  definition  region  for  western  sites. 

In  the  Western  United  States,  it  is  recognized  that  large  earthquakes  are  associated  with 
faults.  For  a  magnitude  8  event,  a  rupture  of  200  miles  is  required  to  release  that  level  of 
energy.  A  200-mile  fault  exhibits  visual  evidence  of  its  existence  and  is  unlikely  to  remain 
undiscovered.  It  is  possible  for  lower  magnitude  events  that  require  considerably  less  fault 
rupture  to  occur  on  faults  lacking  recent  evidence  of  activity,  or  on  faults  that  have  not  been 
identified. 

Where  faults  are  identified  as  sources,  the  area  contained  within  the  source  zone  is 
defined  to  have  relatively  uniform  seismic  potential  in  terms  of  maximum  magnitude  and  event 
recurrence.  A  fault  is  modeled  as  a  line  source  encompassing  a  distance  or  region  surrounding 
the  fault,  such  that  the  activity  of  this  region  can  be  associated  with  events  on  the  fault.  Where 
a  fault  exhibits  variations  of  activity  along  its  length,  it  can  be  divided  into  subelements 
containing  regions  where  activity  is  uniform. 

For  the  procedures  developed  here,  a  fault  consists  of  two  line  segments  defined  by  three 
points.  The  events  to  include  or  associate  with  the  fault  are  defined  by  specification  of  a  distance 
from  the  fault  line,  such  that  all  those  events  within  the  distance  are  grouped  with  the  fault. 
Alternatively,  a  region  can  be  designated  by  four  points  to  bound  the  fault.  Again,  note  a  fault 
can  be  divided  into  pieces  where  activity  or  geometry  so  dictates. 
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In  the  Eastern  United  States,  faulting  may  not  be  readily  identifiable.  Source  zones  are 
specified  as  regions  where  a  zone  of  like  seismicity  is  evident.  The  regional  geology  and 
tectonics  assist  in  defining  the  source  zone  boundaries.  A  source  zone  is  defined  as  a  region  of 
uniform  seismicity,  such  that  an  event  is  equally  likely  to  occur  in  any  portion  of  the  zone.  This 
is  characterized  by  the  concept  of  a  "floating  earthquake,"  an  event  that  can  occur  anywhere  in 
the  zone. 

In  the  development  of  a  site  model,  it  is  important  to  keep  in  mind  that  an  equivalent 
representation  of  a  region  is  being  created  by  a  series  of  fault  line  segments  or  source  zones. 
The  seismicity  must  be  captured  in  terms  of  its  spatial  location  and  in  terms  of  the  level  of 
activity.  Assignment  of  events  to  one  fault  source  as  opposed  to  another  increases  that  fault’s 
contribution  to  the  estimation  of  event  recurrence.  It  is  important  to  capture  all  the  seismicity. 
For  faulting  conditions  where  there  are  a  number  of  parallel  elements,  it  may  not  be  easy  to 
separate  which  events  are  associated  with  which  fault.  Consideration  must  be  given  to  the  dip 
of  the  fault  in  assigning  events,  since  the  epicenter  for  a  sloping  fault  can  actually  occur  in  a 
number  of  kilometers  away  from  the  surface  trace  of  the  fault.  The  large  majority  of  strike  slip 
faults  have  steep  dips  of  70  degrees  or  greater.  On  the  other  hand,  thrust  faults  generally  have 
dips  much  less  than  this,  generally  in  the  range  of  45  to  60  degrees. 

For  the  cases  where  a  fault  is  close  to  a  site  (within  10  miles),  special  considerations 
should  be  given  to  the  location  of  the  fault  line  segments  that  define  the  fault  model.  If  the  fault 
dips  toward  the  site,  the  actual  epicentral  distance  may  be  closer  to  the  site  than  the  surface  trace 
of  the  fault.  In  this  case,  it  would  be  prudent  to  move  the  line  segments  closer  to  the  site  so  they 
represent  the  epicentral  location,  rather  than  the  surface  trace.  For  faults  at  greater  distances, 
the  difference  becomes  less  significant.  The  three-dimensional  effects  of  the  inclined  faulting 
are  captured  without  having  to  resort  to  a  full  three-dimensional  model,  since  most  western 
events  occur  at  5-  to  10-km  depths. 

Once  a  fault  or  region  has  been  defined  as  a  seismic  source,  the  maximum  earthquake 
magnitude  must  be  defined.  In  a  previous  section,  a  plot  was  shown  relating  fault  rupture  length 
to  magnitude.  The  length  of  a  fault  can  be  estimated  from  maps.  An  assumption  can  be  made 
that  a  fault  will  rupture  over  50  to  80  percent  of  its  length.  This  estimate  of  rupture  distance  can 
be  used  to  define  the  fault  magnitude.  Estimates  of  fault  magnitudes  have  been  made  for  some 
Western  United  States  faults.  It  is  essential  to  review  previous  geologic  and  seismological  studies 
for  the  region  to  develop  an  understanding  of  the  site’s  tectonic  setting  and  seismic  potentials. 


COMPUTATION  OF  RECURRENCE  PARAMETERS 

The  procedures  discussed  in  this  section  are  equally  applicable  to  regional  analysis  or  fault 
analysis.  The  subset  of  events  assigned  to  the  source  zone  of  interest  are  used  to  calculate  the 
Richter  A  and  B  coefficients,  Equation  4- la  in  Chapter  4.  This  computation  defines  the 
earthquake  recurrence  as  a  line  on  a  semilog  plot.  The  linear  segment  is  bounded  by  a  maximum 
magnitude  determined  as  discussed  above  and  by  a  minimum  magnitude  below  which  the  data 
becomes  nonlinear.  Typically,  the  value  of  B  is  about  -0.9.  The  general  earthquake  recurrence 
is  thus  intially  defined.  However,  as  will  be  shown  in  the  following  sections,  two  important 
elements  are  added  to  geologic  slip  data  and  characteristic  magnitude. 
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GEOLOGIC  SLIP-BASED  RECURRENCE 


Chapter  2  presented  procedures  for  calculating  recurrence  based  on  the  geologic  slip  rate 
data.  Once  the  seismicity  is  estimated  from  the  historical  data,  the  geologic  data  can  be 
compared.  The  procedure  allows  the  user  to  adjust  the  A  and  B  values  replacing  the  historical 
data  values  with  values  based  on  the  longer  span  geologic  data. 

Should  other  studies  be  available,  the  results  of  these  individual  fault  studies  can  be  used 
here  by  adjusting  the  recurrence  parameters. 


CHARACTERISTIC  MAGNITUDE 

As  discussed  in  Chapter  4,  geologic  data  may  show  the  presence  of  history  of  a 
characteristic  event  at  some  average  return  time.  The  seismicity  defined  by  the  historical  data 
fails  to  capture  this  activity,  so  it  is  important  to  include  it  within  the  set  of  events  developed  for 
the  fault.  Once  the  size  of  the  event  and  the  effective  average  return  time  is  defined,  it  becomes 
a  simple  task  to  randomly  add  these  events  into  the  magnitude  list  of  events.  Again,  if  studies 
with  more  advanced  models  are  available  to  define  temporal  distributions,  that  data  can  be  used 
here. 


COMPUTATIONAL  PROCEDURE 

In  Chapter  4,  various  approaches  were  presented  to  determine  the  probability  of 
earthquake  occurrence.  As  shown  above,  various  amounts  of  data  are  required,  some  of  which 
are  beyond  the  scope  of  an  engineering  investigation.  A  new  approach  was  taken  in  the 
formulation  of  a  Monte  Carlo  simulation  procedure.  The  procedure  uses  the  fault  model  and 
regional  model  discussed  earlier  in  this  chapter,  together  with  the  recurrence  procedure.  As 
stated  above,  the  A  and  B  parameters  combined  with  geologic  slip  rate  data  and  characteristic 
magnitude  form  the  basis  for  the  recurrence  function. 

Once  the  recurrence  function  for  a  fault  is  defined,  the  magnitude  distribution  can  be 
computed.  The  process  is  done  for  each  fault  individually.  A  list  of  5,000  events  representing 
the  largest  magnitudes  expected  to  occur  in  50,000  years  is  computed.  For  each  magnitude,  a 
fault  break  length  is  determined  using  data  by  Coppersmith  (1991).  A  random  epicenter  location 
is  selected  along  the  fault.  The  fault  break  is  then  assigned  to  the  random  epicenter.  Various 
distances  are  computed,  such  as  epicentral  distance,  hypocentral  distance,  and  closest  distance 
of  fault  break  to  site.  The  choice  of  distance  depends  on  the  acceleration  attenuation  equation 
chosen  by  the  user. 

Using  the  magnitude  and  separation  distance,  a  site  acceleration  and  standard  deviation 
are  computed.  A  random  acceleration  is  then  determined.  Associated  with  each  acceleration  is 
the  causative  event  and  distance.  The  process  is  repeated  5,000  times  for  each  fault.  The 
random  fault  data  are  then  combined  for  a  total  site  probability  distribution. 

The  procedure  described  above  has  the  advantage  that  historical  data  are  augmented  with 
available  geologic  slip  data.  Where  characteristic  events  are  defined,  they  may  be  easily 
incorporated  at  the  appropriate  return  time.  The  effective  nonlinear  recurrence  function  attempts 
to  capture  the  temporal  characteristics  of  the  data  without  complex  estimates  of  Markov  or 
Bayesian  parameters. 


63 


REFERENCE 


Coppersmith,  K.J.  (1991).  "Seismic  source  characterization  for  engineering  seismic  hazard 
analyses,”  in  Proceedings  of  the  Fourth  International  Conference  on  Seismic  Zonations,  vol  1, 
1991,  pp  1-60. 

The  procedure  has  been  automated  and  the  Users  Manual  is  contained  in  Appendix  A. 
Much  of  the  detailed  computation  is  performed  automatically.  The  user  must  define  coordinates 
of  faults,  maximum  magnitudes,  geologic  slip  data,  and  exposure  period.  The  program  uses  the 
input  data  and  the  epicenter  data  base  to  compute  recurrence  and  probability  of  acceleration  at 
the  site.  An  example  is  given  in  Chapter  8  as  a  typical  case  study. 
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CHAPTER  7 

RESPONSE  SPECTRA  AND  ANALYTICAL  TECHNIQUES 


SPECTRA 

Engineers  have  found  it  useful  to  examine  the  frequency  content  of  earthquake-induced 
ground  motion.  Strong  motion  accelerograms  have  been  analyzed  as  a  means  of  obtaining 
further  insight  into  ground  motion. 

The  first  technique  available  is  Fourier  analysis.  The  Fourier  spectrum  of  an  acceleration 
shows  the  significant  frequency  characteristics  of  the  recorded  motion.  The  Fourier  spectrum 
is  defined  as: 


T 

F(w)  =  f  z(t)  e'"*  dt  C7*1) 

o 


over  the  interval  0<t<T.  The  acceleration  is  zero  outside  the  limits  O  to  T.  The  Fourier 
amplitude  spectrum  is  given  by  the  square  root  of  the  sum  of  the  squares  of  the  real  and 
imaginary  parts  of  F(w).  Associated  with  the  magnitude  F(w)  is  the  phase  angle  p(w),  defined 
as  the  arc  tangent  of  the  imaginary  part  divided  by  the  real  part.  Units  of  the  spectrum  reflect 
the  time  integration;  thus,  for  an  acceleration  in  ft/ sec2,  the  spectrum  magnitude  will  be  in  feet 
per  second  and  the  phase  angle  will  be  in  radians.  The  Fourier  spectrum  magnitude  and  phase 
angle  represent  only  the  input  motion.  Through  convolution  this  may  be  combined  with  transfer 
functions  for  other  elements  (such  as  soil  structure  interaction)  assuming  elastic  behavior.  The 
Fourier  spectrum  magnitude  and  phase  are  a  complete  record  that  is  unique  and  maintains  the 
total  time  history  record.  The  time  history  record  may  be  recreated  by  reverse  transformation. 
Generally,  only  the  magnitude  of  the  Fourier  spectrum  is  shown  in  reference  illustrations. 

In  earthquake  engineering  it  is  important  to  be  able  to  determine  the  magnitude  of 
maximum  response  of  a  structure.  This  has  given  rise  to  the  response  spectra  technique.  A 
single-degree-of-freedom,  spring-mass-damper  system  can  be  analyzed  and  its  time  history  of 
displacement  calculated  to  determine  relative  displacement  between  the  mass  (the  structure)  and 
the  excited  base  (the  ground).  Relative  velocity  and  relative  acceleration  may  also  be  calculated. 
However,  of  primary  interest  for  engineering  applications  is  the  maximum  absolute  values  of 
structure  relative-displacement,  structure  relative- velocity,  and  absolute  structure  acceleration. 
These  values  SD,  SV,  SA  are  functions  of  the  critical  damping.  Plots  of  SD,  SV,  SA  versus  the 
undamped  natural  period  of  vibration  and  for  various  fractions  of  critical  damping  are  called 


response  spectra.  In  typical  engineering  structures  the  damping  is  small  and  for  harmonic 
excitation  the  following  holds: 


SD  =  (— )SV 
2  n 

SA  =  (— )SV 
T 


(7-2) 


For  earthquake-like  excitations  that  are  not  strictly  harmonic  excitation,  an  engineering 
assumption  is  made  that  the  above  is  still  accurate.  The  following  definitions  are  used: 


PSV 


(^L)SD 

T 


(7-3) 


PSA  s  (-^L)2  SD 


where  PSV  is  the  pseudo-relative  velocity  and  PSA  is  the  pseudo-absolute  acceleration.  The  term 
"pseudo"  is  used  to  recognize  the  assumptions  made  concerning  small  damping  and  harmonic 
motion.  Thus,  SD,  PSV,  PSA,  and  T  make  up  a  set  of  data;  knowing  any  two  makes  it  possible 
to  determine  the  other  two.  This  unique  relationship  makes  it  possible  to  plot  response  spectra 
in  tripartite  form.  The  response  spectra  shows  the  response  of  a  single-degree-of-freedom  system 
(structure)  as  a  function  of  damping  of  the  system  and  period  for  the  given  input  acceleration. 

A  comparison  of  Fourier  and  pseudo-velocity  response  spectrum  reveals  that  both  are  of 
the  same  units.  For  an  undamped  oscillator,  a  similar  mathematical  relationship  exists  between 
the  Fourier  amplitude  spectrum  and  exact  relative  velocity  response  spectrum.  (The  Fourier 
spectrum  may  be  viewed  as  the  maximum  velocity  of  the  undamped  oscillator  in  the  free 
vibration  following  the  earthquake;  however,  the  exact  velocity  response  spectrum  is  the 
maximum  velocity  during  both  the  earthquake  and  the  subsequent  free  vibration.) 

The  Fourier  amplitude  spectrum  is  the  quantity  most  frequently  used  by  seismologists  in 
their  investigations  of  the  earthquake  ground  motion.  The  response  spectrum  is  generally  used 
by  structural  engineers  in  the  design  of  structures.  Generally,  the  pseudo-velocity  response 
spectrum  is  the  upper  envelope  of  the  Fourier  spectrum. 

In  summary,  the  most  representative  measure  of  the  driving  ground  motion  is  the  Fourier 
spectrum,  which  may  be  plotted  in  tripartite  form.  The  response  spectrum  gives  the  response 
of  a  single-degree-of-freedom  structure  to  the  ground  motion.  Amplification  of  motion  occurs 
where  resonant  components  of  motion  interact  with  the  structure.  The  response  spectrum  is  most 
widely  used  in  structural  engineering  with  the  modal  analysis  technique.  This  approach 
determines  the  eigenvalues  and  mode  shapes  for  a  number  of  the  modes  of  structure,  and  using 
the  modal  period  determines  the  acceleration  to  be  used  in  computing  an  equivalent  maximum 
static  force  to  be  applied  to  the  structure.  A  response  spectrum  cannot  be  used  directly  as  input 
to  a  dynamic  structural  analysis  to  generate  a  time  history  response  because  the  response 
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spectrum  does  not  contain  a  measure  of  the  duration  of  excitation.  The  Fourier  spectrum  does 
contain  the  full  earthquake  representation  and  may  be  used  in  conjunction  with  transfer  functions 
to  compute  the  structure  response.  This  technique  is  not  used  in  general  structural  engineering 
because  of  the  difficulty  in  determining  the  transfer  function.  The  more  common  structural 
analysis  techniques  that  produce  dynamic  time  histories  of  response  require  input  acceleration 
histories.  Random  analysis  programs  exist  to  generate  time  histories  from  a  given  spectral 
envelope;  however,  the  duration  must  be  specified.  The  resulting  randomly  generated  signal  is 
not  unique,  and  any  number  may  be  generated  from  a  single  spectral  envelope. 

Time  history  records  of  actual  earthquakes  are  available  from  the  California  Institute  of 
Technology  National  Earthquake  Information  Center.  These  may  be  used  when  required  for  time 
history  input  of  ground  motion.  Care  is  needed  in  selection  since  most  records  are  recorded  in 
structures  that  may  influence  the  recording.  Also,  the  region  around  the  site  may  influence  the 
site  response.  The  propagation  of  seismic  waves  is  influenced  by  the  local  geology  and  soil 
conditions.  The  greater  the  extent  of  softer  soils  over  bedrock,  the  less  the  boundary  effects  of 
the  bedrock  will  have  on  the  site.  The  depth  of  soil  overlaying  bedrock  affects  the  period  of 
vibration  of  the  ground.  This  establishes  a  fundamental  soil  frequency  of  particular  importance 
on  soil-structure  interaction.  Further,  this  is  a  factor  in  determining  the  frequencies  of  waves 
filtered  out  by  the  soil,  thus  directly  affecting  the  time  history  record. 


SITE-INDEPENDENT  SPECTRA 

Based  on  studies  of  response  spectra,  Newmark  (1970)  noted  that  response  spectra  could 
be  related  to  peak  ground  acceleration,  velocity,  and  displacement.  From  this  study  it  was 
possible  to  develop  standard  shapes  for  use  in  structural  analysis.  The  standard  ground  motion 
spectra  (defined  as  Ig,  48  inches/second,  36  inches)  could  be  scaled,  based  on  peak  horizontal 
ground  acceleration  expected  at  the  site.  Figure  7-1  gives  amplification  factors  that  could  then 
be  applied  to  estimate  structural  response. 

McGuire  (1977)  also  demonstrates  ratios  of  damped  spectra.  He  found  that  the  5  percent 
damped  pseudo-velocity  spectra  have  value  approximately  70  to  80  percent  of  the  2  percent 
damped  spectra  and  that  the  10  percent  damped  spectra  have  values  of  approximately  55  to  65 
percent  of  the  2  percent  damped  spectra.  These  ratios  are  constant  throughout  frequency  range 
independent  of  the  type  of  earthquake  distance  from  site  and  magnitude.  This  confirms  the 
Newmark-Hall  (1969)  approach,  Newmark,  et  al.  (1973). 

At  a  frequency  of  about  6  cps,  the  amplitude  acceleration  region  line  intersects  a  line 
sloping  down  toward  the  maximum  ground  acceleration  value  and  intersecting  that  line  at  various 
frequencies,  depending  on  the  damping.  The  intersection  is  at  a  frequency  of  about  30  cps  for 
2  percent  damping.  These  lines  are  designated  as  the  acceleration  transition  region  of  the 
spectra.  Finally,  beyond  the  intersection  with  the  maximum  ground  acceleration  line,  the 
response  spectrum  continues  with  the  maximum  ground  acceleration  value  for  higher  frequencies. 

The  spectra  so  determined  can  be  used  as  design  spectra  for  elastic  responses.  The 
spectra  are  completely  described  when  the  maximum  ground  motion  values  are  given  for  the 
three  components  of  ground  motion,  and  the  damping  is  known.  When  only  the  maximum 
ground  acceleration  is  given,  the  values  used  for  maximum  ground  velocity  and  displacement  are 
taken  as  proportional  to  those  in  the  figure,  or  as  scaled  by  the  same  scale  factor  relative  to  the 
maximum  ground  acceleration  compared  with  it. 
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An  assumption  is  made  that  acceleration,  velocity,  and  displacement  values  are 
proportional  to  one  another,  independent  of  magnitude  of  motion.  The  shape  is  thought  correct 
for  sites  on  firm  ground,  soft  rock,  or  competent  sediments.  However,  for  soft  sediments,  the 
velocities  and  displacements  must  be  increased.  Garcia  and  Roesset  (1970)  performed  studies 
comparing  actual  spectra  with  those  estimated  by  the  Newmark-Hall  procedure.  Results  show 
favorable  comparison  indicating  the  utility  of  the  Newmark-Hall  procedure. 

Vertical  spectra  may  be  estimated  by  taking  two-thirds  of  the  horizontal  spectra  when 
fault  movements  are  horizontal  and  by  taking  horizontal  spectra  when  large  vertical  motions  are 
expected. 

Newmark  (1970)  has  studied  the  response  of  elasto-plastic  systems;  one  set  of  results 
shows  the  response  of  an  elasto-plastic  system  to  the  El  Centro  earthquake.  For  low  frequency 
systems,  the  total  displacement  varies  approximately  inversely  with  the  ductility  factor.  In  a  high 
frequency  (stiff)  structure,  the  mass  acceleration  approaches  the  driving  ground  acceleration. 

The  following  generalizations  were  developed  by  Newmark  (1970): 

1 .  For  low  frequency  systems,  the  total  displacement  for  the  inelastic  system  is  the  same 
as  for  an  elastic  system  having  the  same  frequency. 

2.  For  intermediate  frequency  systems,  the  total  energy  absorbed  by  the  spring  is  the 
same  for  the  inelastic  system  as  for  an  elastic  system  having  the  same  frequency. 

3.  For  high  frequency  systems,  the  force  in  the  spring  is  the  same  for  the  inelastic  system 
as  for  an  elastic  system  having  the  same  frequency. 

Newmark  (1970)  has  outlined  a  method  for  selecting  the  response  spectrum  to  use  (Figure 

7-2): 


"The  elastic  spectrum,  designated  by  the  symbol  n  =  1,  for  displacement 
and  acceleration  (D  and  A)  represents  slightly  amplified  values,  corresponding  to 
an  elastic  response  spectrum  for  the  ground  motion  considered.  The  curve 
marked  D  for  n  =  5  is  the  displacement  spectrum  for  a  ductility  factor  of  5,  and 
the  curve  marked  A  for  n  =  5  is  the  acceleration  or  force  spectrum  for  the  same 
conditions.  These  are  drawn  so  as  to  conserve  displacement  on  the  left-hand  side, 
force  on  the  right-hand  side,  and  energy  in  the  central  part.  An  elastic  analysis 
made  for  the  reduced  acceleration  spectrum  therefore  would  correspond  to  the 
ductility  values  derived  for  the  conditions  described.  The  relations  between  the 
various  bounding  lines  in  (Figure  7-2],  for  an  elasto-plastic  resistance  function, 
are  the  same  and  the  acceleration  is  one-fifth  as  much  for  the  elasto-plastic 
spectrum  as  for  the  elastic  spectrum.  Along  a  constant  velocity  line,  the 
displacement  is  five-thirds  as  great  and  the  acceleration  one-third  as  great  for  the 
elasto-plastic  spectrum  compared  with  the  elastic  spectrum.  Finally,  along  a  line 
of  constant  acceleration,  the  displacement  is  five  times  as  great  and  the 
acceleration  value  is  the  same  as  the  value  for  elastic  response. " 

The  elastic  spectrum  discussed  above  may  be  adjusted  to  approximate  inelastic  behavior 
of  a  structure.  The  displacement  region  and  the  velocity  region  are  divided  by  the  structure 
ductility  to  obtain  yield  displacement  D’  and  velocity  V’  (Figure  7-3).  The  acceleration  region 
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(right  side)  is  relocated  by  choosing  it  at  a  level  which  corresponds  to  the  same  energy  absorption 
for  elasto-plastic  behavior  as  for  elastic  for  the  same  period  of  vibration. 

The  extreme  right-hand  portion  of  the  spectrum,  where  the  response  is  governed  by  the 
maximum  ground  acceleration,  remains  at  the  same  acceleration  level  as  for  the  elastic  case  and, 
therefore,  at  a  corresponding  increased  total  displacement  level.  The  frequencies  at  the  comers 
are  kept  at  the  same  values  as  in  the  elastic  spectrum.  The  acceleration  transition  region  of  the 
response  spectrum  is  now  drawn  also  as  a  straight  line  transition  from  the  newly  located 
amplified  acceleration  line  to  the  ground  acceleration  line,  using  the  same  frequency  points  of 
intersection  as  in  the  elastic  response  spectrum. 

In  all  cases,  the  "inelastic  maximum  acceleration"  spectrum  and  the  "inelastic  maximum 
displacement"  spectrum  differ  by  the  factor  im  at  the  same  frequencies.  The  design  spectrum  so 
obtained  is  shown  in  Figure  7-3.  Both  the  maximum  displacement  and  maximum  acceleration 
bounds  are  shown  for  comparison  with  the  elastic  response  spectrum. 

The  solid  line  DVAA0  shows  the  elastic  response  spectrum.  The  heavy  circles  at  the 
intersections  of  the  various  branches  show  the  frequencies  that  remain  constant  in  the  construction 
of  the  inelastic  design  spectrum. 

The  line  D’V’A’A0  shows  the  inelastic  acceleration,  and  the  line  DVA"A0  shows  the 
inelastic  displacement.  These  two  differ  by  a  constant  factor  n  for  the  construction  shown,  but 
A  and  A’  differ  by  the  factor  2/n  -  1,  since  this  is  the  factor  that  corresponds  to  constant 
energy. 

A  study  by  Newmark  and  Riddell  (1979)  investigated  the  response  of  elasto-plastic 
systems  to  numerous  earthquake  records.  Based  on  this  effort,  the  preceding  work  was  found 
to  be  unconservative  for  damping  larger  than  5  percent  and  for  ductilities  greater  than  3. 


SITE-MATCHED  SPECTRA 

The  data  base  of  strong  motion  records  can  be  a  useful  source  of  seismic  data.  Records 
may  be  selected  to  represent  seismologic,  geologic,  and  local  site  conditions.  Selection  is 
complicated  by  a  number  of  factors.  Ideally,  the  records  should  be  selected  to  match  source-site 
transmission  path,  source  mechanism,  and  local  site  conditions.  These  are  not  readily 
quantifiable.  Thus,  reliance  is  made  on  earthquake  magnitude  site  acceleration  level,  site 
classification,  and  duration  of  motion.  Judgment  is  an  important  factor  in  selecting  and  scaling 
records. 

Appendix  B,  based  on  Ferritto  (1992),  is  used  to  compute  optimized  earthquake  time 
histories  and  response  spectra.  The  program  has  the  data  base  of  about  1 ,000  records  provided 
by  the  National  Oceanographic  and  Atmospheric  Administration.  The  user  may  select  specific 
records  and  obtain  time  histories  and  spectra  or  may  specify  a  ground  acceleration  level,  site 
distance,  and  magnitude  and  the  program  will  search  the  data  base  and  provide  the  user  with  a 
list  of  the  closest  matching  records.  The  user  may  then  combine  a  number  of  spectra  and  obtain 
average,  average  plus  one  standard  deviation,  and  envelope  spectra. 

It  is  suggested  that  site-matched  groups  of  spectra  be  used  to  develop  the  mean  and  mean- 
plus- 1 -standard-deviation  spectra.  These  should  be  compared  with  standard  spectral  shapes  and 
typical  results  for  soft,  intermediate,  or  rock  sites  to  denote  regions  where  the  spectra  may  be 
deficient.  This  is  particularly  important  for  Eastern  sites  since  the  spectra  are  recorded  in  the 
West.  Significant  variations  in  attenuation  have  been  noted  between  Western  and  Eastern  ground 
motion. 


69 


The  same  nominal  ground  motion  level  can  be  produced  by  different  magnitude 
earthquakes  at  different  distances.  As  an  example  of  this,  a  nominal  ground  motion  of  0. 15  g 
was  caused  by  the  following  earthquakes: 


Earthquake 

Magnitude 

Separation 

Distance 

(miles) 

A 

5 

14 

B 

6 

18 

C 

7 

30 

Using  response  spectra  programs,  the  closest  matching  ten  spectral  records  for  each  earthquake 
(intermediate  soil  class)  were  used  to  generate  average  and  maximum  spectra  (Figure  7-4).  Each 
spectrum  represented  the  characteristics  of  the  magnitude  and  separation  distance. 

In  the  response  spectra  program,  the  closest  matching  ten  spectral  records  for  each 
earthquake  (intermediate  soil  class)  were  used  to  generate  average  and  maximum  spectra  (Figure 
7-4).  Each  spectrum  represented  the  characteristics  of  the  magnitude  and  separation  distance. 

The  response  spectra  data  base  and  accelerograms  have  been  categorized  by  soil  site 
conditions:  alluvium,  intermediate,  or  rock.  It  was  of  interest  to  evaluate  this  effect  on  the 
response  of  the  structure.  Ten  spectra  of  each  site  type  for  a  magnitude  6.6  event  at  20  miles 
having  a  0.15  g  nominal  ground  motion  were  combined  to  produce  average  and  maximum 
spectra  (Figure  7-5). 

Two  alternative  procedures  used  to  determine  site  specific  ground  motion  are  as  follows. 
Surface  Motion 

This  technique  utilizes  attenuation  relationships  based  on  surface  motion.  The  computed 
motion  is  then  used  as  a  scaling  value  for  the  response  spectra.  The  specific  response  spectrum 
may  be  based  on  a  group  selected  for  similar  site  properties  or  a  spectral  shape  determined  by 
researchers  to  be  applicable  to  specific  site  conditions. 

Bedrock  Motion 

This  technique  utilizes  an  attenuation  relationship  based  on  bedrock  motion.  The  motion 
may  be  brought  to  the  surface  either  from  empirical  data  or  by  use  of  wave  propagation 
computation  programs. 

An  automated-analysis  technique,  widely  used  today  for  treating  horizontal  soil  layers, 
has  been  developed  by  Schnabel,  Lysmer,  and  Seed  (1972),  based  on  the  one-dimensional  wave 
propagation  method.  This  program,  SHAKE,  can  compute  the  responses  for  a  given  horizontal 
earthquake  acceleration  specified  anywhere  in  the  system.  The  analysis  incorporates  nonlinear 
soil  behavior,  the  effect  of  the  elasticity  of  the  base  rock,  and  variable  damping.  It  computes 
the  responses  in  a  system  of  homogeneous  viscoelastic  layers  of  infinite  horizontal  extent,  subject 
to  vertically  traveling  shear  waves.  The  program  is  based  on  the  continuous  solution  of  the  wave 
equation  adapted  for  use  with  transient  motions  through  the  Fast  Fourier  Transform  algorithm. 
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Equivalent  linear  soil  properties  are  obtained  by  an  iterative  procedure  for  values  of  modulus  and 
damping  compatible  with  the  effective  strains  in  each  layer.  The  following  assumptions  are 
made: 

1.  The  soil  layers  extend  infinitely  in  the  horizontal  direction. 

2.  The  layers  are  completely  defined  by  shear  modulus,  critical-damping  ratio,  density, 
and  thickness. 

3.  The  soil  values  are  independent  of  frequency. 

4.  Only  vertically  propagating,  horizontal  shear  waves  are  considered. 

The  soil  model  is  similar  to  that  developed  by  Seed  and  Idriss  (1970),  using  data  similar 
to  Hardin  and  Dmevich  (1970).  The  absolute  range  of  soil  parameter  variation  may  be  stipulated 
by  merely  inputting  factors  whose  numerical  values  may  be  derived  from  simple  soil  strength 
properties.  These  strength  properties  may  be  the  undrained  shear  strength  of  a  clay  or  the 
relative  density  for  sands.  The  program  requires  the  definition  of  the  soil  profile  down  to 
bedrock  (defined  as  seismic  velocity  2,500  ft/sec)  as  well  as  an  earthquake  time  history  record 
in  digital  form. 

The  motion  used  as  a  basis  for  the  analysis  can  be  given  in  any  layer  in  the  system,  and 
new  motions  can  be  computed  in  any  other  layer.  Maximum  stresses  and  strains,  as  well  as  time 
histories,  may  be  obtained  in  the  middle  of  each  layer.  Response  spectra  may  be  obtained  and 
amplification  spectra  determined. 

Using  the  site  soil  properties,  a  one-dimensional  wave  propagation  analysis  was 
performed.  The  time  histories  of  the  records  selected  were  used  in  the  one-dimensional  analysis. 
The  ground  motion  was  applied  to  the  surface.  The  average  shear  velocity  of  the  site  was  125 
m/sec  with  a  natural  period  of  1.7  seconds.  The  site  responded  as  anticipated,  as  a  deep,  soft 
site  showing  attenuation  of  the  acceleration.  The  surface  motion  (input  at  0.5  g)  was  increased 
to  about  1.0  g  at  a  depth  of  57  meters.  The  site  attenuated  bedrock  motion  by  a  factor  of  about 
two.  This  is  in  general  agreement  with  empirical  data. 

Figure  7-6  shows  input  surface  response  spectrum  and  computed  bedrock  response 
spectrum.  The  bedrock  spectrum  has  its  peak  value  at  a  period  of  0. 13  second.  The  shape  of 
this  spectrum,  location  of  peak  value,  and  relative  magnitudes  are  very  much  in  agreement  with 
standard  spectral  data  based  on  records  recorded  in  rock  sites.  The  soil  layer  responds  as  a  filter. 
Thus,  it  is  also  possible  to  generate  surface  motion  using  a  scaled  rock-site  spectrum  as  input. 
The  magnitude  of  the  spectra  would  be  based  on  attentuation  relations  developed  for  rock. 
Figure  7-7  gives  a  flow  chart  of  two  approaches  for  determining  an  average  surface  response 
spectrum. 

The  approach  of  using  surface- recorded  scaled,  response  spectra  matched  to  the  site 
conditions  is  thought  to  be  a  better  representation  of  actual  conditions  than  the  alternative  of 
attempting  to  compute  ground  motion  propagation  from  bedrock  to  the  surface.  The  limitations 
in  the  accuracy  of  the  attenuation  equation  show  no  statistical  difference  between  peak 
accelerations  recorded  on  rock  and  those  on  soil  at  comparable  distances.  Thus,  the  problem  of 
what  level  of  motion  to  input  must  be  based  on  uncertain  data.  Motion  must  be  artificially 
brought  to  bedrock  by  deconvolution.  Unfortunately,  motions  are  usually  recorded  on  the 
surface  and  not  at  bedrock  depth.  Without  at-depth  experimental  records,  one-dimensional  wave 
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propagation  calculations,  although  very  useful,  may  have  error.  Spectral  shapes  from  such 
calculations  cannot  be  used  with  absolute  certainty.  The  extent  of  site  amplification  is  a 
significant  parameter.  However,  it  cannot  be  computed  from  wave  propagation  analysis  in 
absolute  certainty.  It  should  be  looked  at  as  relating  relative  soil  behavior.  Uncertainty  is 
introduced  by  the  choice  of  material  properties  used  to  characterize  the  site.  The  assumptions 
made  in  one-dimensional  analysis  are  perhaps  more  hidden  and,  thus,  create  a  greater  confidence 
in  the  results. 

It  is  important  to  consider  the  major  assumption  made  in  wave  propagation  analysis:  that 
vertically  propagating  shear  waves  travel  through  horizontal  layers.  For  sites  close  to  the  fault, 
the  inclined  nature  of  the  fault  and  the  close  horizontal  proximity  to  the  energy  source  must  be 
considered.  The  energy  released,  which  is  composed  of  surface  and  body  waves,  cannot  be 
represented  by  a  simplified  one-dimensional  model.  Thus,  it  is  questionable  whether  any 
attentuation  of  motion  would  actually  occur.  The  one-dimensional  model  is  best  suited  for  sites 
at  distances  from  the  source  where  propagation  is  essentially  through  the  more  competent 
subsurface  (bedrock)  layers  refracting  to  the  surface.  This  site  may  indeed  show  attenuation  to 
motion  originating  from  distant  sources.  However,  little  is  known  about  close-in  behavior;  there 
is  not  enough  known  to  justify  a  reduction  in  ground  motion  without  loss  of  confidence  in  the 
results. 


REFERENCES 

Ferritto,  J.M.  (1992).  Optimized  earthquake  time-history  and  response  spectra.  Naval  Civil 
Engineering  Laboratory,  User  Guide  UG-0025.  Port  Hueneme,  CA,  Oct  1992. 

Garcia,  F.,  and  J.M.  Roesset  (1970).  Influence  of  damping  on  response  spectra,  Massachusetts 
Institute  of  Technology,  Department  of  Civil  Engineering,  R70-4.  Cambridge,  MA,  Jan  1970. 

Hardin,  B.,  and  V.  Dmevich  (1970).  Shear  modulus  and  damping  in  soils,  University  of 
Kentucky,  College  of  Engineering,  Soil  Mechanics  Series  Technical  Report  UKY  26-70-CE2. 
Lexington,  KY,  Jul  1970. 

McGuire,  R.  (1977).  "Seismic  design  spectra  and  mapping  procedures,"  Earthquake  Engineering 
and  Structural  Dynamics,  vol  5,  New  York,  NY,  Wiley  Publications,  1977,  pp  211-234. 

Newmark,  N.M.  (1970).  Current  trends  in  the  seismic  analysis  and  design  of  high-rise  structures 
in  earthquake  engineering,  Chapter  16,  R.L.  Wiegel,  editor.  Englewood  Cliffs,  NJ,  Prentice- 
Hall,  1970. 

Newmark,  N.M.,  and  R.  Riddell  (1979).  "A  statistical  study  of  inelastic  response  spectra,"  in 
Proceedings  of  the  Second  U.S.  National  Conference  on  Earthquake  Engineering,  Stanford 
University,  Stanford,  CA,  Aug  1979. 

Newmark,  N.M.,  J.A.  Blume,  and  K.K.  Kapur  (1973).  "Seismic  design  spectra  for  nuclear 
power  plants,"  Journal  of  the  Power  Division,  ASCE,  vol  99,  no.  P02,  Nov  1973. 


72 


Newmark,  N.M.,  and  W.J.  Hall  (1969).  "Seismic  design  criteria  for  nuclear  reactor  facilities," 
in  Proceedings  of  the  Fourth  World  Conference  on  Earthquake  Engineering,  Santiago,  Chile, 
1969. 

Schnabel,  P.B.,  J.  Lysmer,  and  H.B.  Seed  (1972).  SHAKE,  a  computer  program  for  earthquake 
response  analysis  of  horizontally  layered  sites,  University  of  California,  Earthquake  Engineering 
Research  Center,  EERC  Report  No.  72-12.  Berkeley,  CA,  Dec  1972. 

Seed,  H.B.,  and  I.M.  Idriss  (1970).  Soil  moduli  and  damping  factors  for  dynamic  response 
analysis,  University  of  California,  Earthquake  Engineering  Research  Center,  EERC  Report  No. 
70-10.  Berkeley,  CA,  Dec  1970. 


73 


Velocity,  in /toe. 
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(a)  Relative  values  of  spectrum  amplification  factors 

(from  "Seismic  design  spectra  for  nuclear  power  plants," 
by  N.  M.  Newmark,  J.  A.  Blume,  and  K.  K.  Kapur,  in 
Journal  of  the  Power  Division,  ASCE,  vol  99,  no.  P02, 

Nov  1973,  table  5). 


Figure  7-1  Design  Spectra  normalized  to  1.0  g 
(from  Applied  Technology  Council,  1974) 
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Figure  7-2.  Comparison  of  elastic  and  elasto-plastic  spectra. 
(Natan  M  Newmark,  "Current  Trends  in  the  Seismic  Analysis  and 
Design  of  High  Rise  Sructures",  in  Earthquake  Engineering,  edited 
by  Robert  L.  Weigel,  copy  right  1970,  pp  419.  reprinted  by 
permission  of  Prentice  Hall,  Inc.  Englewood  Cliffs,  New  Jersey. 
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Figure  7-3  Design  spectra  for  inelastic  analysis. 
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Figure  7-4a.  Magnitude  5  at  14  miles 
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Figure  7-4b.  Magnitude  6  at  17  miles. 
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Figure  7-4c.  Magnitude  7  at  30  miles. 
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Figure  7-5a.  Magnitude  6.6  at  20  miles,  O.lSg,  alluvium  site. 
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Figure  7-5b.  Magnitude  6.6  at  20  miles,  0.15g,  intermediate  site. 
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.  Magnitude  6.6  at  20  miles,  0.15g,  rock  site. 
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Figure  7-7.  Flow  chart  of  analysis. 


CHAPTER  8 

THE  SAN  DIEGO  AREA  -  AN  EXAMPLE  CASE  STUDY 


SEISMICITY 

The  seismicity  and  regional  geologic  structure  of  the  San  Diego  area  can  be  interpreted 
in  light  of  current  plate  tectonic  theory.  California  is  believed  to  lie  on  the  junction  of  two 
relatively  rigid  plates  of  the  earth’s  crust  that  respond  to  movement  of  subcrustal  material.  The 
main  evidence  of  this  juncture  is  the  San  Andreas  fault.  These  same  forces  that  tend  to  move 
the  portion  of  California  on  the  westerly  side  of  the  San  Andreas  fault  northward  have  resulted 
in  the  formation  of  other  faults,  such  as  the  San  Jacinto,  Whittier-Elsinore  and  Newport- 
Inglewood  faults. 

Distant  faults  that  must  be  considered  significant  to  the  site  region  include  the  Elsinore 
and  San  Jacinto  fault  zones  to  the  northeast  and  the  San  Clemente  fault  zone  to  the  west.  Local 
faults  include  the  Rose  Canyon  and  La  Nacion.  The  San  Andreas  fault  zone  is  not  considered 
very  significant  because  of  its  great  distance  from  the  study  area. 

The  San  Diego  Bay  contains  cretaceous,  tertiary,  and  quaternary  strata,  which  is  generally 
flat  but  locally  folded  and  cut  by  normal  and  right  lateral  faults.  This  area  is  called  the  Rose 
Canyon  zone  (Lamar,  et  al.,  1973).  A  bottom  survey  of  the  bay  revealed  numerous  faults  which 
were  difficult  to  correlate.  The  quaternary  deformations  observed  along  the  Rose  Canyon  fault 
zone  attest  to  the  tectonic  importance  of  the  zone.  Although  no  major  earthquakes  have  occurred 
near  San  Diego  recently,  several  earthquakes  of  about  magnitude  3.5  have  been  recorded  during 
the  past  41  years.  Eleven  took  place  near  the  Rose  Canyon  fault.  The  magnitude  3.5  earthquake 
is  associated  with  a  fault  rupture  length  of  1  km.  The  geologic  structure  of  this  area  shows 
evidence  of  previous  movement.  Surface  traces  of  more  than  24  km  in  length  and  vertical 
separation  of  hundreds  of  feet  are  visible.  Table  8-1  shows  the  key  faults  and  the  maximum 
credible  earthquake. 

San  Jacinto  Fault 

The  San  Jacinto  fault  system  extends  from  its  junction  with  the  San  Andreas  southeast  of 
Palmdale  to  the  Colorado  River  Delta.  Geodetic  data  indicate  an  average  slip  rate  of  0.3  cm/yr. 
Seventeen  large  earthquakes  have  occurred  since  1890  along  the  290-km  long  fault.  The 
magnitudes  determined  were  in  the  range  of  5.7  to  7.1. 

This  is  one  of  the  most  active  faults.  In  the  33-year  period  from  1890  to  1923,  the 
northern  portion  of  this  fault  system  averaged  an  event  each  5.5  years.  Large  earthquake  activity 
in  the  northern  half  of  the  fault  system  has  lapsed  during  the  past  52  years.  Movement  in  faults 
in  the  Imperial  Valley  caused  an  earthquake  in  1915  and  again  in  1940  (25-year  span).  The  last 
event  was  35  years  ago,  and  it  is  thought  that  significant  strain  has  not  been  released  since  that 
event. 
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Whittier-Elsinore  Fault 


This  fault  system  is  composed  of  the  Elsinore  and  Whittier  fault  zones,  Agua  Caliente 
fault,  and  Earthquake  Valley  fault.  Five  recent  earthquakes  of  unknown  magnitudes  have 
occurred  on  this  fault,  the  last  one  in  1935.  No  historical  data  exist  to  construct  a  recurrence 
relationship.  A  slip  rate  of  0.08  cm/yr  was  determined  and  used  to  calculate  a  recurrence 
interval.  It  is  believed  that  sufficient  elastic  strain  to  produce  a  magnitude  6  or  greater 
earthquake  has  accumulated  along  the  fault  in  recorded  historic  time  (several  hundred  years). 

San  Clemente  Fault 

This  fault,  with  a  verified  length  of  176  km,  extends  from  the  eastern  side  of  San 
Clemente  Island  to  the  Cabo  Colonet  area  of  Baja  California,  Mexico.  A  magnitude  5.9 
earthquake  occurred  off  the  southeast  tip  of  San  Clemente  Island  in  1951.  The  maximum 
credible  event  for  a  fault  of  this  length  is  7.7  on  the  Richter  scale.  A  significant  consequence 
of  an  earthquake  on  this  fault  is  the  possible  production  of  a  tsunami  or  seismic  sea  wave,  but 
such  is  not  likely  with  magnitudes  less  than  6.3.  Seven  percent  of  southern  California 
earthquakes  have  submarine  epicenters,  yet  only  two  or  three  locally-generated  tsunamis  are 
known  to  have  occurred  since  1800,  and  none  in  the  San  Diego  area. 

Rose  Canyon  Fault 

The  Rose  Canyon  fault  zone  forms  a  belt  of  fractures  about  a  mile  wide.  The  zone  on 
shore  can  be  traced  southwestward  for  a  distant*  of  more  than  16  km  and  then  projects  under 
San  Diego  Bay  and  continues  to  the  Mexican  border  and  possibly  beyond.  An  investigation  of 
the  bay  (Moore,  1972)  revealed  many  faults.  North  of  La  Jolla  an  offshore  extension  of  the  fault 
exists,  suggesting  that  Rose  Canyon  is  part  of  a  much  larger  northwest  trending  zone  of 
deformation  that  extends  at  least  240  km  from  Santa  Monica,  California  to  Baja,  California  and 
includes  the  New-Inglewood  zone. 1  Geologic  evidence  suggests  that  the  most  recent  movement 
was  less  than  500,000  years  ago.  Fault  displacements  as  recently  as  early  Holocene  time  (10,000 
years  ago)  cannot  be  precluded;  Moore  (1972)  cites  evidence  of  faulting  through  Pleistocene 
deposits.  No  large  earthquakes  have  been  associated  with  the  Rose  Canyon  fault  during  historic 
time.  During  1964,  however,  three  earthquakes  in  the  magnitude  range  of  3.5  to  3.7  were  felt 
in  the  vicinity  of  San  Diego  Bay.  Uplift  has  been  noted  near  the  fault  in  La  Jolla  in  sediments 
that  overlie  the  Linda  Vista  formation  and  are  considered  to  be  late  Pleistocene  (1,000,000  years 
old). 

La  Nacion  Fault 

This  fault  extends  for  24  km  southward  from  La  Mesa.  The  fault  dips  60  to  70  degrees 
toward  the  west  and  consists  of  two  or  more  branches;  several  are  tens  of  feet  apart.  Offsets  of 
Holocene  deposits  along  the  fault  have  been  noted  but  not  confirmed.  There  is  definite  evidence 
that  Linda  Vista  formations  (50,000  years)  have  been  displaced.  Evidence  (McEuen  and 
Pinckney,  1972)  suggests  a  change  in  the  history  of  the  San  Diego  area,  previously  thought  to 


Possible  locus  of  the  1933  Long  Beach  earthquake,  magnitude  6.3. 


86 


be  stable.  Test  borings  approximately  30  meters  deep  show  open  faults  in  sedimentary  rock, 
suggesting  the  area  is  in  tension.  McEuen  and  Pinckney  (1972)  also  note  that  1  km  north  of  the 
Otay  Valley  the  La  Nacion  fault  offsets  Pliocene  formations,  late  Pleistocene  terrace  deposits, 
and  Holocene  alluvium  (dating  10,090  +.  190  years). 


PROBABILITY  ANALYSIS 

This  study  is  intended  as  a  demonstration  of  the  procedure.  The  bounds  of  the  study  area 
are  117.0  to  119.0  W  longitude,  34.5  to  31.0  N  latitude.  The  coordinates  of  the  site  are 
1 17. 125N,  32.708N.  A  set  of  historical  data  was  prepared  for  the  site  containing  over  6,000 
events  with  magnitudes  of  3  or  greater. 

Figure  8-1  shows  the  region  of  interest  with  the  epicenters  plotted.  Figure  8-2  shows  a 
similar  plot  with  only  the  faults  shown.  Figure  8-3  shows  the  computed  recurrence  for  the  La 
Nacion  fault.  All  other  faults  utilized  the  default  data  in  the  program.  Figure  8-4  shows  the 
individual  contributions  of  each  fault  on  the  site  probability.  Figure  8-5  shows  the  total 
probability  of  not  exceeding  the  acceleration  for  a  50-year  exposure.  Figure  8-6  shows  a 
generalized  site-independent  spectra  normalized  to  0.33  g. 


RESPONSE  SPECTRA 

The  site  has  been  characterized  as  an  intermediate  site.  The  closest  ten  records  matching 
a  6.5  event  16  km  from  this  site  producing  0.33  g  acceleration  were  used  to  produce  the  response 
spectra  shown  in  Figures  8-6  and  8-7.  These  are  useful  in  structural  design. 
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Table  8-1 

Fault  Systems  of  Interest  to  NAS  North  Island 


Fault 

Maximum 

Credible 

Magnitude 

Coyote  Creek 

7.0 

Elsinore 

7.5 

Imperial 

7.0 

La  Nacion 

6.8 

Malibu 

7.5 

Newport-Inglewood 

7.0 

Palos  Verdes 

7.0 

Pinto  Mountain 

7.5 

Raymond  Hills 

7.5 

Rose  Canyon 

7.1 

San  Clemente 

7.7 

San  Gabriel 

7.7 

San  Jacinto 

7.5 

Santa  Susana 

6.5 

Sierra  Madre 

6.5 

South  San  Andreas 

7.5 

Superstition  Mountain 

7.0 
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Figure  8-1.  San  Diego  region  showing  epicenters. 
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Figure  8-2.  San  Diego  region  showing  faults. 
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Figure  8*7.  Site  dependent  spectra, 
10  closest  matching  records. 
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CHAPTER  9 
SUMMARY 


NAVFAC  instructions  suggest  that,  in  conducting  site  seismicity  studies  for  key  facilities, 
ground  motion  and  response  spectra  be  defined  on  a  probabilistic  basis.  However,  no  procedure 
was  specified  by  NAVFAC  to  accomplish  this.  This  project  was  initiated  to  provide  a  procedure 
and  the  required  software. 

An  automated  procedure  has  been  developed  to  perform  a  seismic  analysis  using  available 
historic  data  and  geologic  data.  The  objective  of  the  seismicity  study  is  to  determine  the 
probability  of  occurrence  of  acceleration  at  the  site.  To  do  this,  site  coordinates  and  the  study 
bounds  are  specified  in  terms  of  latitude  and  longitude.  A  regional  study  is  first  performed  in 
which  all  of  the  historic  epicenters  are  used  with  an  attenuation  relationship  to  compute  site 
acceleration  for  all  historic  earthquakes.  A  regression  analysis  is  performed  to  obtain  regional 
recurrence  coefficients,  and  a  map  of  epicenters  is  plotted.  The  regional  recurrence  can  be  used 
to  compute  the  probability  of  site  acceleration  for  randomly  located  events  in  the  study  area. 
Such  a  condition  is  used  when  individual  faults  are  not  known  well  enough  to  be  specified. 

Where  individual  fault  areas  can  be  specified,  individual  subsets  of  the  historic  data  are 
used  in  conjunction  with  geologic  data  to  determine  fault  recurrence  coefficients;  these  are  used 
to  compute  the  probability  of  site  acceleration  from  individual  fault  sources.  The  total  risk  is 
determined  for  all  faults  specified.  Confidence  bounds  are  given  on  the  site  acceleration  as  a 
function  of  probability  of  not  being  exceeded. 

The  structural  design  engineer  may  use  either  reponse  spectra  or  time  history  techniques 
in  the  analysis  of  a  structure.  The  data  base  of  recorded  accelerograms  has  been  obtained  and 
a  program  was  prepared  to  search  the  record  of  accelerograms,  given  a  desired  magnitude  event, 
epicenter-site  distance,  acceleration  level,  and  soil  condition,  to  determine  the  closest  matching 
records.  The  program  takes  selected  response  spectra,  and  scales  them,  and  then  computes  the 
mean  and  standard  deviation  spectra  and  the  maximum  envelope  spectrum.  The  spectra  are 
plotted  either  in  tripartite  form  or  in  semilog  form.  The  program  also  is  able  to  scale,  plot,  and 
create  files  of  time  history  accelerograms  for  use  as  input  to  dynamic  finite  element  programs. 

Case  studies  were  conducted  to  evaluate  the  procedure.  Results  compare  favorably  with 
results  by  others.  v 
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User's  Guide 


Seismic  Hazard  Analysis 
Appendix  A 


INTRODUCTION  -  what  is  a  Seismicity  Study 


The  objective  of  a  seismicity  study  is  to  quantify  the  level 
and  characteristics  of  the  earthquake  ground  motion  which  pose  a 
risk  to  a  site  of  interest.  The  seismicity  study  will  produce  a 
probability  distribution  of  expected  site  acceleration  for  a 
given  exposure  period  and  also  give  an  indication  of  the  frequen¬ 
cy  content  of  that  motion.  The  approach  taken  in  this  work  is  to 
use  the  historical  epicenter  data  base  in  conjunction  with  geo¬ 
logic  data  where  available  to  best  estimate  the  earthquake  recur¬ 
rence  of  a  region  or  fault.  This  recurrence  relationship  is  used 
to  determine  the  regional  or  fault  contributions  to  the  overall 
site  acceleration  level.  This  becomes  the  basis  for  definition 
of  response  spectra  suitable  for  use  in  structural  design  and 
analysis.  The  procedure  utilizes  three  parts  to  accomplish  the 
study.  The  first  part  creates  a  subset  of  earthquakes  from  the 
general  data  base  and  plots  all  events  within  a  specified  region. 
The  second  part  utilizes  the  epicenter  data  to  perform  a  regional 
analysis  determining  the  magnitude  recurrence  relationship  for 
the  region  which  may  be  adjusted  for  geologic  data  where  known. 
Additionally  the  program  computes  the  probability  of  acceleration 
at  the  site  location  and  gives  plots  of  recurrence  and  accelera¬ 
tion  data.  The  third  part  analyzes  individual  faults.  It  deter¬ 
mines  fault  magnitude  recurrence  and  probability  of  acceleration 
at  the  site  from  an  event  on  each  fault  specified.  Geologic  data 
may  be  used  to  augment  historical  epicenter  date.  Each  part  will 
be  discussed  in  detail  below. 


GETTING  STARTED 

System  Requirements 

The  Program  is  designed  to  run  on  standard  desk  top  personal 
computers  using  the  MS  DOS  operating  system  version  3  or  higher. 
The  following  are  required: 

MS  DOS  version  3  or  higher 
640  k  system  memory 
Hard  disk 

80287  math  coprocessor 


optional  devices 


plotter 

printer 


The  following  plotters  are  supported: 


EPSON  FX80  printer  or  compatible  (used  as  plotter) 
Hewlett  Packard  Laserjet  printer  or  compatible 
Hewlett  Packard  Plotters 
Houston  Instruments  Plotters 
Tecktronix  4025 
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Loading  The  Program 

Place  the  first  PROGRAM  DISK  in  Drive  A:  switch  to  Drive  A 
then  type: 


A>  INSTALL  <cr> 

This  will  create  a  directory  on  the  C  hard  disk  named  SEISMIC  and 
load  all  of  the  programs  into  that  directory.  When  complete, 
place  the  next  disk  in  Drive  A  and  type  INSTALL  as  above  to 
continue  loading  the  program. 

Place  the  first  DATA  DISK  in  Drive  A:  type  INSTALL  as  above. 
This  will  create  a  directory  named  EPICENTR  and  load  the  data 
files.  Place  the  remaining  DATA  DISKS  in  sequence  in  Drive  A  and 
type  INSTALL  as  above  to  load  all  data. 

Place  the  SAMPLE  DATA  DISK  in  Drive  A:  type  INSTALL  as  above. 
This  will  create  a  directory  named  OUTPUT  and  load  the  data 
files. 

See  the  ADVANCED  USER  section  to  deviate  from  these 
locations  for  program  and  data. 

Configuring  The  Program 

The  program  supports  a  number  of  graphic  output  devices;  the 
user  must  specify  which  device  is  being  used  and  how  it  is 
connected.  He  must  also  tell  the  program  where  the  data  files 
are  being  stored.  To  begin  the  configuration  switch  to  the 
program  directory  and  start  the  program  by  typing  the  following: 

C: 

CD\SEISMIC 

SEISMIC 

The  following  screen  will  appear: 


1  OPTIONS 


Directory 

Configure 


-in  Mu' 


USE 


ARROW  KEYS  THEN  PRESS  EH  TEH 
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The  program  has  five  choices  which  form  the  opening  menu: 
Options,  Earthquake  Selection,  Regional  Study,  Fault  Study,  and 
Exit.  Use  the  LEFT  and  RIGHT  ARROW  keys  which  are  on  the  number 
pad  keys  with  the  4  and  6  to  move  among  the  choices.  Make  sure 
the  NUM  LOCK  key  is  off  to  permit  the  arrow  keys  to  function. 
When  the  OPTIONS  choice  is  selected  a  window  opens  giving  two 
choices:  Directory,  and  CONFIGURE.  Use  the  down  arrow  key  to 
choose  the  CONFIGURE  choice  then  press  ENTER  (or  RETURN) .  The 
following  screen  appears: 


Enter  Plotter  DEVICE  1M18ER 
see  User's  lamia 1  e.j.  LPT1  =  1 

1 

Enter  Plotter  HODEL  NUUBER 

see  User’s  Nanual  e.g.  HP  Lasejet  =60 

60 

Enter  Drive  and  Directory  for  EPICENTER  files 
example  CAEQUAKES 

FAEQUAKES 

Enter  Drive  and  Directory  for  OUTPUT  files 
example  CAOUTPUT 

FASEISHIC 

Any  Changes?  Y  /  N 


The  DEVICE  NUMBER  refers  to  the  port  to  which  the  hard  copy 
plotting  device  is  connected;  see  Table  1  for  configuration 
options  and  check  the  manual  for  the  hard  copy  plotting  device. 
Enter  the  DEVICE  NUMBER  then  press  ENTER  or  DOWN  ARROW. 

The  MODEL  NUMBER  can  be  obtained  from  Table  2a  for  the 
devices  supported  by  the  program.  Laser  printers  produce  high 
quality  plots  rapidly  and  are  recommended.  Table  2b  gives  a 
matrix  MODEL  NUMBERS  for  various  compatible  printers  which  can 
be  used  to  obtain  plots.  Table  3  gives  the  recommended 
configuration  for  specific  devices.  Enter  the  MODEL  NUMBER  then 
press  ENTER  or  DOWN  ARROW. 

Type  the  directory  name  where  the  epicenter  data  files  are 
located.  If  you  used  the  INSTALL  routine  to  accept  the  default 
directory  creation  then  type: 

C:\EPICENTR 

Note  the  back-slash  (\)  key  and  the  spelling  of  EPICENTR,  8 
letters  without  the  E.  If  you  chose  to  locate  the  programs  else 
where  enter  the  following: 

Drive  letter : \Directory 
example  D:\EPIC 
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Tab  la  1.  Davica  Nuabar 


Dow lea 

Output  Davies 

Nuabar 

Porollol  Port 

0 

PRN:  (PUN 

:  Is  oqulvslout  to  LPT1: 

) 

1 

LPT1: 

2 

LPT2: 

3 

LPT2: 

-coosola* 

99 

-aarial  porta* 

davica 

bond 

parity 

Odata 

- i 

fstop 

rsto 

bits 

bits 

300 

C0K1: 

300 

N 

8 

1 

301 

corn: 

300 

0 

7 

302 

cam: 

300 

E 

7 

1 

1200 

C0H1: 

1200 

H 

8 

1 

1201 

cam: 

1200 

0 

7 

1 

1202 

cam: 

1200 

E 

7 

1 

4000 

cam: 

4800 

N 

8 

1 

4801 

CONI: 

4800 

0 

7 

1 

4802 

C0M1: 

4800 

E 

7 

1 

9600 

C0H1: 

9800 

N 

8 

1 

9801 

CONI: 

9800 

0 

7 

1 

9802 

CONI: 

9800 

E 

7 

1 

parity: 

N-Nons 

E-Evan 

O-Odd 

COM2: -Add  SO  to  valua  for  C0M1: 
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Tabla  2a.  Nodal  Nuabar 


Nodal 

Nuabar 

Printar-Plottar-Scrsan 

Davlca  Idantlf ication 

0 

Epson  FX-80  Printer,  single  density. 

1 

Epson  FX-80  Prlntar,  donbla  dansity. 

2 

Epson  FX-80  Prlntar,  donbla  spaad,  dual 
dansity. 

3 

Epson  FX-80  Prlntar,  quad  dansity 

4 

Epson  FX-80  Prlntar,  CXT  Graphics  I. 

5 

Epson  FX-80  Prlntar,  plottar  graphics. 

6 

Epson  FX-80  Prlntar,  CXT  Graphics  II. 

10 

Epson  FX-100  Prlntar,  singla  dansity. 

11 

Epson  FX-100  Prlntar,  donbla  dansity. 

12 

Epson  FX-100,  donbla  spaad,  dnal  dansity. 

13 

Epson  FX-100  Prlntar,  quad  dansity. 

14 

Epson  FX-100  Prlntar,  CXT  Graphics  I. 

IS 

Epson  FX-100  Prlntar,  plottar  graphics. 

16 

Epson  FX-100  Prlntar,  CRT  Graphics  II. 

20 

HP  74704  Graphics  Plottar. 

30 

HP  74754  Graphics  Plottar. 

40 

Epson  LQ- 1500  Prlntar,  singla  dansity. 

41 

Epson  LQ-1500  Prlntar,  donbla  dansity. 

42 

Epson  LQ-1500,  donbla  spaad,  dnal  dansity. 

43 

Epson  LQ-1500  Prlntar,  quad  dansity. 

51 

Houston  Inatruoant  ONP-51  HP  or 

DHP-52  HP  Plottar,  0.001"  stop  slxa. 

52 

Houston  Instrunant  DHP-51  HP  or 

DHP-52  HP  Plottar,  .005”  stap  size. 

60 

HP  26864  LasarJat  Prlntar  or  LasarJat 

PLUS  prlntar,  using  4  slxa  papar 
(8.5"  x  11")  (216  an  x  280  an). 

Drawing  rasolution:  75  dots  par  inch. 

61 

HP  26864  LasarJat  Prlntar,  using  B5  siza 
papar  (7.2N  x  10.1”)  (182  ■  x  257  an). 

Drawing  rasolut ion:  75  dots  par  inch. 

62 

HP  26864  LasarJat  Prlntar,  using  4  slxa 
papar  (8. 5"  x  11")  (216  an  x  280  an). 

Drawing  rasolution:  150  dots  par  inch. 

63 

HP  26864  LasarJat  Prlntar,  using  B5  slxa 
papar  (7.2"  x  10.1")  (182  an  x  257  an). 

Drawing  rasolution:  150  dots  par  inch. 

64 

HP  26864  LasarJat  Prlntar,  using  4  slxa 
papar  (8.5"  x  11")  (216  an  x  280  an). 

Drawing  rasolution:  300  dots  par  inch. 

65 

HP  26864  LasarJat  Prlntar,  using  B5  siza 
papar  (7.2"  x  10. l")  (182  an  x  257  an). 

Drawing  raaolutlon:  300  dots  par  inch. 

contlnuod 
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Table  2a.  (continued) 


•0 

IP  75808,  HP  7S8SB,  or  HP  75888  Drafting 

Plottar  oalag  also  A/M  to  D/Al  papar. 

IP  7S50A  Graphics  Plottar  using  slsa 

A/M  to  B/A3  papar. 

IP  7440A  Co lor Pro  plottar  using  sixa 

US/M  papar. 

85 

IP  75858  or  IP  75888  Drafting  Plottar 
using  alma  S/AO  papar. 

90 

Tektronix  4025. 

99 

IBM  color  graphics  non i tor  (CUT). 

Tab  la  2b.  Dot  Matrix  Printer  Usage  by  Modal 


Printer 

Nodal 

0 

D 

2 

D 

□ 

D 

D 

10 

u 

12 

13 

□ 

15 

□ 

Spoon  FX-80 

* 

a 

a 

a 

D 

D 

D 

■ 

■ 

B 

fl 

fl 

■ 

fl 

Spoon  HX-80 

* 

a 

a 

a 

I 

1 

■ 

fl 

1 

m 

■ 

m 

1 

H 

IBM  Printer 

• 

a 

a 

a 

1 

1 

9 

m 

1 

I 

1 

B 

■ 

Centronics  GUP 

a 

a 

a 

a 

■ 

I 

9 

fl 

1 

m 

m 

■ 

1 

m 

Okldata  92 

a 

a 

a 

a 

I 

■ 

1 

n 

fl 

j 

m 

m 

I 

Spoon  XX-80 

a 

a 

a 

a 

□ 

I 

Q 

fl 

1 

fl 

Spoon  FX-100 

■ 

■ 

9 

fl 

fl 

fl 

B 

B 

B 

B 

Spoon  MX- 100 

I 

B 

B 

I 

■ 

Okldata  93 

_ 

9 

■ 

9 

fl 

fl 

fl 

B 

■ 

II 

B 

*  ■  Tha  prlntar  eaa  aaa  this 


Table  3.  g«co— nd»d  Configuration 


Output  device 

Darien 

Modal 

Epson  FX-80 

0 

5 

Epson  MX- 80 

0 

1 

IBM  Pr intar 

0 

1 

Cantronica  QLP 

0 

1 

Okidata  92 

0 

1 

Epson  RX-80 

0 

1 

Epson  FX-100 

0 

15 

Epson  MX- 100 

0 

11 

LQ* 1300 

0 

41 

Okldata  93 

0 

11 

HI  DMP-51 

9800/9650 

51 

HI  DHP-S2 

9600/9650 

31 

HP  7 440 A 

9600/9650 

80 

HP  7470A 

9600/9650 

20 

HP  7475A 

9600/9650 

30 

HP  7S30A 

9600/9650 

80 

HP  73801 

9600/9650 

80 

HP  7383B 

9600/9650 

80/85 

HP  7386B 

9600/9650 

80/83 

HP  2686A 

9600/9650 

60/61 

Tektronix  4023 

4800/4850 

90 

IBM  color  graphics 
nonltor 

*99 

99 
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No  spaces  are  allowed  in  the  name.  For  further  information  see 
your  DOS  manual  on  sub-directory  creation  and  naming.  Press  ENTER 
or  DOWN  ARROW  to  advance  to  the  next  item. 

Type  the  directory  name  where  the  output  files  are  located. 
If  you  INSTALLed  the  SAMPLE  DATA  you  created  a  directory  called: 

C : \OUTPUT 

If  you  did  not  then  a  directory  has  not  been  created.  Select  and 
type  a  name  for  the  output  file  location  such  as: 

C:\OUTPUT 

Remember  this  name;  we  will  use  it  again.  Press  ENTER 

Upon  completion,  you  are  asked  if  you  wish  to  make  changes. 
If  everything  is  correct  type  N  for  no  or  else  Y  for  yes  and 
repeat  the  data  entry.  Use  DOWN  ARROW  keys  to  skip  over  accept¬ 
able  answers  and  change  only  what  is  in  error. 

If  you  did  not  use  the  INSTALL  for  SAMPLE  DATA  do  the  following. 
When  N  for  NO  CHANGES  is  pressed  the  program  returns  to  the 
OPTIONS  WINDOW.  Use  the  RIGHT  ARROW  key  to  EXIT.  When  out  of 
the  program  create  the  OUTPUT  LOCATION  directory  if  one  does  not 
exist.  Type  the  following: 

CD\ 

MD\OUTPUT 

CD\SEISMIC 

SEISMIC 

The  opening  screen  will  reappear  as  shown  above.  This  time  press 
ENTER  to  accept  the  DIRECTORY  choice.  A  screen  like  the  following 
will  appear  showing  the  earthquake  epicenter  files  on  disk. 


DIRECTOR?  F:\EPIC\».EPC 


11  file(s)  found 


E2.EPC 

E8.EPC 

El.EPC 

ES.EPC 

E3.EPC 

E10.EPC 

E4.EPC 

Ell.EPC 

ES.EPC 

Eh. ETC 

E7.EPC 

Filenane  :  E2.EPC 


Use  cursor  keys  to  select  file  then  press  <ENTER> 
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Press  ENTER  to  continue  and  a  directory  of  the  files  in  the 
output  location  directory  will  appear: 

•IRECTOR*  r:\SDAIAS*.* 


16  file(s)  found 


ONE .PL! 

TWO. PL 

TEST.3PT 

1V0.PLT 

TWO.OU 

RESPONSE.  PLT 

THREE. PL T 

TEST.30T 

ONE .OUT 

TEST.10T 

TUO.OUT 

IEST.1PT 

THREE. OUT 

TEST. EQS 

SQUARES 

ONE. PL 

Filenane  :  ONE. PL T 


Use  cursor  keys  to  select  file  then  press  <ENTER> 


Press  ENTER  to  return  to  the  OPTIONS  window  menu.  We  are  now 
ready  to  begin  a  problem. 

EARTHQUAKE  SELECTION 

The  epicenter  data  base  has  been  prepared  from  the  National 
Oceanic  and  Atmospheric  Administration's  data  base.  Each  file 
covers  a  specified  region  and  contains  date,  latitude,  longitude 
and  magnitude  data.  This  section  is  used  to  create  a  subset  of 
earthquakes  for  a  specific  region  from  the  main  set  of  earthquake 
files.  The  program  searches  a  rectangle  specified  by  maximum  and 
minimum  latitudes  and  longitudes  to  find  all  events  within  the 
box  above  the  specified  magnitude.  To  become  familiar  with  the 
program,  the  user  may  advance  to  the  EQ  SELECTION.  The  follow¬ 
ing  screen  appears: 
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ED  SELECTION 


Specify  Area 
Devise  Data 
Begin  Analysts 


View  Desalts 
Print  Results 
Plot  Results 
Save  Results 


M*  . 


*  •  . *ru 


USE  ARROW  KEYS  THEN  PRESS  ENTER 


Use  the  DOWN  ARROW  key  and  then  choose  REVISE  DATA  to  see  the 
example  problem.  Press  RETURN  to  accept  each  value  unchanged. 

Then  select  VIEW  RESULTS  from  the  EQ  SELECTION  window  to  see  the 
output  file  on  screen.  Selecting  PRINT  RESULTS  or  PLOT  RESULTS 
will  print  or  plot  the  sample  problem. 

Specify  Data 

This  choice  permits  data  entry  for  a  new  problem.  Selecting 
this  choice  will  overwrite  previous  data  examples  unless  they 
were  saved  to  a  named  file.  You  will  overwrite  the  sample  data 
but  it  can  be  copied  from  the  original  disk  again  if  needed.  The 
following  screen  will  appear: 


Enter  Maximum  Longitiude  (Degrees) 
example  128.0 

121 

Enter  Nininun  Longitude  (Degrees) 
example  112.0 

113 

Enter  Naximum  Latitude  (Degrees) 
example  40.0 

41 

Enter  Nininun  Latitude  (Degrees) 
example  34.0 

37 

Enter  Site  Longitude  (Degrees) 
example  115.0 

115 

Enter  Site  Latitude  (Degrees) 
example  38.0 

38 

Enter  Nininum  Magnitude  Cutoff 
example  3.0 

3.0 

Any  Changes?  V  /  N 
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Enter  the  longitude  and  latitude  for  a  rectangle  bounding  the 
problem  and  the  site  longitude  and  latitude.  The  study  area  must 
be  large  enough  to  include  distant  events  which  can  influence  the 
site.  Geologic  data  should  be  consulted  to  look  for  boundaries 
of  tectonic  provinces.  The  site  should  be  near  the  center  of  the 
study  region  unless  otherwise  required  by  a  tectonic  boundary. 

The  region  will  form  the  boundary  for  selecting  events  to  be  used 
in  the  study  and  are  relevant  to  establishing  the  site  seismic 
potential.  For  regional  studies  consideration  should  be  given  to 
selecting  an  area  large  enough  to  establish  the  regional  tectonic 
setting.  Enter  the  MINIMUM  MAGNITUDE  CUTOFF  value,  typically  3.0. 
Events  below  this  level  may  not  have  been  recorded  and  their 
absence  distorts  the  relationship  for  recurrence. 

The  program  computes  the  acceleration  at  the  site  location 
from  each  epicenter  location  where  a  magnitude  is  specified  in 
the  subset  of  earthquakes.  This  computation  uses  an  acceleration 
distance  attenuation  equation.  The  following  attenuation  rela¬ 
tionships  for  acceleration  are  included: 

1.  McGuire 

2 .  Geomatrix 

3.  Campbell 

4.  Idriss 

5 .  Donovan 

6.  Joyner  and  Boore 

The  user  is  has  the  option  of  selecting  which  to  use.  Since  each 
equation  analyzes  the  data  differently  and  since  there  is  a  large 
uncertainty  in  ground  motion  results  may  differ  depending  upon 
the  equation  selected.  The  user  is  encouraged  to  read  the  refer¬ 
ences.  The  following  screen  is  used  to  enter  the  acceleration 
equation  choice  and  the  hypocenter  depth  if  needed,  usually  10 
miles. 


1  He  Guire 

2  Geonatrix 

3  Canpbell 

4  Idriss 

5  Donovan  and  Bomstein 

6  Joyner  and  Boore 
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Enter  depth  to  hypocenter  (wiles) 
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ENTER  CHOICE  FOR  ACCELERATION  EQUATION 
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Once  the  acceleration  equation  data  entries  are  completed,  the 
user  is  asked  to  select  epicenter  files  to  be  searched  .  Files 
have  been  divided  into  regions  to  keep  the  search  time  to  a 
minimum.  A  screen  similar  to  the  following  is  shown  from  which 
the  user  may  move  through  the  list  using  the  ARROW  KEYS, 
PAGEUP/DOWN  KEYS  and  then  press  ENTER  to  select  his  choices. 
There  is  no  limit  to  the  number  of  choices. 


Calif 

Calif 

Calif 

Calif 

Calif 

Calif 

Calif 

Calif 

Calif 

Calif 

Calif 


114-119 

114-128 

114-121 

114-122 

114-123 

114-123 

114-123 

114-124 

114-125 

114-125 

114-125 


31- 32 

32- 33 

33- 34 

34- 35 

35- 36 

36- 37 

37- 38 

38- 39 

39- 40 

40- 41 

41- 42 


Press  ESCAPE  when  done.  The  program  then  shows  the  choices  and 
asks  for  confirmation. 


VOU  HAVE  SELECTED  THE  FOLLOWING  RECORDS: 


7  Calif  114-123  37-38 

8  Calif  114-124  38-39 

9  Calif  114-125  39-48 

10  Calif  114-125  40-41 

11  Calif  114-125  41-42 


IS  THIS  CORRECT?  V'N 
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Revise  Data 


If  the  user  wishes  to  revise  a  number  he  may  choose  the 
REVISE  DATA  choice  and  will  be  given  the  same  questions  as  in  the 
SPECIFY  DATA  section  with  the  previous  choices.  Press  enter  or 
DOWN  ARROW  to  accept  the  value  unchanged.  Overwrite  the  revised 
value  completely  to  alter  a  number. 

Begin  Analysis 

This  choice  begins  the  actual  data  search. 

View  Results 

The  VIEW  RESULTS  choice  permits  the  user  to  see  the  output 
file  from  an  analysis  on  the  screen.  It  must  be  run  after  an 
analysis  has  been  performed  and  data  exists. 

Print  Results 


PRINT  RESULTS  prints  the  output  files.  The  output  consists 
of  a  list  of  epicenters  with  the  computed  site  acceleration  for 
that  event,  and  a  histogram  of  the  distribution  of  acceleration. 

THE  FOLLOWING  IS  AN  EXAMPLE  OUTPUT 


MAX 

LONGITUDE 

73.000 

MIN 

LONGITUDE 

69.000 

MAX 

LATITUDE 

44.000 

MIN 

LATITUDE 

41.000 

MIN 

MAGNITUDE 

.000 

MAX 

MAGNITUDE 

9.000 

SITE  LONGITUDE 

72.000 

SITE  LATITUDE 

42.000 

ACCELERATION  EQUATION 

.000 

LIST  OF  SITE  EPICENTERS 


YEAR 

LATITUDE 

LONGITUDE 

MB 

MS 

MO 

ML 

AVM 

DISTANCE 

ACCEL 

1976. 

41.66 

69.97 

.00 

.00 

.00 

3.00 

3.00 

106.776 

.003 

1979. 

43.98 

69.80  3 

.80 

.00  4.00 

4.10 

3.97 

177.364 

.004 

1977. 

41.84 

70.70 

.00 

.00  3.10 

.00 

3.10 

67.614 

.006 

1974. 

41.70 

71.50 

.00 

.00 

.00 

2.50 

2.50 

32.976 

.008 

1976. 

41.56 

71.21 

.00 

.00 

.00 

3.50 

3.50 

50.660 

.011 

where 

MB 

MS 

MO 

ML 

AVM 

DISTANCE 

ACCEL 


Body  wave  magnitude 
Surface  wave  magnitude 
Other  magnitude 
Local  magnitude 
Average  magnitude 
Distance  event  to  site 

Site  acceleration  computed  estimate,  g's 
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Plot  File 


This  choice  creates  an  epicenter  plot  for  the  region. 

Figure  1  is  an  example  plot. 

Save  Results 

This  choice  permits  the  user  to  write  the  input,  output  and  plot 
files  to  a  named  file  of  the  user's  choice.  This  prevents  the 
files  from  being  overwritten. 
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73. 


.00 


44.00 


41.00 


4 
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Note:  numbers  are  approximate  magnitudes 
Figure  1.  Epicenter  plot  for  region. 
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REGIONAL  STUDY 


Moving  the  RIGHT  ARROW  key  to  REGIONAL  STUDY  reveals  the 
fol loving  screen: 


REGIONAL  STUDY 


Specify  Regie* 
Revise  Data 
Begin  Analysis 


View  Results 
Print  Results 
Plot  Results 
Save  Results 
Response  Plot 


USE  ARROW  KEYS  THEN  PRESS  ENTER 


This  section  is  used  to  conduct  a  regional  seismicity  study 
typically  for  eastern  sites  where  faulting  is  not  known  It  uses 
the  epicenter  subset  created  previously  .  This  section  will 
estimate  the  level  and  characteristics  of  the  earthquake  ground 
motion  which  pose  a  risk  to  the  site  of  interest.  He  will  use 
the  historical  epicenter  data  base  in  conjunction  with  geologic 
data  where  available  to  best  estimate  the  probability  of  site 
acceleration  levels.  This  becomes  the  basis  for  definition  of 
response  spectra  suitable  for  use  in  structural  design  and  analy¬ 
sis. 


For  a  regional  analysis  the  epicenter  data  base  is  used  to 
define  the  regional  magnitude  recurrence  based  on  a  log-linear 
fit  of  the  data  and  the  Richter  A  and  B  recurrence  coefficients 
are  determined. 


Log  (N)  -  A  +  B  M 

The  program  allows  the  user  to  specify  minimum  cutoff  magnitude 
to  enhance  the  fit.  Generally  the  epicenter  data  base  is 
deficient  on  small  events  since  events  less  than  magnitude  3  may 
not  be  large  enough  to  be  recorded  at  distant  seimograph 
stations.  Thus  a  cutoff  of  3  is  usually  used  to  insure  the  fit 
of  the  data  is  through  the  linear  portion  of  the  data.  Figure  2 
illustrates  the  case  in  which  the  initial  data  did  not  use  a 
cutoff  minimum  magnitude.  The  line  should  be  fit  through  the 
linear  portion.  The  program  may  be  repeated,  with  the  user 
specifying  the  A  and  B  coefficients  based  on  his  modification  to 
the  computed  fit  of  the  data. 
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100.000 


PROBABILITY  NOT  EXCEEDING  ft 


Having  established  the  magnitude  recurrence  relationship,  a 
"floating"  earthquake  analysis  is  performed  using  a  Monte  Carlo 
simulation  routine  A  series  of  events  representing  a  5000  year 
exposure  are  randomly,  spatially  assigned  with  magnitude  shaped 
by  the  recurrence  relationship.  A  maximum  cutoff  magnitude  may 
be  included  to  fit  the  tectonics  of  the  region.  A  probability 
distribution  for  a  specified  exposure  period  of  acceleration  at 
the  site  location  is  computed.  The  simulation  process  is  random 
in  location  and  includes  the  error  associated  with  the  earthquake 
attenuation  equations.  A  histogram  of  magnitude  and  acceleration 
is  computed. 

Use  the  DOWN  ARROW  key  and  then  choose  REVISE  DATA  to  see  the 
example  problem.  Press  RETURN  to  accept  each  value  unchanged. 

Then  select  VIEW  RESULTS  from  the  window  to  see  the  output  file 
on  screen.  Selecting  PRINT  RESULTS  or  PLOT  RESULTS  will  print  or 
plot  the  sample  problem. 

Specify  Region 

From  the  REGIONAL  STUDY  window  select  SPECIFY  REGION  and 
press  ENTER.  The  following  questions  will  appear: 


Enter  Maximum  Long  it  hide  (Degrees)  121 

example  120.0 

Enter  Minimum  Longitude  (Degrees)  116 

example  112.0 

Enter  Maximum  Latitude  (Degrees)  41 

example  40.0 

Enter  Nininun  Latitude  (Degrees)  36 

example  34.0 

Hinimum  Earthquake  Magnitude  Cutoff  3 

example  3.0 

Maximum  Earthquake  Magnitude  Cutoff  8.5 

example  8.0 

Enter  Site  Longitude  (Degrees)  117 

example  115.0 

Enter  Site  Latitude  (Degrees)  38 

example  38.0 


Any  Changes?  Y  /  N 
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Enter  the  longitude  and  latitude  of  a  rectangle  to  bound  the 
study  area.  This  may  be  a  smaller  than  the  region  chosen  for  the 
EARTHQUAKE  SELECTION  discussed  previously.  Consider  the 
tectonics  of  the  region  in  selecting  bounds  for  the  study.  Be 
sure  to  include  sufficient  distance  from  the  site  so  that  all 
events  which  can  cause  significant  ground  motion  at  the  site  are 
included.  Enter  the  MINIMUM  MAGNITUDE  CUTOFF,  usually  3  since 
events  less  than  3.0  may  not  have  been  recorded  and  distortion  of 
the  recurrence  estimation  might  result.  Enter  the  MAXIMUM 
MAGNITUDE  of  the  region. 

As  discussed  above,  enter  the  acceleration  equation  and 
hypocenter  depth  if  needed,  usually  10  miles.  Enter  the  EXPOSURE 
PERIOD  of  the  study  in  years. 


1  He  Guire 

2  Geonatrix 

3  Canpbell 

4  Idriss 

5  Donovan  and  Bomstein 

6  Joyner  and  Boore 

6 

ENTER  CHOICE  FOR  ACCELERATION  EQUATION 

Enter  depth  to  hypocenter  (niles)  B 

Enter  Exposure  Per  ion  (Years}  50 

exanple  58.8 
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At  this  point  you  will  be  asked  whether  to  use  the 
earthquake  data  base  or  to  use  regional  recurrence,  A  and  B 
values.  Enter  Y  to  use  the  epicenter  data.  If  N  is  entered 
enter  the  values  of  A  and  B  as  shown  here: 


Do  you  wish  to  confute  recurrence  data 
fro*  the  epicenter  data  base  (Y)es  or  (H)o 
If  you  answer  No,  you  mist  enter 
Richter  A  and  B  values  for  the  region 
Your  Choice?  ¥  /  N 


Enter  Richter  A  value 
exanple  4.2 

49 

Enter  Richter  B  value 
exanple  -.85 

-.3 

Revise  Data 

The  user  may  revise  data  once  entered  by  selecting  REVISE 
DATA  from  the  REGIONAL  STUDY  window.  The  same  questions  given  in 
SPECIFY  DATA  are  asked  with  the  previous  responses.  Advance 
through  the  data  using  the  DOWN  ARROW.  Overwrite  the  revised 
value  completely. 

Begin  Analysis 

This  choice  begins  the  actual  data  search. 

View  Results 

The  VIEW  RESULTS  choice  permits  the  user  to  see  the  output 
file  from  an  analysis  on  the  screen.  It  must  be  run  after  an 
analysis  has  been  performed  and  data  exists. 
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Print  Results 


PRINT  RESULTS  prints  the  output  files.  The  output  consists 
of  a  recurrence  data,  the  regional  recurrence  coefficients,  and 
the  site  probability  distribution. 


SAMPLE  OUTPUT 


M/DEG  LONGITUDE 

50.906 

M/DEG  LATITUDE 

69.057 

MAX  LONGITUDE 

73.000 

MIN  LONGITUDE 

69.000 

MAX  LATITUDE 

44.000 

MIN  LATITUDE 

41.000 

MIN  MAGNITUDE 

3.000 

SITE  LONGITUDE 

71.000 

SITE  LATITUDE 

42.500 

EXPOSURE  TIME 

50.000 

ACCELERATION  EQUATION  . 000 

MAGNITUDE  NUMBER  EVENTS 


.10 

10 

.20 

10 

.30 

10 

.40 

10 

8.70 

0 

8.80 

0 

8.90 

0 

9.00 

0 

GE  M  NUMBER  EVENTS/YR  GE  M 

.2564 
.2564 
.2564 
.2564 

.0000 

.0000 

.0000 

.0000 


YEARS  COVERED  39.00 

RICHTER  FIT  EQ  L0G10(N)=A  +  B*M 


A  =  3.698  B  =  -.909 


REGIONAL  PROBABILITY  DISTRIBUTION 


CONFIDENCE  ACCELERATION 

MAGNITUDE 

DISTANCE 

LOWER  CONF 

UPPER  COl 

.9900 

.2276 

5.5186 

13.8088 

.9705 

.9965 

.9802 

.2086 

4.8364 

14.0230 

.9537 

.9918 

.9704 

.1973 

5.2427 

16.7578 

.9390 

.9864 

.9608 

.1796 

5.6075 

20.4107 

.9254 

.9805 

.9512 

.1500 

4.9374 

13.8097 

.9125 

.9742 

.9418 

.1461 

5.0632 

13.8910 

.9002 

.9677 

.9324 

.1425 

5.4726 

13.8134 

.8883 

.9610 

.9231 

.1397 

6.0602 

14.2851 

.8768 

.9541 

.9139 

.1371 

5.0184 

13.8089 

.8656 

.9485 

.9048 

.1356 

5.2686 

13.9865 

.8547 

.9402 

.8958 

.1301 

4.2273 

14.7429 

.8440 

.9331 

.8869 

.1300 

4.7795 

16.7417 

.8335 

.9260 

.8781 

.1282 

4.6073 

14.9659 

.8233 

.9188 
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8694 

.1262 

5.2475 

15.1447 

.8133 

.9117 

8607 

.1239 

5.5734 

15.8349 

.8034 

.9044 

8521 

.1237 

5.3539 

18.6373 

.7939 

.8972 

sore 


HISTOGRAMS  COVERING  5000.0  YEARS 


NO.  P(X) 

F(X) 

HISTOGRAM  MAGNITUDE 

86  .086 

.086 

3.307 

********************************************* 

109  .109 

.195 

3.450 

********************************************* 

73  .073 

.268 

3.593 

********************************************* 

77  .077 

.345 

3.736 

************************************ '  ■ 

68  .068 

.413 

3.879 

************************************ 

51  .051 

.464 

4.022 

********************************************* 

66  .066 

.530 

4.165 

********************************************* 

33  .033 

.563 

4.308 

********************************* 

34  .034 

.597 

4.451 

********************************** 

36  .036 

.633 

4.594 

************************************ 

35  .035 

.668 

4.737 

*********************************** 

36  .036 

.704 

4.880 

************************************ 

33  .033 

.737 

5.023 

********************************* 

31  .031 

.768 

5.166 

******************************* 

20  .020 

.788 

5.309 

******************** 

more 

NO.  P(X) 

F(X) 

HISTOGRAM  ACCELERATION 

563  .563 

.563 

.007 

********************************************* 

214  .214 

.777 

.021 

************************************  Bore 

86  .086 

.863 

.036 

********************************* ***» 

47  .047 

.910 

.050 

********************************************* 

18  .018 

.928 

.064 

****************** 

13  .013 

.941 

.079 

************* 

9  .009 

.950 

.093 

********* 

9  .009 

.959 

.107 

********* 

6  .006 

.965 

.121 

****** 

6  .006 

.971 

.136 

****** 

4  .004 

.975 

.150 

**** 

2  .002 

.977 

.164 

** 

3  .003 

.980 

.179 

*** 

2  .002 

.982 

.193 

** 

2  .002 

.984 

.207 

** 

1  .001 

.985 

.221 

* 

0  .000 

.985 

.236 

2  .002 

.987 

.250 

** 

store 


Plot  File 

This  choice  creates  a  plot  of  recurrence  and  site 
acceleration  probability  for  the  region.  Figure  2  is  an  example 
plot. 
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Save  Results 

This  choice  permits  the  user  to  write  the  input,  output  and 
plot  files  to  a  named  file  of  the  user's  choice.  This  prevents 
the  files  from  being  overwritten. 

Response  Plot 

This  choice  creates  a  standard  shaped  site  independent 
response  plot  on  your  plotting  device.  Figure  3.  The  user  is 
requested  to  provide  the  base  ground  motion  level  in  g's. 
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VELOCITY  (IN/SEC) 


DAMPING  2.  5,  7.  AND  10  PERCENT 


PERIOD  (SEC) 


Figure  3.  Typical  response  spectra. 
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FAULTS  STUDY 


For  western  sites  where  fault  locations  are  known  in  store 
detail,  a  study  can  be  performed  by  specifying  the  location  of  a 
fault  in  terms  of  coordinates  of  several  points  defining  line 
segments.  All  earthquakes  within  a  specified  distance  or 
boundary  from  the  fault  line  under  study  are  made  a  subset,  and 
the  recurrence  of  the  fault  is  calculated.  A  probability 
analysis  is  then  performed  to  calculate  expected  site 
acceleration  and  causative  earthquake  magnitude  and  epicentral 
distance  for  that  fault.  The  program  randomly  selects  the 
epicenter  of  an  earthquake  somewhere  along  the  specified  length 
of  the  fault.  Using  the  fault  recurrence  data  in  terms  of 
Richter  coefficients  and  maximum  earthquake  magnitude  associated 
with  the  fault,  the  program  determines  and  earthquake  magnitude 
and  the  length  of  fault  break  (assumed  centered  on  the 
epicenter) ,  and  then  calculates  the  distance  of  the  site  to  the 
fault  break  (the  epicentral  distance  and  the  hypocentral 
distance) .  These  distances,  along  with  the  acceleration- 
magnitude  attenuation  relationship  with  its  uncertainty  defined 
by  the  standard  deviation  for  that  level  of  motion,  give  the  site 
acceleration.  The  process  is  repeated  using  a  Monte  Carlo  scheme 
to  produce  a  list  of  site  accelerations  and  related  causative 
magnitude  and  epicentral  distance  thus  defining  the  site's 
probability  distribution.  It  is  important  to  note  that  the 
program  is  random  for  each  and  every  fault  in  the  following: 

a)  Location  along  fault  length,  2  dimensional 

b)  Magnitude  shaped  by  recurrence  coefficient 

c)  Acceleration  level  using  mean  and  standard  deviation 
relationship  of  magnitude  -  distance 

Provisions  are  included  to  use  recurrence  data  from  slip 
analysis  or  regional  seismicity  in  lieu  of  fault  specific  data. 

The  program  determines,  tabulates,  and  plots  the  probability  of 
not  exceeding  various  levels  of  acceleration  in  the  time  period 
specified.  These  data  are  available  for  all  individual  faults,  and 
then  are  combined  for  all  faults  acting  together.  The 
determination  of  total  risk  to  the  site  is  of  importance  for 
establishing  design  levels. 

A  single  recent  large  event  may  release  strain  built  up  over 
hundreds  or  thousands  of  years.  As  such,  it  might  indicate  a 
period  of  less  activity  in  the  immediate  future.  However,  since 
the  data  base  is  relatively  short,  the  return  time  for  this  event 
might  be  erroneously  indicated  as  much  less.  For  example,  if 
this  were  a  500-year  event  and  occurred  during  a  50-year  data 
base  its  return  might  be  estimated  at  0.02  rather  than  0.002. 

The  plotted  data  points  and  line  of  best  fit  are  determined  by 
the  computer  analysis  using  regression  analysis  techniques. 

These  should  be  reviewed  and  judgment  used  to  adjust  this  type  of 
datum  point  that  will  clearly  plot  significantly  higher  than  the 
linear  portion  of  the  recurrence  data. 

The  procedure  is  intended  to  be  repeated  several  times. 
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during  which  the  engineer  can  compare  geologic  data  and  lines  of 
best  fit  from  historic  data  and  converge  on  the  best  estimate 
using  his  judgment. 

The  FAULTS  STUDY  window  reveals  the  following  choices: 


FAULTS  STUDY 


Devise  Date 
Beg  in  Analysis 


View  Results 
Print  Results 
Plot  Results 
Plot  Faults 
Save  Results 
Response  Plot 


Edit  Faults.Lst 


USE  ARROV  KEYS  THEM  PRESS  ENTER 


Use  the  DOWN  ARROW  key  and  then  choose  REVISE  DATA  to  see  the 
example  problem.  Press  RETURN  to  accept  each  value  unchanged. 

Then  select  VIEW  RESULTS  from  the  window  to  see  the  output  file 
on  screen.  Selecting  PRINT  RESULTS  or  PLOT  RESULTS  will  print  or 
plot  the  sample  problem. 

Specify  Data 

This  is  the  data  entry  section  for  a  new  problem.  Selecting 
this  choice  shows  the  following  screen: 
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Enter  Ibxinun  Longitiude  (Degrees)  J2B 

example  120.8 

Enter  Nininun  Longitude  (Deyr sen)  1M 

example  112.8 

Enter  Naximum  Latitude  (Degrees)  42 

nxanple  40.8 

Enter  Minimum  Latitude  (Degrees)  36 

example  34.8 

flininun  Magnitude  Cutoff  3 

exanple  3.0 

Enter  Site  Longitude  (Degrees)  117 

exanple  115.8 

Enter  Site  Latitude  (Degrees)  39 

exanple  38.8 

Enter  Exposure  Period  (Years)  50 

exanple  50. 8 


Any  Changes?  Y  /  N 

Enter  the  longitude  and  latitude  of  a  rectangle  to  bound  the 
study  area.  This  may  be  a  smaller  than  the  region  chosen  for  the 
EARTHQUAKE  SELECTION  discussed  previously.  Consider  the 
tectonics  of  the  region  in  selecting  bounds  for  the  study.  Be 
sure  to  include  sufficient  distance  from  the  site  so  that  all 
events  which  can  cause  significant  ground  motion  at  the  site  are 
included.  Enter  the  MINIMUM  MAGNITUDE  CUTOFF,  usually  3  since 
events  less  than  3.0  may  not  have  been  recorded  and  distortion  of 
the  recurrence  estimation  might  result.  Enter  the  EXPOSURE 
PERIOD  of  the  study  in  years. 
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Enter  the  acceleration  equation  to  use  as  discussed  above. 


1  lie  Guire 

2  Geonatrix 

3  Canpbell 

4  Idriss 

5  Donovan  and  lomstein 

6  Joyner  and  Boore 

6 

ENTER  CHOICE  FOR  ACCELERATION  EQUATION 


The  user  is  then  given  a  menu  of  known  faults  in  the  California 
area.  This  menu  is  NOT  complete;  but  is  meant  as  a  vehicle  to 
reduce  data  entry.  The  user  may  add  to  this  list;  this  will  be 
discussed  later. 


SAN  CLENENTE 
PALOS  VERSE 
NEWPORT- INGLEWOOD 
ROSE  CANYON 
UNRATED 
ELSINORE 
WHITTIER 
SAN  JACINTO 
COYOTE  CREEK 
SUPERSTITION  NT 
SUPERSTITION  HILL 
INFER IAL 
BANNING 
-  SAN  ANDREAS 
XY  LOCAL 
CRISTIANITOS 
ALISO 


Use  the  ARROW  KEYS  or  press  the  first  letter  of  the  fault  name 
desired  to  move  throught  the  list.  Press  ENTER  to  select  a  fault; 
up  to  15  faults  may  be  selected. 
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The  question  is  then  asked  whether  to  enter  additional  faults. 
If  the  answer  is  yes  the  following  screen  is  shown: 


Enter  Fault  Nane 
exanple  San  Andreas 


0  i>.l:  ji*fc  - : i-i : 


Enter  Haxinun  Magnitude  iHTl 

exanple  7.8 

Enter  Characteristic  Magnitude  IE 

exanple  7.5 

Enter  Return  Tine  in  Tears  JO 

exanple  258.8 

Enter  Designator  1  or  2  or  3  13 

1  for  2  line  segnent 

2  for  box  segnent 

3  for  specification  of  Richter  A  and  B 


Enter  Depth  to  Hgpocenter  (niles)  HE! 

exanple  5.8 


Ang  Changes?  V  /  N 


Enter  the  name  of  the  fault  and  the  faults  MAXIMUM  MAGNITUDE. 

The  CHARACTERISTI  MAGNITUDE  is  the  magnitude  which  a  geological 
evidence  shows  ha  a  frequent  history  of  occuring  as  a  major 
event  and  the  RET  RN  TIME  requested  is  the  return  time  of  that 
event.  This  evenc  is  added  to  the  seismicity  computed  by  the 
historical  data,  '"'he  CHARACTERISTIC  MAGNITUDE  may  exceed  the 
MAXIMUM  MAGNITUDE  which  is  the  cutoff  for  the  recurrence 
calculation  based  on  the  historical  data  or  the  input 
coefficients. 
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The  fault  may  be  described  by  three  choices: 

1.  A  2  line  segment  where  events  a  specified  distance  from 
the  fault's  line  segments  are  included  in  the  subset  of 
earthquake  events  used  to  calculate  the  recurrence  of  the  fault. 

2.  A  2  line  segment  where  earthquake  events  within  a  4  sided 
region  are  used  to  calculate  the  recurrence  of  the  fault. 

3.  A  2  line  segment  where  the  recurrence  of  the  fault  is 
specified  in  terms  of  the  Richter  A  and  B  values. 

Enter  your  choice  for  the  specific  fault  being  defined. 

If  you  select  1: 

The  following  questions  will  be  asked 


Enter  Feint  1  Longitude  (Degrees) 
exanple  129.0 

128 

Enter  Point  1  Latitude  '’Jegrees) 
exanple  49.9 

49 

Enter  Point  2  Longitude  (Degrees) 
exanple  118.9 

tie 

Enter  Point  2  Latitude  (Degrees) 
exanple  38.9 

37.5 

Enter  Point  3  Longitude  (Degrees) 
exanple  117.9 

115.3 

Enter  Point  3  Latitude  (Degrees) 
exanple  37.9 

37 

Distance  fron  fault  to  include  Events  (niles) 
exanple  19.9 

19 

Any  Changes?  ¥  *  8 


Enter  Fault  Dip,  degrees 
exanple  45.8 

Enter  Type  of  Fault  Designator  1  or  2  or  3 
t  Tor  Strike  Slip 

2  for  Oblique 

3  for  Thrust 
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These  define  the  line  segments  and  distance  away  from  the  fault, 
see  Figure  4. 

If  you  select  2: 


The  following  questions  will  be  asked  as  above. 


Enter  Point  1  Lone  it  indie  (Degree*) 
example  128.8 

12B 

Enter  Point  1  Latitude  (Degrees) 

example  48.0 

48 

Enter  Point  2  Longitude  (Degrees) 
example  118.8 

118 

Enter  Point  2  Latitude  (Degrees) 
example  38.0 

37.5 

Enter  Point  3  Longitude  (Degrees) 
example  117.0 

116.3 

Enter  Point  3  Latitude  (Degrees) 
example  37.0 

37 

Additionally  coordinates  for  a  4  sided  region  must  be  entered. 
Start  with  the  uppermost  right  point  and  proceed  CLOCKWISE.  See 
Figure  5. 

CUTER  4  POINTS  FOR  BOX  STARTING  AT  UPPER  RIGHT  AND  GO  CLOCKWISE 

Enter  Point  1  Long  it iude  (Degrees)  188 

example  128.0 

Enter  Point  1  Latitude  (Degrees) 
example  40.0 

40 

Enter  Point  2  Longitude  (Degrees) 
example  118.0 

118.5 

Enter  Point  2  Latitude  (Degrees) 
example  39.5 

38 

Enter  Point  3  Longitude  (Degrees) 
example  118.5 

128.5 

Enter  Point  3  Latitude  (Degrees) 
example  37.8 

37.5 

Enter  Point  4  Longitude  (Degrees) 
example  128.5 

128 

Enter  Point  4  Latitude  (Degrees) 
example  37.5 

48 
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If  you  select  3: 


Enter  the  points  for  the  2  line  segments  as  above. 


Enter  Point  1  Longitiude  (Degrees) 

120 

exanple  128.0 

Enter  Point  1  Latitude  (Degrees) 

40 

example  40.0 

Enter  Point  2  Longitude  (Degrees) 

118 

exanple  118.0 

Enter  Point  2  Latitude  (Degrees) 

37 .5 

exanple  38.0 

Enter  Point  3  Longitude  (Degrees) 

116.3 

exanple  117.0 

Enter  Point  3  Latitude  (Degrees) 

37 

exanple  37.0 

Any  Changes?  H  /  N 


Then  enter  the  A  and  B  values. 


Enter  Richter  A  for  Fault 
exanple  4.0 

4.9 

Enter  Richter  B  for  Fault 
exanple  -.89 

-.9 
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Revise  Data 


The  user  may  revise  data  once  entered  by  selecting  REVISE 
DATA  from  the  FAULTS  STUDY  window.  The  same  questions  given  in 
SPECIFY  DATA  are  asked  with  the  previous  responses.  The  user  may 
revise  the  data  specified  in  the  predefined  faults  list. 


DO  ¥0U  WISH  TO  REVISE  DATA  FOR  THE  FAULTS 
YOU  SELECTED  FROfl  THE  PREDEFINED  LIST 
(VIES  /  (M)0 


Selection  of  N  for  NO  allows  the  user  to  accept  or  revise 
the  selection  of  faults  from  the  predefined  list  of  faults  in  the 
menu.  The  data  for  these  faults  can  not  be  revised  for  this 
choice.  However,  the  user  may  revise  the  data  for  the  faults  he 
entered.  Selection  of  Y  for  YES  treats  the  predefined  faults  as 
if  they  were  user  entered  and  allows  the  user  complete  freedom  to 
change  all  values.  In  editing  faults  the  user  may  skip  a  fault 
by  pressing  the  ESCAPE  KEY  after  the  fault  name  appears  rather 
than  the  ENTER  or  DOWN  ARROW  KEY.  This  advances  to  the  next 
fault  leaving  the  previous  fault's  values  unchanged. 

NOTE  AND  WARNING 

THE  PREDEFINED  FAULTS  ARE  REASONABLE  FIRST  ESTIMATES.  THEY 
WILL  NOT  SUIT  ALL  CASE  STUDIES.  IN  PARTICULAR,  THE  USER  SHOULD 
PAY  ATTENTION  TO  THE  DEFINITION  OF  THE  REGION  SURROUNDING  THE 
FAULT  LINE  TO  INCLUDE  EARTHQUAKE  EPICENTERS  FOR  THE  FAULT 
RECURRENCE  RELATIONSHIP.  THE  INTENT  IS  TO  INCLUDE  ONLY  THOSE 
EVENTS  WHICH  CAN  BE  ATTRIBUTED  TO  THAT  FAULT  AND  EXCLUDE  THOSE 
FROM  OTHER  FAULTS.  BASED  ON  YOUR  SPECIFIC  PROBLEM  IT  MAY  BE 
NECESSARY  TO  CHANGE  THE  SELECTION  MODE  FROM  A  STANDARD  DISTANCE 
AWAY  FROM  THE  FAULT  TO  A  QUADRILATERAL  DEFINITION.  AFTER  A 
PRELIMINARY  ANALYSIS  FIRST  RUN  SOLUTION,  THE  USER  SHOULD  EXAMINE 
EACH  RECURRENCE  CURVE  AND  INCORPORATE  GEOLOGIC  DATA.  THE  SLOPE  OF 
THE  RECURRENCE  LINE  SHOULD  BE  CHECKED  TO  INSURE  IT  PASSES  THROUGH 
THE  LINEAR  PART  OF  THE  EVENTS  DATA.  THE  USER  SHOULD  REVISE  THE 
FAULT  DATA  FOR  THE  ALPHA  AND  BETA  COEFFICIENTS  COMPUTED  FROM  THE 
INITIAL  PRELIMINARY  ANALYSIS.  THIS  ITERATION  CONTROLS  THE 
QUALITY  OF  THE  ANALYSIS.  THE  ACCURACY  OF  THE  ANALYSIS  IS  A 
FUNCTION  OF  THE  EFFORT  SPENT  BY  THE  USER  TO  DEVELOP  THE  MODEL. 

THE  PROGRAM  HAS  THE  CAPABILITY  OF  PRODUCING  HIGHLY  ACCURATE 
RESULTS . 
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Begin  Analysis 

This  choice  begins  the  actual  data  search. 

View  Results 

The  VIEW  RESULTS  choice  permits  the  user  to  see  the  output 
file  from  an  analysis  on  the  screen.  It  must  be  run  after  an 
analysis  has  been  performed  and  data  exists. 


Print  Results 

PRINT  RESULTS  prints  the  output  files.  The  output  consists 
of  fault  recurrence  data,  the  fault  recurrence  coefficients, 
and  the  fault  and  total  site  probability  distribution. 

SAMPLE  OUTPUT 


M/DEG  LONGITUDE 

57.710 

M/DEG  LATITUDE 

69.057 

MAX  LONGITUDE 

119.000 

MIN  LONGITUDE 

115.000 

MAX  LATITUDE 

34.500 

MIN  LATITUDE 

32.000 

MIN  MAGNITUDE 

3.000 

SITE  LONGITUDE 

117.360 

SITE  LATITUDE 

33.300 

EXPOSURE  TIME 

50.000 

ACCELERATION  EQUATION 
SAN  CLEMENTE 

.000 

INDIVIDUAL  FAULT  STUDY 

FAULT  COORDINATES 

118.700  32.200 

118.300  32.800 

117.800  32.650 


FAULT  MAX  CREDIBLE  EARTHQUAKE  7.70 


FAULT  EPICENTERS 

LONGITUDE 

LATITUDE 

MAGNITUDE 

DIST 

ACC 

117.800 

32.583 

4.500 

55.645 

.018 

117.833 

32.800 

4.000 

44.015 

.017 

117.833 

32.716 

4.400 

48.699 

.020 

117.866 

32.800 

3.100 

45.221 

.009 

more 

AVE 

ACCELERATION 

0110 

MAX 

ACCELERATION 

0378 
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EARTHQUAKE  RECURRENCE 


MAGNITUDE 

NUMBER  EVENTS  GE  M 

NUMBER  EVENTS/ YR  GE  M 

.10 

45 

.6164 

.20 

45 

.6164 

.30 

45 

.6164 

.40 

45 

.6164 

8.70 

0 

.0000 

8.80 

0 

.0000 

8.90 

0 

.0000 

9.00 

0 

.0000 

YEARS 

COVERED 

73.00 

RICHTER  EQUATION  L0G10(N)  =  A  +B 

M 

• 

A  =  1.640 

B=  -.601 

FAULT  REGION  ACCELERATION  PROBABILITY 

CONFIDENCE  ACCELERATION  MAGNITUDE 

DISTANCE  LOWE  R  CONF 

UPPER  CONF 

.9900 

.3014 

7.2624 

30.0969 

.9705 

.9965 

.9802 

.2364 

6.7600 

30.0901 

.9537 

.9918 

.9704 

.2084 

6.0210 

34.3689 

.9390 

.9864 

.9608 

.1984 

7.3362 

30.2334 

.9254 

.9805 

.9512 

.1799 

6.7060 

30.1197 

.9125 

.9742 

.9418 

.1795 

5.9095 

30.5953 

.9002 

.9677 

.9324 

.1606 

6.6048 

53.9805 

.8883 

.9610 

.9231 

.1571 

7.4193 

55.4437 

.8768 

.9541 

.9139 

.1548 

6.9925 

31.2901 

.8656 

.9485 

.9048 

.1306 

6.7041 

34.9663 

.8547 

.9402 

.8958 

.1292 

7.5633 

30.1013 

.8440 

.9331 

.8869 

.1287 

5.6727 
_ 5"^ 

30.0893 

.8335 

.9260 

A 

SIMILAR  DATA  TABULATION 

IS  GIVEN 

FOR 

THE  TOTAL 

PROBABILITY  ACCELERATION 

FROM  ALL 

OF  THE  FAULTS 

Plot  File 

This  choice  creates  a  plot  of  recurrence  and 
probability  data  for  each  fault  and  the  total  site  probability  of 
acceleration.  Figure  6  is  an  example  plot. 
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MAGNITUDE 


Note:  The  user  should  review  the  fit  of  the  recurrence  data 
shown  by  the  line.  This  line  should  pass  through  the  linear 
portion  of  the  data.  Sometimes  when  the  number  of  events  is 
limited  or  the  minimum  event  is  not  specified  the  line  may 
miss  that  region.  The  user  should  then  rerun  the  problem 
revising  the  fault  data  by  entering  the  A  and  B  values  from 
a  new  fit  correctly  through  the  data. 


Figure  6a.  Recurrence  for  a  fault. 


PROBABILITY  NOT  EXCEEDING  A 


Plot  Faults 


This  choice  creates  a  plot  of  the  faults  entered  in  the 
SPECIFY  DATA  or  REVISE  DATA  option.  The  epicenters  in  the 
EQUAKES  file  are  plotted  also.  This  choice  requires  that  the 
data  be  specified  for  the  faults  and  the  epicenter  file  be 
created.  However  the  analysis  need  not  be  performed,  so  this  may 
be  used  to  check  the  events  near  faults. 

The  plotting  procedure  draws  the  grid  of  the  region 
specified  and  then  draws  the  faults  and  plots  the  epicenters. 
Often  the  fault  extends  beyond  the  region  specified.  IT  IS  NORMAL 
IN  SUCH  CASES  TO  SEE  AN  ERROR  MESSAGE  STATING  UNPLOTTED  VECTORS 
OR  CLIPPED  VECTORS  for  the  elements  which  go  beyond  the  plot 
boundaries . 


Save  Results 

This  choice  permits  the  user  to  write  the  input,  output  and 
plot  files  to  a  named  file  of  the  user's  choice.  This  prevents 
the  files  from  being  overwritten. 

Response  Plot 

This  choice  creates  a  standard  shaped  site  independent 
response  plot  on  your  plotting  device.  Figure  3.  The  user  is 
requested  to  provide  the  base  ground  motion  level  in  g's. 

Edit  Faults . Lst 

This  option  permits  the  user  to  change  the  values  of  the 
data  used  to  define  the  faults  in  the  predefined  faults  list 
which  appears  as  a  menu  from  which  the  user  may  select  faults  for 
an  analysis.  Appendix  A  explains  how  to  use  this  option  and  also 
explains  how  to  calculate  the  recurrence  coefficients,  A  and  B 
from  a  recurrence  plot.  The  appendix  gives  a  list  of  faults  and 
the  recurrence  plots  for  all  the  faults  defined  in  the  FAULTS. LST 
data  file  of  predefined  faults.  The  user  may  add  to  this  list. 

NOTE  recurrence  plots  for  faults  with  predefined  recurrence 
coefficients  are  not  plotted  to  save  time;  refer  to  Appendix  A 
for  the  plots. 
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PROGRAM  LIMITATIONS 

The  EARTHQUAKE  SELECTION  routine  can  process  20,000 
epicenters  to  the  file  EQUAKES. 

The  REGIONAL  STUDY  can  input  20,000  epicenters. 

The  FAULTS  STUDY  can  input  20,000  epicenters  and  30  faults. 
Each  fault  can  have  up  to  5,000  epicenters  associated  with  it. 

To  keep  the  computational  process  manageable  within  the 
limits  of  a  desktop  computer  tradeoffs  had  to  be  made  in  Monte 
Carlo  process.  As  can  be  seen  the  95  percent  confidence  bounds 
increase  as  the  probability  of  not  exceeding  the  acceleration 
increase.  This  program  is  intended  to  give  an  estimation  of 
earthquake  motions  up  to  the  90  percent  probability  of  not  being 
exceeded  in  50  years  with  high  confidence.  This  program  is  not 
intended  to  be  used  to  predict  5,000  year  or  10,000  year  events. 

The  user  constructs  a  model  of  the  seismicity  of  a  region  to 
which  epiceters  are  assigned,  recurrence  relationships  computed 
and  ground  motion  probabilities  determined.  It  should  be  obvious 
that  if  the  user  fails  to  include  sufficient  fault  definition  the 
model  will  lack  those  details  and  their  contribution  to  the 
site's  expected  motion.  Regional  studies  can  not  be  substituted 
for  fault  studies  with  the  same  level  of  accuracy.  While  this  is 
a  simple  but  accurate  analytical  tool  in  the  hands  of  a  trained 
engineer  it  can  be  misused  by  the  uninformed. 

WARNING 

CHECK  YOUR  RESULTS  TO  INSURE  LESS  THAN  20,000  EVENTS  WERE  SELECT¬ 
ED. 

The  program  will  display  the  number  of  events  selected  and  the 
number  of  events  read.  If  5000  events  were  read  reduce  your 
area  if  possible  or  break  the  analysis  into  2  parts.  Failure  to 
do  this  will  result  in  the  omission  of  all  events  beyond  the  5000 
event  limit  thus  reducing  the  coverage. 

FOR  THE  ADVANCED  USER 

Data  and  Output  Files 

The  following  shows  the  input  data  files  created  and  the 
output  results  files  for  each  section  of  the  program. 


CONFIGURATION 
EARTHQUAKE  SELECTION 


INPUT 


ONE. IN 
EPICENTR. LST 
epicenters  .EPC 


OUTPUT 

DEVICE. CFG 

ONE . OUT 
ONE . PLT 
EQUAKES 
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REGIONAL  STUDY 


TWO. IN  TWO. OUT 

EQUAKES  TWO.PLT 


FAULTS  STUDY 


THREE. IN 
FAULTS . LST 


THREE. OUT 
THREE. PLT 


Files  with  the  extension  . EPC  are  epicenter  files.  Files 
with  the  extension  .IN  are  data  files  created  by  the  program. 
Files  with  the  extension  .LST  are  files  supplied  to  the  user 
which  the  user  can  modify;  these  will  be  discussed  later.  Files 
with  the  extension  .OUT  are  output  files  to  be  printed.  Files 
with  the  extension  .PLT  are  plot  files. 


When  files  are  saved  the  above  files  are  copied  to  disk  with 
a  new  name  specified  by  the  user  as  follows: 


EARTHQUAKE  SELECTION 


REGIONAL  STUDY 


FAULTS  STUDY 


ONE. IN 

XXXXXXXX. 1IN 

ONE. OUT 

XXXXXXXX. 10U 

ONE. PLT 

XXXXXXXX . 1PT 

TWO. IN 

XXXXXXXX. 2 IN 

TWO. OUT 

XXXXXXXX. 20U 

TWO.PLT 

XXXXXXXX. 2 PT 

THREE. IN 

XXXXXXXX. 3IN 

THREE. OUT 

XXXXXXXX. 30U 

THREE. PLT 

XXXXXXXX. 3PT 

The  program  is  set  to  use  the  standard  names  in  the  left  column. 
To  restore  a  saved  data,  output  or  plot  file  for  use  copy  the 
file  to  the  standard  name  as  follows: 


COPY  XXXXXXXX . 1IN  C : \SEISMIC\ONE . IN 
COPY  XXXXXXXX. 1PT  ONE. PLT 

Note  that  the  INPUT  files  are  stored  in  the  PROGRAM  DIRECTORY 
and  the  OUTPUT  and  PLOT  files  are  stored  in  the  DATA/RESULTS 
DIRECTORY.  The  program  will  re-run  or  revise  the  last  case 
executed  with  the  standard  names.  To  revise  or  re-run  a  case, 
only  the  INPUT  file  (.IN)  need  be  copied  to  the  standard  name. 

Epicenter  Files 

The  epicenter  data  base  is  broken  down  into  separate  files 
covering  small  regions  to  minimize  search  time.  The  user  may  add 
to  the  list  of  epicenter  file  by  editing  the  file  EPICENTR. LST 
which  contains  the  number  of  files  and  for  each  epicenter  file 
in  the  data  base  the  file  description  and  file  name. 


Number  of  Epicenter  Files 
Description  ,  xxxxxxx.EPC 
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The  epicenter  data  files  contain  data  in  the  following 
format : 


Year 

Col 

5 

to 

Col 

8 

F4.0 

Latitude 

Col 

20 

to 

Col 

24 

F5.3 

Longitude 

Col 

26 

to 

Col 

31 

*6.3 

BM  Body  Magnitude 

Col 

36 

to 

Col 

38 

FI.  2 

SM  Surface  Magnitude 

Col 

54 

to 

Col 

55 

FI.  1 

OM  0'  her  Magnitude 

Col 

61 

to 

Col 

63 

FI.  2 

IK  Local  Magnitude 

Col 

83 

to 

Col 

85 

FI.  2 

Note  no  decimal  points  used. 
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User’s  Guide 

Optimized  Site  Matched  Spectra 
Appendix  B 


INTRODUCTION 


Earthquake  resistant  design  and  analysis  utilize  response  spectra 
and  acceleration  time  histories  as  quantification  of  the 
earthquake  loading.  This  program  will  permit  the  user  to  select 
a  set  of  earthquake  records  from  a  data  base  of  records  collected 
by  the  California  Institute  of  Technology  and  currently  being 
distributed  by  the  National  Oceanographic  and  Atmospheric 
Administration.  The  program  allows  straight  selection  of  one  or 
more  records  or  allows  collection  of  a  set  of  up  to  10  records 
matched  to  site  classification,  acceleration  level,  site  to 
source  distance,  and  earthquake  magnitude  each  specified  by  the 
user. 

SYSTEM  REQUIREMENTS  AND  INSTALLATION 

The  program  is  intended  to  run  on  a  personal  computer  with  the 
following: 

MS  DOS  3.2  or  later 
64 OK  memory 
hard  disk 

80287  math  coprocessor  chip 

The  program  will  run  on  a  Zenith  248  computer  with  the  above.  The 
program's  earthquake  data  base  comprises  about  40  megabytes  of 
data.  It  is  not  practical  to  distribute  this  on  disk.  Currently 
the  data  base  resides  on  5.25  -  disk  cartridge  in  compressed 
form. 

To  obtain  paper  plots  of  the  data  one  or  more  of  the  following 
can  be  used: 


Epson  FX80  printer 
Hewlett  Packard  Laser jet  printer 
Hewlett  Packard  Graphics  Plotters 
Houston  Instruments  Plotters 
Tecktronix  4025 

The  laser  printer  is  a  rapid  means  for  producing  high  quality 
plots  and  is  recommended.  Specific  models  will  be  given  in  the 
configuration  section.  The  programs  are  distributed  on  floppy 
disk  but  are  intended  to  be  run  from  a  hard  disk.  Installation 
consists  of  placing  the  program  disk  in  Drive  A,  and  typing 

INSTALL 

This  will  run  a  batch  file  which  will  create  a  directory  on  the  C 
drive  (hard  disk)  named  "EQ”  and  will  copy  all  the  programs  into 
that  directory.  If  you  wish  to  create  a  directory  in  another 
location  or  with  another  name  you  may  do  so  by  using  the  DOS 
command  MD  for  make  directory  and  then  copy  the  program  disk 
files  to  that  directory. 
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The  following  are  the  files  which  are  required  to  run  the 
program: 


EQUAKE.EXE 

Main  MODULE 

E QUAKE 1 . EXE 

Overlay  1 

EQUAKE2.EXE 

Overlay  2 

EQUAKE 3 . EXE 

Overlay  3 

EQUAKE4 . EXE 

overlay  4 

EQ1.SCN 

Part  of  EQUAKE.EXE 

EQ2.SCN 

Part  of  EQUAKE.EXE 

EQ3.SCN 

Part  of  EQUAKE.EXE 

EQ4.SCN 

Part  of  EQUAKE.EXE 

EQ5.SCN 

Part  of  EQUAKE.EXE 

EQ6.SCN 

Part  of  EQUAKE.EXE 

RECORD. IDX 

Earthquake  record  index 

RECORD. DTA 

earthquake  record  identifiers 

Data  filet  on  the  Bernoulli  Box  cartridge  are  labeled  as 

xxxx.THS  For  time  history  data 

xxxx . RES  For  response  spectra  data 

one  additional  file  is  created  during  the  configuration  for  the 
individual  computer  setup. 

DEVICE. CFG 

This  file  comes  with  the  default  parameters  for  data  storage 
locations  and  a  laser  printer  connected  to  port  LPT1. 


CONFIGURATION 

To  begin  the  program  move  to  the  EQ  directory  and  type 

EQUAKE 

This  will  show  the  following  opening  screen: 


Use  the  RIGHT  ARROW  key  to  advance  to  OPTIONS  and  press  RETURN; 
the  following  screen  will  appear. 


=  Record  Natch  Time  History  Response  Spectra 

Options 

Exit 

Configure 

USE  ARROW  REVS  THEM  PRESS  ENTER 


Press  RETURN  again  and  the  following  questions  will  be  seen: 


Enter  Plotter  DEVICE  NUMBER  1 

see  User's  Manual  e.g.  LPI1  =  1 

Enter  Plotter  MODEL  NUMBER  60 

see  User’s  Manual  e.g.  HP  Lasejst  =60 

Enter  Location  (or  Tine  History  Records  1  to  500  E*. 

example  E: 

Enter  Location  for  Tine  History  Records  501  to  end  E: 

example  E: 

Enter  Location  for  Response  Spectra  Records  1  to  500  E:VRES1-500V 
example  E:\RES1-500S 

Enter  Location  for  Response  Spectra  Records  501  to  end  E:\RES500-S 
example  E’.SRES500-S 

Any  Changes?  V  /  N 
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The  DEVICE  NUMBER  refers  to  the  port  to  which  the  hard  copy 
plotting  device  is  connected;  see  Table  1  for  the  configuration 
options  and  check  your  manual  which  came  with  the  hard  copy  plot 
device.  The  MODEL  NUMBER  can  be  obtained  from  Table  2a.  Table  2b 
gives  a  matrix  of  MODEL  NUMBERS  for  various  compatible  printers. 
Table  3  gives  the  recommended  configuration  for  specific  devices. 

Enter  the  location  for  the  data  file  supplied  on  cartridge  by 
giving  the  DRIVE  designation  and  DIRECTORY.  To  verify  the 
location  you  may  leave  the  program  and  switch  to  the  Bernoulli 
Box  and  perform  a  directory  check  by  typing  •’DIR" .  To  leave  the 
program  press  RETURN  to  accept  the  default  values  then  press 
ESCAPE  several  times. 

Once  the  information  is  entered  correctly  press  "N"  for  NO 
CHANGES  and  you  may  return  to  the  opening  screen. 
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Table  1.  Device  Nuaber 


Device 

Output  Device 

Nuaber 

Parallel  Port 

o 

PUN:  (PRN 

:  la  equivalent  to  LPT1: 

) 

1 

LPT1: 

2 

LPT2: 

3 

LPT3: 

•console- 

99 

•aerial  porta- 

device 

baud 

parity 

#data 

fa  top 

rata 

bits 

bits 

300 

B3| 1 

300 

N 

6 

1 

301 

1 

300 

0 

7 

1 

302 

1  ■ 

300 

E 

7 

1 

1200 

1 

1200 

N 

8 

1 

1201 

1200 

0 

7 

1202 

1200 

E 

7 

1 

4000 

CONI: 

4800 

N 

8 

1 

4801 

CONI: 

4800 

0 

7 

1 

4802 

4800 

E 

7 

1 

■  ' 

N 

8 

1 

9601 

1 

9600 

0 

7 

1 

9602 

9600 

E 

7 

1 

n 

parity: 

NaNone 

E*Even 

O-Odd 

C0H2:*Add  SO  to  value  for  C0M1: 
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Tibi*  2*.  Modal  Nuabar 


Modal 

Nuabar 

Pr intar - P lottar -Screen 

Davica  Idantification 

0 

Epson  FX-80  Printer,  single  density. 

1 

Epson  FX-80  Printer,  double  density. 

2 

Epson  FX-80  Printer,  double  speed,  dual 
density. 

3 

Epson  FX-80  Printer,  quad  density 

4 

Epson  FX-80  Printer,  CRT  Graphics  I. 

5 

Epson  FX-80  Printer,  plotter  graphics. 

6 

Epson  FX-80  Printer,  CRT  Graphics  II. 

10 

Epson  FX-100  Printer,  single  density. 

11 

Epson  FX-100  Printer,  double  density. 

12 

Epson  FX-100,  double  speed,  dual  density. 

13 

Epson  FX-100  Printer,  quad  density. 

14 

Epson  FX-100  Printer,  CRT  Graphics  I. 

15 

Epson  FX-100  Printer,  plotter  graphics. 

16 

Epson  FX-100  Printer,  CRT  Graphics  II. 

20 

HP  7470A  Graphics  Plotter. 

30 

HP  7475A  Graphics  Plotter. 

40 

Epson  LQ-1500  Printer,  single  density. 

41 

Epson  LQ-1500  Printer,  double  density. 

42 

Epson  LQ-1500,  double  spaed,  dnal  density. 

43 

Epson  LQ-1500  Printer,  quad  density. 

51 

Houston  Instrument  DMP-51  HP  or 

DMP-52  MP  Plotter,  0.001"  step  size. 

52 

Houston  Instrument  DMP-51  MP  or 

DMP-52  MP  Plotter,  .005"  step  size. 

60 

HP  2686A  LaserJet  Printer  or  LaserJet 

PLUS  printer,  using  A  size  paper 
(8. 5"  x  11")  (216  am  x  280  urn). 

Drawing  resolution:  75  dots  per  inch. 

61 

HP  2686A  LaserJet  Printer,  using  B5  size 
paper  (7.2"  x  10.1")  (182  um  x  257  am). 

Drawing  resolution:  75  dots  per  inch. 

62 

HP  2686A  LaserJet  Printer,  using  A  size 
paper  (8.5"  x  11")  (216  am  x  280  am). 

Drawing  resolution:  150  dots  per  inch. 

63 

HP  2686A  LaserJet  Printer,  using  B5  size 
paper  (7.2"  x  10. l")  (182  am  x  257  an). 

Drawing  resolution:  150  dots  per  inch. 

64 

HP  2686A  LaserJet  Printer,  using  A  size 
paper  (8.5"  x  11")  (216  am  x  280  an). 

Drawing  resolution:  300  dots  per  inch. 

65 

HP  2686A  LaserJet  Printer,  using  B5  size 
paper  (7.2"  x  10. l")  (182  an  x  257  am). 

Drawing  resolution:  300  dots  per  inch. 

cont inued 
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Table  2a.  (continued) 


80 

IP  75808,  HP  75851,  or  HP  75668  Drafting 

Plottor  ualag  also  A/A4  to  D/Al  papar. 

HP  7S50A  Graphics  Plotter  using  slsa 

A/A4  to  8/83  paper. 

HP  7440A  ColorPro  plotter  using  size 

US/A4  paper. 

85 

HP  75858  or  HP  75668  Drafting  Plotter 
using  sise  E/AO  paper. 

90 

Tektronix  4025. 

99 

IBM  color  graphics  aonitor  (CRT). 

Tabla  2b.  Dot  Matrix  Printer  Usage  by  Modal 


Printer 

Model 

0 

1 

2 

0 

D 

5 

6 

10 

a 

12 

13 

14 

15 

iij 

Epson  FX-80 

* 

* 

* 

* 

D 

D 

G 

■ 

■ 

B 

fl 

n 

fl 

fl 

Epson  MX -60 

* 

* 

* 

* 

■ 

■ 

m 

H 

I 

m 

m 

■ 

IBM  Printer 

* 

* 

* 

• 

■ 

■ 

■ 

m 

m 

®g 

■ 

n 

m 

S: 

Centronics  GXP 

* 

* 

• 

• 

1 

■ 

■ 

m 

1 

m 

1 

m 

m 

f§ 

Okidata  92 

# 

* 

• 

* 

I 

■ 

■ 

I 

m 

fil 

1 

§g 

m 

i 

Epson  RX-80 

* 

* 

* 

• 

□ 

■ 

□ 

1 

I 

Epson  FX- 100 

■ 

■ 

■ 

M 

El 

D 

D 

3 

B 

B 

Epson  MX* 100 

1 

■ 

M 

El 

O 

I 

IK 

m 

Okidata  93 

■ 

■ 

■ 

B 

fl 

fl 

B 

■ 

1 

*  ■  Tha  printar  can  uaa  thla  nodal  nunbar. 


* 
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Table  3.  Racoa»aadad  Configuration 


Output  davica 

Device 

Model 

Epson  FX-80 

0 

5 

Epson  MX -80 

0 

1 

IBM  Pr intar 

0 

1 

Centronics  CLP 

0 

1 

Okldata  92 

0 

1 

Epson  RX-80 

0 

1 

Epson  FX-100 

0 

15 

Epson  MX- 100 

0 

11 

LQ-1500 

0 

41 

Okidata  93 

0 

11 

HI  DMP-51 

9600/9650 

51 

HI  DMP-52 

9600/9650 

51 

HP  7440A 

9600/9650 

80 

HP  7470A 

9600/9650 

20 

HP  7475A 

9600/9650 

30 

HP  7550A 

9600/9650 

80 

HP  7580B 

9600/9650 

80 

HP  7585B 

9600/9650 

80/85 

HP  7586B 

9600/9650 

80/85 

HP  2686A 

9600/9650 

60/61 

Tektronix  4025 

4800/4850 

90 

IBM  color  graphics 
aonitor 

,99 

99 
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RECORD  SELECTION 


From  the  opening  screen  repeated  here  the  user  has  two  choices. 
Records  may  be  selected  directly  by  identification  number  and 
title,  see  Appendix  A,  or  the  user  may  choose  to  select  records 
which  most  closely  match  a  set  of  parameters  he  establishes. 


]  =  Record  Natch  Tine  History  Response  Spectra  Options  Exit 

Select  Records 
Optimize 

USE 


ARROU  KEYS  THEN  PRESS  ENTER 
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Let  us  choose  SELECT  RECORDS  to  directly  select  records  fron  a 
list.  After  a  brief  delay  in  which  the  record  list  is  read  by  the 
conputer  the  user  is  given  a  display  of  record  numbers  and 
titles.  Use  the  DP  and  DOWN  ARROW  keys,  the  PAGE  UP  and  PAGE  DOWN 
keys  and  the  HOME  and  END  keys  to  find  a  record  then  press  ENTER. 
Repeat  selecting  one  or  more  records  to  a  limit  of  up  to  10 
records.  Upon  completion  of  the  selection  press  ESCAPE. 


Choose  up  to  10  files.  Press  [Esc  1  when  done . 


A001  EL  CENTRO  SITE  IMPERIAL  <M 
A001  EL  CENTRO  VALLEY  IRRIGATIO 
A002  FERHDALE  CUT  HALL  CON? 
A002  FERHDALE  CIII  HALL  COUP 
A00Z  FERHDALE  CITY  HALL  COHP 
A003  PASADENA  -  CALTECH  AIHENAE 
A003  PASADENA  -  CALTECH  ATHENAE 
A003  PASADENA  -  CALTECH  ATHENAE 
A004  TAFT  LINCOLN  SCHOOL  TUNNEL 
AB04  TAFT  LINCOLN  SCHOOL  TUNNEL 
A004  TAFT  LINCOLN  SCHOOL  TUNNEL 
A00S  SANTA  BARBARA  COURTHOUSE 
A005  SANTA  BARBARA  COURTHOUSE 
A005  SANTA  BARBARA  COURT  HOUSE 
A006  HOLLYWOOD  STORAGE  BASEMENT 
A0B6  HOLLYWOOD  STORAGE  BASEMENT 
A00S  HOLLYWOOD  STORAGE  BASEMENT 
A007  HOLLYWOOD  STORAGE  P.E.  LOT 
A007  HOLLYWOOD  STORAGE  P.E.  LOT 


*1 


You  will  then  be  shown  you  choices  and  asked  to  confirm  the 
selection;  press  MY"  for  yes  or  "NH  for  no. 


yOU  HAVE  SELECTED  THE  FOLLOWING  RECORDS: 
142  C848  8244  ORION  BLOB.  1ST  FLOOR 


IS  THIS  CORRECT?  V N 


If  the  second  choice,  OPTIMIZE,  was  selected  from  the  opening 
menu,  the  user  would  be  given  the  following  questions: 


Enter  Earthquake  Magnitude  6.5 

exanple  6.5 

Enter  Epicentral  Distance  in  NILES  18 

exanple  18.8 

Enter  Site  Acceleration  Level  in  G's  8.25 

exanple  8.25 

Include  ALLUVIUM  SITES  (Ties  or  (N)o  T 

Include  INTERMEDIATE  SITES  (Vies  or  (N)o  N 

Include  ROCK  SITES  (Vies  or  (Nlo  N 

Include  UNCLASSIFIED  SITES  (Vies  or  (Nlo  N 

Any  Changes?  V  /  N 


ENTER  REQUESTED  VALUES 
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Enter  the  EARTHQUAKE  MAGNITUDE ,  SOURCE  to  SITE  DISTANCE  in  miles, 
and  the  SITE  ACCELERATION  in  g*s.  Next  indicate  by  typing  "Ynto 
include  or  "N"  to  omit  records  from  ALLUVIUM  SITES,  INTERMEDIATE 
SITES,  ROCK  SITES  and  those  sites  which  are  UNCLASSIFIED.  Upon 
completion  type  "N"  for  no  changes  or  "Y"  to  repeat  the 
questions.  At  this  point  the  computer  will  use  the  entered  data 
to  compare  the  records  in  the  data  base  and  order  them  in  terms 
of  the  closest  matching  records  to  the  prescribed  parameters. 
Equal  weight  is  given  to  the  three  paramaters  of  magnitude, 
distance  and  acceleration.  The  computer  will  give  a  list  of  the 
records  in  order  with  the  closest  matching  record  first.  Select 
from  the  one  or  more  records  from  the  list  up  to  a  limit  of  10 
records  by  pressing  ENTER.  Upon  completion  press  ESCAPE.  The 
user  will  then  be  shown  his  choices  and  asked  to  confirm  them  by 
pressing  HY"  or  may  re-select  by  pressing  "N" 


Choose  up  to  10  files.  Press  [Esc!  when  done. 


TXMB  HISTORIES 


BEFORE  TIME  HISTORIES  CAN  BE  PROCESSED,  SPECIFIC  RECORDS  MUST  BE 
SELECTED.  SEE  THE  PREVIOUS  SECTION  FOR  RECORD  SELECTION 
From  the  main  menu,  selecting  TIME  HISTORY  produces  the  following 
sub-menu : 


5  Record  Natch 

Tine  History 

Response  Spectra  Options 

Exit 

% 

Scale 

Screen  Plot 

< 

Hard  Plot 

Urite  File 

USE 


A8R0U  KEYS  THEN  PRESS  ENTER 


The  user  may  choose  to  plot  or  write  the  time  history  data  as  is 
or  to  SCALE  it.  Scaling  is  accomplished  by  multiplying  the 
acceleration  record  amplitude  by  the  ratio  of 

Value  Entered  by  IZseg 
Peak  Acceleration  of  the  Record 

If  SCALE  is  selected  from  the  menu  the  following  question  will  be 
asked: 


Enter  Peek  Acceleration  in  9'* 

1 


Enter  the  value  of  the  peak  acceleration  to  which  the  record  is 
to  be  scaled  per  the  above  ratio. 

The  user  may  obtain  screen  plots  of  the  time  histories  of  each  of 
the  earthquake  records  previously  selected  or  may  opt  for  hard 
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copy  plots.  NOTE  the  plotter  MOST  BE  CONFIGURED  PRIOR  TO 
SELECTION  OF  HARD  PLOT,  SEE  SECTION  ON  CONFIGURATION.  See 
Appendix  B  for  an  example.  A  last  option  is  to  write  the  time 
history  data  to  a  file  whose  name  is  specified  by  the  user.  The 
data  from  the  selected  records  is  then  written  in  ASCII  format  to 
that  file.  This  file  may  be  used  as  the  load  file  for  a  finite 
element  structural  program.  Note  the  increment  between  data 
points  is  0.02  seconds  and  the  values  are  written  in  g's. 

RESPONSE  SPECTRA 

BEFORE  RESPONSE  SPECTRA  CAN  BE  PROCESSED,  SPECIFIC  RECORDS  MUST 
BE  SELECTED.  SEE  PREVIOUS  SECTION  FOR  RECORD  SELECTION 
From  the  opening  menu  select  RESPONSE  SPECTRA  and  the  following 
sub-menu  will  appear. 


Record  Natch  Tine  History  Response  Spectra  Options  Exit 


Sale 

Screen  Plot 
Herd  Plot 
Hrite  File 


USE  ARROll  KEYS 


THEN  PRESS  ENTER 


As  with  time  histories,  the  earthquake  records  selected  can  be 
scaled  using  the  SCALE  option  as  discussed  above.  Both  screen 
plots  and  hard  copy  plots  can  be  obtained.  Screen  plots  include 
plots  of  acceleration  versus  period  for  five  damping  values. 
Envelope  value  plots  and  plots  of  mean  and  mean  plus  one  standard 
deviation  are  given.  Hard  plots  of  the  same  are  given  plus 
tripartite  plots  are  also  given.  See  Appendix  B  for  an  example. 
The  data  may  also  be  written  to  a  file  whose  name  can  be 
specified  by  the  used. 
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DISTRIBUTION  LIST 


NAVFACENGCOM  /  CO,  ALEXANDRIA,  VA 
NAVFACENGCOM  CHESDIV  /  CO,  WASHINGTON,  DC 
NAVFACENGCOM  LANTDIV  /  CO,  NORFOLK,  VA 
NAVFACENGCOM  NORTHDIV  /  CO,  LESTER,  PA 
NAVFACENGCOM  PACDIV  /  CO,  PEARL  HARBOR,  HI 
NAVFACENGCOM  SOUTHDIV  /  CO,  CHARLESTON,  SC 
NAVFACENGCOM  SOUTHWESTTOIV  /  CO,  SAN  DIEGO,  CA 
NAVFACENGCOM  WESTDIV  /  CO,  SAN  BRUNO,  CA 
PWC  /  PENSACOLA, 

PWC  /  GREAT  LAKES,  IL 

PWC  /  CO,  PEARL  HARBOR,  HI 

PWC  /  CO,  FPO  AP 

PWC  /  CO,  SAN  DIEGO,  CA 

PWC  /  CO,  OAKLAND,  CA 

PWC  /  CO,  FPO  AP 

PWC  /  PHILIPPINES,  FPO  AP 

PWC  /  YOKOSUKA,  FPO  AP 


